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ABSTRACT 


, A comprehensive literature review covering the period 
from 1964 through midyear 1970 on heat .pipe technology is 
presented. A brief citation of early heat pipe work is 
followed by a presentation of heat pipe phenomenology in 
which the mechanism of operation, external boundary condi- 
tions, operational limits, the influence of noncondensable 
gases, and startup behavior, are discussed. Experimental 
investigations directed at determining the suitability of 
various substances for use as working fluids and wicks are 
described. In addition, numerous experimental studies 
dealing with operational characteristics of heat pipes are 
evaluated. Several possible areas of heat pipe application 
including heat transfer, temperature control, heat flux 
conversion and control of thermal conductance are examined. 

A discussion of basic heat pipe theory together with numer- 
ous modifications and simplifications concludes the review. 

An experimental program is described which consisted of 
three phases. The first phase involved the design and con- 
struction of an apparatus for the measurement of wicking 
properties such as the maximum wicking height and the 







transient velocity characteristics in a porous sample. ' The 
second phase of the experimental program involved the design, 
construction, and testing of a conventional cylindrical lieat 
pipe. Wall temperatures were measured along the pipe as a 



* * 


function of both heat input and p,ij»c orientation. It was 


found that the temperatures at the evaporator end were 
extremely sensitive to orientation. 

The third and final phase of the investigation con- 



sisted of experimentation with a coplanar heat pipe. The 
objective of this part of the study was to determine to 


what degree the components in a two component heat pip*’ 
separate. Both single component and water-methanol mixtures 
were used for the working fluid. The temperature distribu- 


tion was measured under rhe wicks and in the vapor space 
for various combinations of power input, orientation, and 
methanol mass fraction. It was found that stratification 
had an influence on the temperature field within the vapor. 
The temperature profiles measured with water-methanol mix- 
tures exhibited a consistent trend when compared to the 
profiles obtained with pure water. The behavior of the 



profiles suggests that a partial separation of components 
occurred. The phenomenon of complete component separation 




CHAPTER I: INTRODUCTION 

In rccont years it has become increasingly important 
to develop methods for the efficient transport of thermal 
energy from one location to another. The advent of the 
space age has moreover stimulated research for heat transfer 
devices which are light weight and have relatively long life 
expectancies. Such a device, although not for space appli- 
cation, was first proposed by Gaugler [1] of the General 
Motors Corporation in 1944. Unfortunately for Gaugler, the 
thermal transport problems of that time could be solved 
using more conventional heat transfer methods and devices, 
thus effectively concealing the true potential of his inven- 
tion for some 20 years. In 1962 Trcfethcn [23 submitted a 
report to tho General Electric Company in which he suggested 
the possible use of a passive thermal device for spacecraft 
applications. This device was to consist of a hollow tube 
with a porous 1 iner covering the inside surface. Energy 
would be transferred from one end to the other by means of 
?. capillary induced, continuous mass cycling. Ko. experi- 
mental verification of this concept was attempted, however, 
and the suggestion was quietly buried in company files. 

In 1964 Grover, et al. [33 of tho Los Alamos Scientific 
Laboratory independently rediscovered a device similar. to 
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Caugler and Trefethen’s and coined the name "heat pipe" to 
' describe it. Grover and his co-workers Were working in the 
area of spacecraft power generation at the time hut they 
immediately recognized the potential of the heat pipe in 
other areas. A heat pipe is defined as a closed structure 
containing some working fluid which transfers thermal energy 
from one part of the structure to another part by means of 
vaporization of a liquid, transport and condensation of the 
vapor, and the subsequent return of the condensate from the 
condenser by capillary action to the evaporator. Because 

energy is transferred by the flow of a pure saturated vapor, 

■ ' * * * • - . *«. r ’ • ■ ' - 

a heat pipe is usually very nearly isothermal. However, 

heat pipes often operate nonisothermally if the saturated 
vapor is "contaminated" with a second gas. The second gas 
may either be introduced into the pipe intentionally for the 
purpose of control, or may be accidentally admitted through 
leaks or improper filling procedures. If the second gas 
(vapor) participates in the evaporation processes, then the. 
pipe is called a "two component" heat pipe. If the second 
gas is unable to liquify and participate in the heat trans- 
fer, it is usually referred to as a noncondensable gas. - 
Many investigators have had experience with noncondcnsahlc 
gases in heat pipes (usually unintentionally) but little 
has been done to . inves tigatc the operation of two or multi- 
component heat pipes. This lack of activity has provided 
the motivation for the present research program. ' ’ ' ■ ■. 
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v Problem Formulation 
The first paper which mentioned two and multicomponent 
heat pipes was published by Cotter in 1965 ( t ) . Itc per- 
formed an elementary one dimensional analysis and concluded 
that for a heat pipe operating at steady state, the various 
components would separate and form a series of segments each 
containing a pure component which would operate as $ii 
‘'independent" heat pipe within the mul ti -component heat 

* ’ 1 • * „ ' , % . y ' 

pipe. The temperature of each segment would adjust itself 
so that the pressure within the vapor is nearly’ con^akhtt 
throughout the entire pipe, thus forming a series of tern-, 
perature plateaus. The most volatile component, of course, 
would be located at the condenser end of the heat pipe. 

S. Katzoff [5] analyzed the wall temperature distribution 
for a copper heat pipe containing a mixture of ethanol and 
methanol. He assumed that the energy needed for the leg 
with the more volatile (methanol) component is transferred 
past the ethanol-methanol interface primarily by conduction 
along the wall' and wick structure. If this heat flow is 
sufficient to maintain circulation in the methanol segment, 
then the situation described by Cotter would exist, i.e., 
the two legs would operate independently and the temperature 
distribution in the pipe would comprise two plateaus. In - 
the transition region between the two legs the composition;.', 
of both the liquid and vapor would, vary continuously, but r ’ 


. i 
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at each location the relative compositions would have to be 
such that the liquid and vapor arc in equilibrium./ Katzoff 
speculated that if on the other hand the conduction past the 
interface is not enough to maintain the methanol leg, then 
it may be possible for an ethanol-methanol mixture to exist 
along the entire pipe in both the liquid and the vapor. 

This would be possible if tho composition varied both 
radially and axially instead of just axially. ; 

Tien 1 9 2 1 measured the ax*al temperature distribution 
along the outside of a water-ethanol heat pipe. He deter- 
mined the pressure inside of the pipe and compared - this; with 
the pressure which had to prevail if pure ethanol occupied ; 
the condenser and pure water occupied the evaporator. He ._ 
concluded that separation into pure components in a heat 
pipe is extremely difficult, if not impossible, to achieve. 
Instead he found that if the initial composition was rich in. 
ethanol, the data attested to the existence of a water- " 
ethanol mixture in the evaporator while nearly pure ethanol 
(i.e.., the azeotropic, mixture) occupied the condenser see- ' 

* * . ■ ■ - ■ y • 

tion. All of Tien’s, data, however, were .obtained wir.h the 
pipe operating vertically with the evaporator below the .1 
condenser. If the pipe contains excess liquid such an t. • 

arrangement is usually - referred to as a reflux condenser * 
because gravity forces, instead of capillary forces can ■ • • 

always return the condensate to the evaporator. In ^uch 
cases results would have to be viewed with some caution " ' 


rhC^vri 






because the wick structure may or may not significantly 
alter the liquid vapor equilibrium conditions. 


' It is clear that present knowledge about even the funda- 
mentals of multi -component heat pipe operation is rather 
limited. It was the objective of this study to extend this 
knowledge for the case of two component heat pipes. The 




aims of the research reported in this study were: 

1) to' perform a survey and evaluation of the heat 
pipe literature 

2) to experimentally study and describe the basic 

physical behavior of a two component heat pipe, 
and •' ' -v ' 1 1 

3) to analytically predict the qualitative tempera- 

ture and concentration profiles in the vapor core 
of a two component heat pipe. ' 

Because an understanding of a single component heat 
pipe and some of the problems associated with it is a pre- 
requisite for the investigation of a two component heat - 
pipe, chapter II contains a discussion of the general 
phenomenology of heat pipes. In chapter III a review of. 
all literature pertaining to heat pipes is presented. The 
articles to be discussed include those published up to - 
approximately March 1970, In addition to these, about 




twelve articles to be presented at the 1970. ASME Space - 
Technology and Heat Transfer Conference , Los Angeles , ' 

California, June 21-24 will be referenced in an appendix to 
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this stud/. Due to the large number of references cited, . 
chapter III is somewhat lengthy. Readers primarily interested 
in the experimental investigation may find it expedient to 
proceed directly to chapter IV.* The analytical formulation 
for a two component heat pipe is discussed in chapter IV. 

In chaptor V the various apparatus used during the course of 
the current research are described in detail. Experimental 
results are presented in chapter VI together with the pre- 
dictions. from chapter IV. A number of appendices in which 
, ' . i ; 

supporting information is contained concludes the study. 


Chaptors II and III together with a modified intro- 
ductory chapter are being published as a monograph entitled, 
"The Heat Pipe," in volume seven (7) of the "Advances in 
Heat and Mass Transfer," New York, Academic Press, 1971. 




Descriptio n and Types oT Heat Pipes 

As evident from the definition given in the introduc- 

» ‘ 

tion, all heat pipes have a number of common features. 
First, all heat pipes incorporate what is usually referred 



tn as cn evaporator. This is the part of the heat pipe 
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through which thermal energy from some external source is 
introduced into its walls and from there subsequently trans- _ - 
ferred to the working fluicT. Second, all heat pipes include 
a condenser section. The working fluid condenses here and ' 
ultimately transfers its heat of condensation to an external 
sink. Many heat pipes contain also an adiabatic section 
located between evaporator and condenser.. The adia.atic 
section, besides providing a passage for the fluid, serves 
no functional purpose other than separating the heat source 
and heat sink to make the heat pipe compatible with any • 
given external geometric requirement. , V 

In addition to the longitudinal sections, i.e,, 
evaporator, condenser, and adiabatic section, a heat pipe . ~ 

may also be subdivided for the purpose of discussion into . •' ’. ’•••*• 



three radial components. The outermost shell is usually 
referred to simply as the "container." The container's /'in- 
sole- mechanical purpose is to enclose the functioning parts 
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>f the heat' pine and to lend it structural rigidity. Since 


the internal pressure is often different from the environ- 
mental pressure, the container must be capable of with- 


standing pressure differences without bulging or bursting. 
This constraint, along with cost and manufacturing considera- 


tions has led to the wide use of cylindrical "pipes’* as con- 


taining structures. In addition to fluid and ’’pressure". ' 
containment, the container also acts as an important part of 
the heat flow path from the source to the sink (sec Figure 
2.1). Hence the container wails should be thin to minimize 


their therma-i resistance. This feature is in direct opposi- 


tion to the thick wall requirement for .pressure containment 


and hence an opportunity for an optimization presents itself. 


The next radial clement is usually referred to as the wick, 


For ease of the present discussion, this may be regarded 
simply as a porous material filled with small random inter- 


connected capillary channels. Various types of wicks and 


their properties arc discussed in greater detail later in 
this study. The wick returns the liquid from the condenser 
to the cvapo»3tor utilizing the surface tension forces of 
the liquid. Although it is not a requirement, the wick is 
usually firmly attached to or pressed against the inside . _ 


wall of the container. Since the wick is in general 


saturated with a low conductivity working fluid (except in 


the case of liquid metals) the wick-fluid matrix represents^ 


usually the major resistance along the heat flow path.-; - It 
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Figure 2.1 Heat Pipe Components 


is therefore necessary to consider thermal properties as 
well as liquid transport properties when selecting a suita- 
ble wick. The interior space of the heat pipe is usually 
referred to as the vapor core which provides a passage for 
the vapor as it flows from the evaporator to the condenser. 

It should now be evident that the heat pipe definition 
entails no geometric constraints in regard to its structure; 
and in fact a large number of heat pipes of many different 
shapes have been built and tested. Several conventional and 
unconventional heat pipes are depicted in Figures 2.2 
through 2.4. A . heat pipe , as originally conceived by 
Grover, et al. is illustrated in Figure 2.2-A [3). This 
particular geometry exhibits two features which early 
investigators felt were important for efficient . heat pipe 









operation, i.e., a relatively large length to diameter ratio, 

and a porous wick material which covers the inside surface 

of the structure. Figures 2.2-B and 2.2-C illustrate two 

typical heat pipe configurations which also have a large 

length to diameter ratio, but which provide for capillary 

transport of the liquid in grooves and crevices forming an 

integral part of the containing structure, contrary to the 

porous wicking material sketched in Figure 2. 2-A, which is 

only held against the inside wall. The heat pipe shown -rn 

Figure 2.2-D has a very small length to diameter ratipu 
% ' i ■ • . : 

Heat pipes having such proportions are often called “vapor 

chambers" or "vapor chamber fins." The device shown*© 

Figure 2.2-E also fits our definition of a heat pipe although 

the liquid and vapor flow paths are separated mechanically 

'i 

whereas in the mor^ conventional heat pipe only the liquid- 
vapor interface separates the flow’ paths. Finally, the . 
device depicted in Figure 2.2-F has recently been introduced 
[ 6 ] as a "rotating heat pipe." Here the liquid return is 
caused by centrifugal forces of the rotating . contrivance. 
Although the rotating device appears promising fdTmany 
thermal transport problems, it does not fit our definition 
of a heat pipe and hence will not be discussed in detail. . 

Evidently the variations in heat pipe shapes are 7 ’ 
unlimited; both Katzoff [SI and Conway and Kelley [7] have ' • 
considered a doughnut shaped heat pipe. Several investi- 
gators [S, 8 , 9] have proposed and designed flexible heat 
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pipes. RCA has built and operated heat pipes with a variety 
of geometries. Among these are the two illustrated’- in 
Figures 2.3-A and 2.3-B (10]. Tho configuration displayed 
in 2.3-A will effectively transport heat around a 90° bend. 
The five-pronged device in Figure 2.3-B allows the use of 
any combination of prongs as evaporators and of the remain- 
ing prongs as condensers. For case of manufacturing j both 
types have cylindrical cross sections although this is not 
a requirement. In addition to the above geometries, 
several investigators (10, 11] have built and tested a so- 
called radial heat pipe. As illustrated in Figure 2.4, the 
radial heat pipe provides for thermal energy transport from 
a heat source to a concentric heat sink. The wick lining 
the inner walls of the annulus in this case is complimented 
by spokes consisting of additional wick materials. Here, 
as with most other heat pipe geometries, the relative posi- 
tions of the condenser and evaporator may be interchanged 
in order to accommodate any particular thermal transport 
problem. The variety of geometries depicted above arc by 
no means inclusive of all possible configurations' and are 
presented only to illustrate the extreme versatility of the 
heat pipe for heat transfer problems. 

Functioning of the Heat Pipe 
At first glance, the operation of- a heat pipco.appcars • 
exceedingly simple. Thermal energy is transferred from the 
evaporator to the condenser by continuous mass cycling and 
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phase change of a suitable working fluid.. The mechanise) of 
phase change with the accompanying absorption or release of 
the latent heat of transformation has long been recognized 
as an efficient heat transfer process. Many gadgets, c.g., 
the coffee percolator and the reflux condenser, combine heat 
transfer by phase change with a gravity induced mass cycling. 
Boilers in many cases utilize a mechanical pump to continu*. 

t * 

I . ^ V • • ’ 

ously circulate and replenish the working fluid. In a heat 
pipe, however, the working fluid is continuously cycled by 
the surface tension fo.rccs of the fluid itself. It is this 
unique method of mass transfer which has both stimulated a 
growing interest in the Heat pipe, and has also proved to he 
one of the major impediments for a successful heat pipe 
operation. To better understand the functioning and the 
limitations of heat pipes let us consider in more detail 
the physical effects occurring in a heat pipe. . 

The steady state operation of a heat pipe may be repre- 
sented schematically as shown in l-'igure 2.S. The inside 
wall of the container is lined with a porous capillary 
structure which is saturated with some working fluid. A 
sufficient amount of fluid must be supplied in the container 
in order to fill (saturate) all the pores of the capillary 
structure. The penalty for having a slight excess of fluid 
is small compared to the possibility of heat pipe failure 
which might arise from a deficiency of fluid. 'The vapor in 


the core of the pipe is essentially f.t the saturation 
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Figure 2.5 Heat Pipe Schematic 


pressure corresponding to the liquid surface temperature. 

In actuality the saturation pressure of a vapor in equi- 
librium with a liquid surface depends also on the radius of 
curvature of the surface. The vapor pressure is greater 
than that acting on a plane surface. if the liquid surface 
is convex, and less if the meniscus is concave. The effect 
is usually too small to warrant consideration and is not 
significant until the meniscus radius is of the order of 
one micron [12]. Since the typical capillary pores in most 
heat pipes arc larger than one micron, no noticeable error 
is introduced by neglecting this effect. : 




The heat transfer from the source to the sink is . * 

:/ .affected mainly by six simultaneous and interdependent 
processes: 1) heat transfer from the source through the 

container wall and wick- liquid matrix to the liquid-vapor 
interface; 2) evaporation of the liquid at the liquid-vapor 
interface in the evaporator; 3) transport of the vapor in 
the core from the evaporator to the condenser; 4) condensa- 
tion of the vapor on the liquid-vapor interface in the 
condenser; 5) heat transfer from the liquid-vanor interface 
through the wick-liquid matrix and container wall to the 
' sink; and 6) return flow of the condensate from the con- 
denser to the evaporator caused by capillary action in the 
wick. Let us now consider each of these processes separate!, 
and in more detail. 

Heat transfer from the source to the liquid-vapor inter 
face in the evaporator is essentially a conduction process. 
For low conductivity fluids, e.g., water or alcohols, the 
thermal energy is conducted through the. wick-liquid matrix 
almost entirely by the porous wick material since the wick 
has a higher thermal conductivity than the fluid. For high 
conductivity liquid metals, however, the heat is conducted 
both through the wick structure and by the liquid in the 
pores. Heat transfer by convection is very small because , 
the pores arc too small for any significant convection cur- 
rents to develop. The temperature drop associated with con- 
duction across the wick-liquid matrix depends on the working 
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fluid,- wick materials, wick thickness and the net radial 
heat flux. . This temperature drop may range from a few 
tenths to several hundred degrees Fahrenheit and is one of 
the major. temperature gradients along the heat path. 

Once the thermal energy has been transferred to the 
vicinity of the liquid-vapor interface, evaporation of the 
liquid can take place. As the liquid evaporates, the net 
mass flow away from the surface causes the liquid-v«y>or 
interface to recede into the Kicking structure. The result- 
ing concave shape of the meniscus, shown in Figure '2. S, is 
responsible for the functioning of the heat pipe. A simple 
force balance on a single pore shows that for a spherical 
interface the pressure of the vapor, exceeds the liquid 


pressure by an amount equal to twice the surface tension 
divided by the meniscus radius. This pressure difference 
is the basic driving force for both the liquid and vapor 
flows. It is opposed mainly by the gravitational and viscous 
forces acting on the liquid during circulation. The assumed 
form of the liquid-vapor interface sketched in Figure 2.5 is 
probably quite realistic for relatively low heat fluxes. As 
the heat flux' increases, however, the meniscus recedes even 
further into the wick and assumes a more complex shape [13] 
which may eventually interfere with the liquid flow in the 
capillaries. Once the liquid has absorbed the latent heat ‘ 






flow is caused by a small pressure difference prevailing in 
the vapor. core. This pressure difference is caused by the 
slightly higher temperature (saturation pressure) in the 
evaporator as compared to the temperature (and hence lower 
saturation pressure) in the condenser. This temperature . .. 
drop is often used as a criterion for successful heat pipe 
opc‘.-j‘on, and i f the difference is less than 1 or 2°F, the 
heat pipe is often said to be operating in the "heat pipe 
regime" i .e. ,• isothermally (3, 14, 15]. As the vapor flows. i; 
toward the condenser, additional mass is added from the down- 
stream portions of, the evaporator and consequently the mass 
flow rate and velocity in the axial direction continue to - 
increase throughout the cvapcratoi . Inverse conditions pre- 
vail in the condenser section of the heat pipe. 

The vapor flow in the evaporator and condenser of a 
heat pipe is dynamically identical to pipe flow with infec- 
tion or suction respectively through a porous wall ... The w 
flow may be either laminar or turbulent depending on the 
operating conditions of the heat pipe. As the vapor flows • 
through the evaporator (and the adiabatic section) the pres- " 
sure continues to decrease due to both viscous and accelera- 
tion effects. Once the condenser section is reached and the 
vapor begins to condense on the liquid-wick surface, a par- 
tial dynamic recovery in. the decelerating flow tends to in- 
crease the pressure in the direction of fluid motion;.! 1 1: 

V. on Id be mentioned that the driving pressure in the vapor 
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core is somewhat smaller than the vapor pressure difference 
of the fluid in the evaporator ar.d condenscrl This is so, 
because the vapor pressure of the liquid in the evaporator 
must exceed the pressure in the adjacent vapor in order to 
maintain a continued evaporation process. Likewise, the 
pressure of the condensing vapor must exceed the vapor 
pressure of the adjacent liquid in order to maintain con- 
tinued condensation. •. - ■ 

As the vapor condenses the liquid saturates the pores 
in the condenser. The meniscus has a very large radius of 
curvature, and in fact it may be considered essentially 
infinite. Any excess working fluid in the pipe collects on 
the condenser surface thus virtually insuring a plane inter- 
faced The heat of condensation is conducted through the 
wick-liquid matrix and container wall to the heat sink. 

If excess liquid is present, the temperature drop from the 
interface to the outside of the container will be larger 
than the corresponding temperature drop in the evaporator. 

In fact, some investigators [14, 1 6 J feel that the thermal 
resistance in the condenser is one of the major parameters 

- -v* ‘ ' ‘ * • 

to be considered in heat pipe design. 'v ’ 

Finally, the condensate is "pumped” through the wick to 
the evaporator by capillary action. The liquid flow is gen- 
erally regarded to be laminar and assumed to be dominated: by 
viscous forces. The pressure along the liquid flow path •- 
decreases due to both the viscous losses and the increase- in 
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elevation if the heat pipe is' operated in a'gravity field. • 
Operation of the heat pipe in presence of gravity with the 
condenser above the evaporator actually defeats the purpose 
of the wick since gravity can be used to return the conden- 
sate along the inside of the container Wall with less "• 
Viscous loss than liquid flow in the wick would cause. .In 
this mode of operation the heat pipe is said to have degen- 
erated into a reflux condenser or therisosyphon.' Therefore, 
in this monograph, the pressure loss due to gravity will 
always be considered greater than zero (evaporator above 
condenser) or equal to zero (horizontal orientation of heat 
pipe simulating a gravity free environment). " 

Because the vapor temperature, or operating temperature 
as it is sometimes called, of the heat pipe is essentially 
determined by the coupling of the heat pipe to the heat 
source and the heat sink, a brief discussion of possible 
source-sink combinations and their effect on the heat pipe 
operation is warranted. The vapor temperature adjusts 
itself in such a way that the temperature drop across the • 
wick-liquid matrix and the container wall in evaporator and 
condenser is adequate to transfer the given heat flow from 
heat source to heat sink. In other words, the absolute 
vapor temperature is established in response to the tempera- 
tures imposed on both the evaporator; and condenser by the ; 
source and sink. . The temperatures at the "outside wall of K’; 
the heat pipe may be either, ."fixed" or '’floating” depending r, 
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■ ’ the typo of constraints imposed by sere* or »*; . | 

' , r „ "flooring" temperature is usually the.,' : 

the" evaporator, f hcatflux boundary condition . ' ; 

result of - This u easily accomplished by ■ ■ : . 

upon the heat sourc . * R p co ils or 

• . incsnters induction coils, W. con 

employing resistance * " At the condenser, 

■ . . r he heat source. At tnv 

radiative heating , j y ‘ radiative 

a "Coating- temperature is -only t V . 

V fixed tenocrature can be «-‘»»tain . 

ric the heat pipe .th constant t-rature haths or 
uti'iting the heat o£ evaporation or ccndensa .. . ■ 

lie^uorhihg fluid lor heat addition or removal ^ _ 

lv at constant temperatures. - 

• consider possible qualitative temperature V 

Let us now consider. p iV for - - 

: " heat flow path (see Figure 2.1) 

profiles along _ tio ns Figure 2. 6-A depicts. the . 
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Figure 2.6 Qualitative Temperature Profiles 


Along Meat Flow Path 
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of source and sink is commonly found in laboratory testing 
of . low temperature heat pipes (resistance heating and water 
jacket cooling). As shown by profiles a and' c, the soiree 
temperature and the vapor temperature increase with increas- 
ing heat flux. Conversely, if the heat flux is maintained 
constant and the sink temperature is increased as in pro- 
files a and b, the vapor temperature and the source tempera- 
ture again increase, but now the temperature gradients in 
the evaporator and condenser remain the. same. Figure 2.6-C 
describes the profiles which are obtained if the source • 
temperature is fixed and the sink temperature is allowed to 
float. Here we see from profiles a and b that the vapor.; 
temperature must drop in order to accommodate a larger heat 
flux for a fixed source temperature and conversely it must 
rise for a given heat flux if the source temperature is 
increased. A significant omission from Figure 2.6 is the 
case where both the source and the sink temperature are 
allowed to float, a situation often encountered during the 
laboratory testing of high temperature liquid metal heat 
pipes (induction or RF coil heating and radiative cooling). 
For this case the operating temperature of the heat pipe 
adjusts itself to a value at which the total neat input . 
equals the total heat rejection. Since this self-adjustment 
depends on the exact type of source and sink used, and 



possibly on certain properties of the container wall itselfy ' •; 
e.g. , electrical resistivity or emissivity, an exact/, ... • : ' J 




statement about the operating temperature cannot be made. 

In general,' increased heat fluxes cause an increase in the 
operating temperature of the heat pipe: r 

The proper functioning of heat pipes depends also upon 
a continuous circulation of the working fluid; consequently, 
it is not surprising that virtually all of the limitations 
(limits) for successful heat pipe operation are associated 
in one way or another with the interruption of this mass 
circulation; The limit discussed and analyzed most often .. *•••' 
in the literature is the so-called "wicking limit.” This ' , “ 

condition is reached when a given heat flux causes the 
liquid in the liquid-wick matrix to evaporate faster than ~ 

it can be supplied by capillary pumping in the wick (wicking 
action). Once this condition occurs, the liquid-vapor 
meniscus continues to recede into the wick until all of the . 
liquid has been depicted. The wick in the evaporator 
becomes dry and the container temperature increases without 
bound until a "burnout" condition is reached which usually 
results in destruction of the pipe. Of course, the burnout 
condition can only be reached if the source is of the float- 
ing temperature type. For a fixed temperature source, the, 
heat pipe would simply. cease functioning once the capillary 
limit is achieved, and no mechanical damage would occur. 

Two additional limits of heat pipe operation are common- ’ 
ly associated with heat pipe "startup" and low operating . ” ■ ' 

temperature conditions. One of these is called "sonic limit". 




a condition found in heat pipes in which the source tempera 
ture is kept constant while the sink temperature is lowered. 
The vapor density decreases and the vapor velocity increases 
correspondingly until the velocity becomes sonic. • The v**P or 
flow "chokes” at the evaporator exit. Just as it does when' 
the sonic condition is reached at the throat of a convergent 
nozzle [17). Once choking occurs, a further decrease of the 
sink temperature, in analogy to a reduction of the exit 
pressure in a nozzle, docs not; result any longer in an 
increase of the total heat f low. • - ' ' . 

A so-called "entrainment limit" is reached when the 
vapor velocity is high enough and the vapor stream shears 
off droplets from the liquid interface entraining and carry- 
ing them to the condenser. Quite frequently the droplets 
can be heard as they impinge upon the end cap of the heat 
pipe [IS]. The premature depletion of the working fluid 
from tho wick means that less liquid can reach the evaporator 
where it is needed for successful heat pipe operation. The 
entrainment limit depends to a large extent on the surface 
pore size of the wick material and also on the surface ten- 
sion of the working fluid. The use of small pore sizes and 
fluids having large surface tensions is perhaps the most 
effective way of avoiding liquid entrainmeirt. 

The relative position of the operating limits in a heat 
flow Q, versus temperature T plot, is illustrated in Figure 
2.7. Mere Q represents the total axial heat transfer rate 
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Figure 2.7 Limits to Heat Pipe Operation [18] 


and T is the average vapor temperature in the heat pipe- 
Successful heat pipe operation is possible only under con- 
ditions existing below the curve ABODE. The shape of the 
area under the curve may vary drastically depending on the 
wick material and working fluid used: however, the basic 
shape of each limit curve should remain as shown. Numerical 
examples for these various limits will be presented later in 
the text of this thesis. 

The presence of a noncondensable gas in a heat pipe may 
have a detrimental effect on heat pipe performance. The non 
condensable gas can be added intentionally for the purpose 
of control, or it can be the result of improper filling pro- 
cedures, container leaks, or chemical reactions between the 
working fluid and the container or wick material. 'Neglect- 
ine the control aspect for the present, the roost common 
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noncondcnsablcs ‘are air (from leaks) and hydrogen (from . > 
chemical reactions). - During heat pipe operation the. non-" 
condensable gas is swept to the condenser and forms a stag- 
nant gas layer.. The temperature in this zone adjusts itself 
in such a way that the total pressure in the vapor core 
remains approximately constant throughout. Heat is trans-. 

.* V ■ I- • . - - ' ’ _ •*'_ 

ferred through this zone to the liquid-wick surface primar- 
ily by conduction. Because this mode of heat transfer is 
extremely slow compared to that taking place in a normal 
condensation process, the zone containing the stagnant gas 
is effectively eliminated as a functioning part of the heat 
pipe. The result is an effective shortening of the pipe, 
thus reducing its total axial heat transfer capability. The 
length of the noncondensablc gas zone depends on the operat- 
ing temperature and pressure in the system. For increased 
axial heat fluxes and the corresponding increased pressures, 
the zone will contract and allow more of the condenser to 
become operative again. Conversely, for decreased operating 
pressures, the gas zone expands and reduces the area avail- 
able for condensation. If the quantity of noncondensable., 
gas or the operating pressures are such that the entire. con- 
denser section comes to lie within the stagnant gas zone, 
heat pipe operation ceases. The possibility -.of using a ' 
noncondensable gas for the variation of the effective con- 
denser area as a control technique will be expanded upon ;.v 
later. " • ‘ . 1 • r ". "A ■ 
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The> preceding discussion has dealt exclusively with 
steady state heat pipe behavior. Of equal importance for 
the practical use of these devices is an understanding of 
the transient operating conditions through which a heat pipe 
passes during startup. Cotter [19J conceived three basi?c 
modes of startup which may be recognized by the shape of the 
developing temperature profiles. The three modes are experi- 
enced when the evaporator is boated uniformly over its 
entire length at constant heat flux which may, however, be 
varied with time. The condenser is cooled uniformly cither 
by radiation or heat conduction to a sink hept at uniform 
temperature. The various modes of startup are illustrated 
in Figure 2.8, the abscissa of which represents the distance 
along the heat pipe axis and the ordinate the vapor tempera- 
ture. The uniform startup in Figure 2.8-A.takcs place when 
the vapor density is high .at the ambient temperature so that 
the working fluid begins to reflux throughout the pipe immedi- 
ately in response to an increase in the heat flux. This type 
of startup procedure may be accomplished very rapidly without 
detrimental effects to the pipe. The frontal startup in 
Figure 2.8-B is encountered when the vapor density is very 
low at ambient temperature, a case often observed when start- 
ing liquid metal heat pipes from room temperature. In this . 
case the vapor density is so low that the molecule mean free 
path exceeds, the vapor core diameter. As the heat flux is.- 
increased, the vapor density in the evaporator section rises- 
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A. UNIFOAM STARTUP 
(VAPOR DENSITY HIGH 
AT AMBIENT TEMPERATURE) 


B. FRONTAL STARTUP 

(VAPOR DENSITY LOW AT 
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C. FRONTAL STARTUP 

(NON CONDENSABLE GAS 
PRESENT) 


D. STARTUP FAILURE 
(RUNAWAY HOTSPOT IN 
EVAPORATOR) 


k.s : -€ x ^p * 


,.| .» .... .-s 


i . . ■ *.••. 

!" •" - 


I HEAT , PIPE 

EVAPORATOR < CONDENSER 




Figure 2.8 Transient Temperature Profiles of Heat Pipes 
During Various Startup Modes |'19] 
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and, the' molecule mean Free path becomes small compared to 
the vapor core diameter. The vapor in the evaporator sec- 


tion enters the continuum flow regime while the vapor -in the 
condenser remains in the free molecule flow regime with, of 
course, a transition region located in between. This mode 



} . 

of startup is further complicated by compressible flow 
effects since transonic va t .or velocities are achieved. 


Finally, the vapor may condense into liquid droplets in the 
vapor core since the vapor is nearly saturated in the 
evaporator but subcoolcd in its expanding flew toward the 
condenser. The frontal startup in Figure 2. 8-C illustrates 
a situation which can be expected if a significant amount of 
noncondcnsable gas is present, a case in which the evaporator 
heats up relatively uniformly. As the vapor temperature and 
hence pressure increase the noncondensablc gas. is moved 
toward the condenser where it collects in a fairly well 
defined zone. The temperature in this zone adjusts itself 



so that the total pressure in the vapor core is approximately 


constant. As the heat flux is increased and the vapor pres- 
sure and temperature increase, the noncondcnsable zone is 
further compressed thus causing the temperature profiles 
displayed in Figure ."..8-C. This mode of startup may also be 
accomplished very rapidly. The startup modes described .■ 
above arc somewhat idealized and various intermediate modes, 
may be observed depending on the vapor density at initial ” 
temperature and the amount of noncondcnsable gas present. 
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Figure 2.S-D illustrates a common failur*: encountered during 
startup. A hot spot is formed in the evaporator and the 
temperature increases without bound until failure results. 
The failure is usually a consequence of cither the attain- 
ment of the wicking or boiling limit during the startup 
sequence of heat pipes having a floating temperature source. 




’ CHAPTER III: LITERATURE SURVEY 


General Literatur e 

A very extensive research effort has been devoted to 
heat pipes since 1964 when Grover and his. co-workers (3] 
at the Los Alamos Scientific Laboratory, Los Alamos, Xcw 
Mexico, first reported the successful operation of a heat 
pipe. A 347 stainless steel container lined with five 
layers of 100-mesh 304 stainless steel screen saturated 
with 40 gm of sodium became the prototype of all subsequent 
heat pipes. Five chromel -alumel thermocouples were welded 
along the 9-cm-long pipe. The temperature distribution was 
measured for various input power levels ranging from SO to 
600 watts. Of particular interest among the experimental 
results reproduced in Figure 3.1 are the constant tempera- 
ture plateaus extending from the heated end of the pipe 
revealing the zones which were refluxing. The temperature 
drops occurring at the unheated end were attributed to : 
stagnant hydrogen gas formed by the impure sodium at elevated 
temperatures. The measured temperature gradient in the’ 
refluxing region amounted to less than .05°K/cm. If the 
heat pipe were considered a solid rod, it would have an 
effective thermal- conductivity in excess of 24,000 Btu/hr 
ft°F. It is not surprising that a device capable of such 
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Figure 3.1 Temperature Profile Xiong 

Heat Pipe, Grover ct al. £33 

performance stimulated considerable interest among large V 
numbers of researchers. Subsequently, research and develop 
ment programs were initiated simultaneously in many univer- 




sity and industrial laboratories leading to considerable 
duplication of the research efforts. Much of this effort 
was focussed on the determination of basic material proper- 
ties, especially Kicking properties. Because of the large 
number of papers published in a relatively short period of 
time, it is very difficult to present a chronological dis- 
cussion of the references. Instead, the literature will be 
dealt with by subject in the following order: material 

properties, operating characteristics of heat, pipes , heat 
pipe applications, heat pipe control and heat pipe theory. 
It should not. be overlooked that, several very useful review 


v., 


- ^Y- i 


i - 'Mji 








V? ' 


«,-v;v*4 




P ■ 1 












• -,_*v > > " 




; » Zx \\. ’• . »**«-.*»'*, 


articles [20, 21, 22] have appeared in the literature, and . 
while they did not add new experimental or theoretical 
information, they probably encouraged further research and 
hence contributed to the growing field of heat pipe tech- 
nology. * 

Material Tests 


..‘ it*-, V* 4 


- Working Fluids 

The choice of a working fluid for a heat pipe applica- 
tion is dictated to a large degree by sevoral physical pro- 
perties of the fluid and by the chemical compatibility cf 
the fluid with the container and the wick. Deverall and 
Kemme [23] were the first investigators to formulate the * 
requirements for suitable heat pipe fluids: 1) high latent 

heat of vaporisation, 2) high thermal conductivity, 3) . low 
viscosity, 4) high surface tension, 5) high wetting ability, 
and 6) suitable boiling point. Parker and Hanson [24] have 
pointed out that the vapor pressure curve dictates the tem- 
perature range of applicability for a given fluid. In gen- ... 
eral t a fluid should be used in a steeply sloped region of 
its vapor pressure--temperature curve so that the temperature 
change associated with the given pressure drop is minimized. 

In addition, the vapor pressure should be reasonably high,- 
since a low vapor pressure would result xu low vapor densities 
and high pressure drops in the vapor flow. ■_ • 
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as a function of temperature [8, 2C, 27, 28, 29,39); In 
addition, Frank et al. f 30 ] have published extensive property 
group plots for water, sodium, and cesium while references 
[51, 32) yield values of liquid metals for possible heat pipe 
application. Langston and kunr [13] have presented a table 
comparing the value of N for several low temperature fluids 
including the freons, alcohols, and glycols. Basiulis and 
Dixon [53] have assembled some property data on.potential 
w rking fluids which are electrically, insulating. 

The choice of a particular fluid, of course, depends on 
the specific application; however, a few .general conclusions 
can be drawn. For high temperature., high heat flux heat 
pipes, the liquid metals are definitely superior to non- 
metallic liquids due to their vapor pressure characteristics, 
high surface tension, and high latent heat. The outstanding 
flu»d for low temperature work appears to be water due to , 
primarily its high surface tension and latent heat. No . 
valid statements can be made at this time regarding the 
suitability of cryogenic fluids for low temperature heat 
pipe applications. • . 

Wicks 

Prior to the recent interest in heat pipes , the major- 
ity of work on flow through porous materials came from such 
diversified fields as soil mechanics, petroleum engineering, 
water purification and ceramic engineering. The particular 
type of flow in porous bodies under study wn:-. usually either 



a gravity induced flow or a forced or pressurized flow. 

Many publications in these fields as well as the fundamental 
theories on flow through porous media are presented and dis- 
cussed in a book by Schcidcggcr [34]. Instead of gravity or 
mechanical work, the heat pipe on the other hand, utilizes 
capillary induced fluid flow for its operation. This par- 
ticular feature, along with several other prerequisites to 
be discussed subsequently, has led to considerable research 
efforts aimed at developing or finding existing wick materi- 
als suitable for employment in heat pipes. 

As mentioned earlier, the primary requisite for a heat 
pipe wick is that it acts as an effective capillary pump. 
That is, the surface tension forces developed between the 
fluid and the wick structure must be sufficient to overcome 
all viscous and other pressure drops in the pipe and still 
maintain the required fluid circulation. Rccausc the heat 
pipe may often be required to operate in a gravity field 
with the evaporator located above the condenser, the wick 
should be capable of lifting the working fluid to heights 
equal to or greater than the maximum difference in elevation 
between the evaporator and condenser. The requirements are 
of opposing nature since on one hand large pore sizes are 
called for to minimize the viscous loss in the wick and on 
the other hand small pore sizes arc needed to provide for 
sufficient capillary pumping and maximum lift height. As 
a result, seme sort of pore size optimization procedure 



appears warranted and in fact numerous authors have addressed 
themselves to this problem and their work will be treated 
later. 

The property data on potential Kicking materials, which 
have been accumulated to date, originate primarily from Kick- 
ing height measurements and the measurement of the permea- 
bility, which is defined as a proportionality constant 
between the flow rate and the pressure drop in a porous body. 
The techniques used and the values obtained for these proper- 
ties will be discussed later. In addition to the operating 
-characteristics, several mechanical features must be con- 
sidered when examining potential heat pipe wicks. Of spec al 
importance is the reproducibility of a wick structure so 
future heat pipe investigators may rely on data generated 
during earlier investigations. The wick should be mechani- 
cally stable and should be rigid enough so that its flow 
properties do not change in response to w»ick sagging ox 
stretching. The ease of wick fabrication and the cost arc 
also important and it is conceivable that these considera- 
tions could someday be the major criteria for wick selection, 
if heat pipes are ever mass-produced. 

A wide variety of wicks have been success ful ly employed 
in heat pipes. As mentioned earlier, the first wick, and 
also probably the most widely used to date, consists of 
several layers of fine mesh screen. Various methods have V 
been used to guarantee mechanical contact between the screen 


and the container. Meal f 2 8 ) rolled the screen on a mandril 
and upon insertion into the pipe removed the mandril.. The 
screen was held against the heat pipe wall by its own resil- 
ience, but Neal found that the resilience varied from screen 
to screen and hence the performance of the pipe was not 
reproducible. Pcverall and Kcmme [23] forced a steel ball 
through a heat pipe, after the screen layers had been 
inserted, and apparently achieved . good contact between the 
wick and the wall. However, no attempt was made to check 

i. 

the reproducibility of the. wick structures* Keramc [17] has 
constructed rigid screen wicTcs using the following procedure. 
Several layers of stainless steel screen were wrapped around 
a copper tube. The structure was placed into another copper 
tube and drawn through a die to compress the screen layers; 
the copper wafe removed chemically. The screen tube was then 
heated to 1000°C in a vacuum oven to bond the structure. 
Finally the baked, rigid screen tube ’.as inserted into a 
heat pipe where the screen and the crescent annulus between 
screen and wall formed an effective wick structure with 
reproducible qualities. Numerous other techniques for 
screen wick preparation have either been used or suggested. 
McKinney [35] has employed a coiled spring to hold a screen 
wick firmly against the heat pipe wall. Katzoff [5] con- 
structed wicks in which.- a single layer of screen is metal- 


lically bonded to the wall. One construction technique 
consisted of electroplating a thin coating of indium or tin 








onto both, the screen and the wall, pressing the two firmly 
together and baking them in an oven (200.°C for indium and 
275°C for tin). Excellent results were also obtained by 
diffusion bonding stainless steel screens to stainless steel 
plates. The bonding was effected by pressing the screen 
against the plate with a pressure of 15 psi while baking 
them at 1100°C for 2 hours in a vacuum oven. 

In addition to using screens as wicks, several investi- 
gators have employed screens only as a retaining structure. 
Heat pipes constructed at North Carolina State University 
[36, 37, 38, 39, 40 , 413 have utilized wicks consisting of 
various types of beads packed in an annulus between a retain- 
ing screen and the heat pipe wall. Kicks of this type have 
been successfully constructed using beads of monel, glass, 
and stainless steel of various diameters. 

: Several different textile fabrics have also been em- 
ployed as wicks. llaskin [25] used a rayon cloth as wick 
for a nitrogen heat pipe. The cloth was held firmly against 
the heat pipe wall by sliding together two halves of a 
slotted, diagonally cut retaining tube. Shlosinger [42, 45] 
and Shlosinger, et al . [44] selected a quartz fiber cloth as 
a wick for their experiments with flexible heat pipes, in 
which the cloth was pressed against the w'all with springs. 
Attempts to bond the cloth to the heat pipe wall met only 
with limited success. It was found that ordinary rubber 
cement produced acceptable bonds for . operating 


temperatures up to S0°C if the rubber was allowed to cure 
and become "tacky" before the wick was pressed against the 
wall. Experiments with higher temperature silicon rubber 
failed because low molecular weight silicon compounds were 
released in the curing process and effectively waterproofed 
the wick. Good bonding was achieved by applying heat seal- 
able film materials, such as thernosett ing and thermoplastic 
sheets of the polyester resin and polyethylene- type. 

Other investigators experimented with commercially 
available porous metals as potential heat pipe wicks. Neal 
‘[28] constructed a heat' pipe using a sitttTred copper fiber 
wick. Unfortunately the wick was not bonded to the heat 
pipe wall and poor results were obtained. Porous metal 
wicks arc extremely difficult to machine and normal cutting 
techniques such as hand-sawing, shearing and grinding tend 
to close the surface pores along the cut. The use of a 
filler material which could be removed after cutting is not 
generally rjeonmended, since it is difficult to completely • 
remove the filler, which changes the wetting characteristics 
of the porous material. Langston and Kunz [13] have micro- 
scopically examined porous metals, which were cut by electro 
discharge machining . (EDM) and electrochemical (ECM) machin- 
ing techniques. They found that EDM tended to erode the 
surface whereas ECM cut the porous material very cleanly 
with a minimum of» pore distruction. A technique to avoid, 
the machining of porous metals completely was recently 


developed { 4 S ] . The technique entails the application of a 
mixture of particular matter, binder, and solvent to a 
surface. As the solvent evaporates the surface tension of 
the fluid draws the particles together, compacting them, 
yet leaving open pores. Of special significance is the fact 
that with this technique wicking structures can be applied 
to geometrically complex surfaces. Wick structures have 
been successfully fabricated from powders of AljOj., SiC* . 
aluminum, copper and nickel. 

With the exception of screen wicks, the wicking struc- 
tures which have received the most attention are tire so- 
called low resistance and composite wicks. Bohdanskv et al. 
[46] first suggested the possibility of using channels, cut 
into the interior surface, running axially the length of the 
tube. Bussc and his co-workers [47] constructed several beat 
pipes employing such an integrated type of wick and found 
that the structure was very stable and the pore size was 
easily controlled. Kemme [17] has advocated the use of 
"composite" wicks with a- fine pore size at the liquid-vapor 
interface to provide good capillary pumping and a larger 
pore size underneath for the return flow of -the liquid. He 
fabricated such a composite wick from several layers of 
screen of different mesh sizes with the finer screens ■ 
installed on the inside to provide capillary pumping forces 
ana the coarse screen located in the annulus between the 
fine screen and the wall to serve as flow passage-, another 



into a cylinder and inserted into a heat pipe. The tri- 
angular passages formed between successive corrugations 
served as low resistance fluid return paths. Calimbas and 
Mulctt [493 modified the basic screen wick by placing nickel 
ribbon spacers between the layers of screen to create a 
series of concentric annuli. McSwecney [503 made a wick by 
using a scries of 1/16" diameter rods wrapped with 10 mil 
wire to space them apart. The rods were held against the 
pipe wall by a coarse screen. Turner and Harbaugh [513 and 
Turner [523 attempted to construct a noncircular heat pipe 



such that the corners of the tube would supply the capillary 
pumping force. His "configuration pumped heat pipes" were 
not very successul due to structural deficiencies causing 
the pipes to bulge and distort under even minor pressure 
differences between the interior and the surroundings. 

A wick property of major importance is the maximum; 
height to which a wick lifts a given v;orking fluid. Two 
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reasons must be cited to underline its importance for heat 
pipe operation. First, the maximum lift height places a 
constraint on the dimensions of a heat pipe if it is to be 
operated in a gravity field with the evaporator located 
above the condenser. Second, the measurement of the maximum 
lift height represents an efficient method for the evaluation 
of the capillary pumping capability of a potential wick. The 
capillary pumping pressure must be greater than the sum of 
all viscous pressure losses and gravity losses if the pipe 
is to function at all. Once the maximum lift height has 
been determined, the minimum effective radius of. the capil- 
lary structure can be calculated with equation (.5-2) 

. 2o cos 6 


max 


R r 


(3-2) 
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Equation (3-2) is obtained from a simple static force 
balance performed on a meniscus located a distance, h, above 
a free surface in a cylindrical tube. After r . has been 
determined, the capillary pumping pressure' is obtained from 
equation (3-3) 
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which is derived from a simple force balance. The idea of 
designating a pore radius or diameter is of course an ideal- 
ised approach because in general capillary pores consist of 
irregular non-circular channels.-. This irregularity.- in', 
channel size has prompted investigators at North Carolina 
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State University [36, 37, 38, 39, 40, 41} to define both a 
"rising" and a "falling" capillary ccjui librium height. In 
general the height to which a liquid will rise in a wick 
(rising height) is less than the height to which the liquid 
will fall in the same wick after it has first been completely 
soaked (falling height). Katzoff r 5 J has indicated that the 
difference between the two levels is typically of the order 
of 2S per .cent and attributes this to wetting difficulties 
associated with the rising fluid. . f .' : y 

Several investigators have experimentally measured 
capillary equilibrium levels in various wicks. Ferrell and 
his co-workers [ 40 ] measured the falling equilibrium height, 
of water in packed beds consisting of stainless steel parti- 
cles (40 to 100 mesh) and glass beads (80 to 100 mesh). 
Further work [ 41 } indicated that the equilibrium height as 
a function of particle diameter could be* closely predicted 
by assuming that the beads were arranged in a cubic array. : 
Phillips and llindorman [S3] measured the maximum capillary 
pressure for a 200 mesh screen of stainless steel, bronze 
and nickel using water, methanol and benzene respectively. 

In addition, several metal "foams" and ''felts" were tested 
with these fluids. (The porosities of the samples ranged 
from 89 % to 96 %. ) Capillary pressures were determined by 
the standard technique of measuring the equilibrium height 
and the pressure necessary to force an air bubble through a 
saturated wick. Both methods gave comparable results and 
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the latter approach proved to he easier and less time con- 
suming. For the screens, capillary pressures ranged from 
.12 psi to .36 psi and for Che porous metals from .02 psi 
to .15 psi. K 3 tzoff [£] measured the lift capability of six 
screens and found the values of the minimum effective menis- 
cus radius ranging between .75 and .90 times the spacing of 
the wires in the screen. Ernst. (54) studied Katzoff’s data 
and concluded that the effective meniscus radius can be 
expressed by (d ♦ d,), where d is the mesh opening half 
width and d ( is the radius of the* wire in the mesh.. Langston 
and Kunz [13, 55] measured the equilibrium height olt water 
and Freon 113 in 23 wick samples. The samples consisted of 
three classes of porous materials: sintered metal screens/ 

sintered meta] powders, and sintered metal fibers (the 
porosity of the samples varied from 47.73. to 91.81). Equi- 
librium heights greater than 16" were found for several of 
the sintered .nickel fiber samples with water as test fluid. 

In a study not dircctiy concerned with heat pipes Ginwala 
ct al. 156) examined 17S potential wicking materials which'.'-., 
included cellular types, textile and synthetic fibers, fil- 
ter papers, inorganic fibers, porous ceramic and refractory 
products, porous and fibrous metals, etc. Equilibrium 
heights were measured for the fifteen most promising mate- - 
rials. Maximum rise heights were obtained with the. Silica 
Vitreous fibers and filter papers while acceptable heights 
were observed in Viscous Rayon. , 





A property of great: importance when selecting a wick 
is its permeability. Permeability (K) has the dimensions 
of length squared and is defined by Darcy's law given by 
the following relationship: 

- P £ g x sinaj (3-4 

The permeability is dependent upon the dimensions and the 
geometry of the passages in the Kicking material and can be 
determined experimentally by passing a liquid through a 
wicking material and measuring the pressure drop in the 
direction of the flow (see Figure 3.2). The pressure drop 
and the measured liquid mass flow rate along with the area 
normal to the flow and the fluid properties arc then used 
in equation (3-4) to evaluate the permeability. The flow 
through the porous body may be either forced or gravity 





induced .since static pressure variations due to gravity are 

• * . . i ' . 

accounted for in the second tern in brackets. Unfortunately, 
no correlations between permeability and more easily measured 
wick properties havo been found which apply to general wick 
configurations. Attempts to correlato permeability with 
porosity have been unsuccessful [34] since porosity in no 
way accounts for the coupling of one poro to another. ' F or 
wick materials whose geometry is easily identified, as for - 

example spheres in a cubic array, correlations can readily be 

\. 

found as shown by Carnesalc, ct al. [36] and Ferrell and 
Allcavitch [41] who measured the permeability of several beds 
cf packed spheres of different diameter. For comparable pack 
ing the porosity of all beds was 40% and a correlation 
between permeability and particle diameter was obtained. 

Numerous investigators have measured the permeability 
of more complex wicks for which no correlations exist. 
Cinwala, et al. [563 measured the flow rate of distilled 
water flowing through various types of felts, fibrous mate- 
rials, and cellular materials for three different pressure 
heads. Their experiments indicated that the flow rate 
under constant pressure head decreased with increasing time 
for all Kicking materials. This behavior was thought to be 
caused by an accumulation of gases ar.d microscopic particles 
in the wick initially dissolved and suspended in the liquid. 
Langston and Kunz [13,. 55] experimentally determined the 
permeabilities for a number of sintered metallic materials < 



fabricated from fcltedi fibers, powders, and screens. Special 
care was taken to degas the fluids before thoir use. The 
permeability was found to be independent of the nature of the 
fluid, time, flow rate and fluid temperature. A special 'est 
in which water was intentionally aerated before being passed 
through the sample was performed to ovaluate tho effect of 
dissolved gases in the fluid. The results showed that* with 
a high degree of air saturation, the permeability of the wicks 
decreased by about 18t in about SO hours of operation. Using 
both equilibrium height and permeability data* values for a 
capillary pumping parameter defined as the product of maximum 
lift height and permeability were then evaluated for all sam- 
ples. The magnitude of this pumping parameter is a direct^ 
measure of the efficiency with which a material night func- 
tion as a pipe wick. The data attested that sintered metallic 
fibers as a group make tho best heat pipe wicks while sin- 
tered powders were the next best and screens were the worst. 
Phillips £293 and Phillips, et al. £533 measured the permea- 
bility of sintered metal screens, fibers, and foams using 
forced flow, gravity flow, and condenser flow in an operating 
heat pipe. For the force flow test the permeability was found 
to decrease for increasing flow rates. This is contrary to 
Langston and Kunz's results which demonstrate that permea- 
bility is independent of the fluid flow rate. The disparity 
may perhaps be explained by the fact that for comparable 


samples Phillips ct al. used a much wider range of flow 
rates than Langston and Kune did (by a factor of almost 
four). Their data showed for a 96t porosity nickel 

foam wick a variation in permeability of approximately 11% 
over the entire range of flow rates. Gravity induced flows 
were used to measure the permeability of very thin wicks 
(such as one or two layers of screen or sintered metal 
samples less than .0S0 inches thick) as a function of menis- 
cus radius and flow rate. The meniscus radius was main- 
tained constant along a sample by adjusting the flow rate 
or angle of inclination in such a way that the viscous pres- 
sure drop was exactly countered by the increase in static , 
head, This procedure assured pressure constancy throughout 
the liquid. The meniscus radius was varied by changing the 
gas pressure on the vapor side of the vapor- liquid interface 
Tor thin wicks the data revealed the permcahi 1 i ty as a 
strong function of meniscus radius and for some samples 
(especially single layers of sciiecn) also as a strong 
function of the liquid flow rate. Permeability values 
measured by forced and gravity induced flow techniques were 
in general not in agreement. The gravity induced flow 
technique yielded permcahility values up to 2 1/2 times . 
those measured by the forced flow method. The disagreement 
was attributed to the manner in which the samples were *> 
mounted. Tor the forced flow* tests, the wiexs were com- 
piece iy enclosed and the measured permeability was based . 


on flow solely in the porous structure. In case of the 
gravity flow tests a fluid fillot formed along the edge of 
the wick which allowed some of the liquid to bypass the 
wick, leading to an erroneous and high value for the perme- 
ability. Several attempts were also nude to measure permc- 
i 

ability of a wick in an operating heat pipe. The pressure 
distribution in the wick was measured at five locations in 
the condenser and the fluid flow rate was determined from 
tuc product of the measured heat transfer rate and the 
known latent heat of vaporization of the liquid. Unfortun- 
ately the results were erratic and consequently their repro- 
ducibility for identical experiments was very poor. 

Farran and Starncr (S71 measured the capillary pressure 
and the permeability of a compressible wide of braided Si0 2 
fibers. They noticed that the capillary pressure (referring 
to Figure 3.2) is defined ns, ' 


AP c " l ’v, * p . 
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could be evaluated from fP,-P 2 ) if the lower end of the wick 
was submerged and the meniscus had an infinite radius of 
curvature thus equalizing the pressure in both phases, 

P 2 » P (it is assumed that P yj * P VJ ). Equation (3-4) 
transforms to 
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or Inserting the relation for the mass flov rate in terms of 
the fluid voloclty leads to 
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Two techniques were selected to evaluate the pressure differ- 
ence, iP c * and the permeability. The first technique .(dis- 
played schematically in Figure 3.3-A) consisted of the 
measurement of the steady state mass flow in the wick. The 
fluid was removed from the top of the wick by evaporation^ 

The second technique involved the measurement of the trans- 
ient rise of a liquid in a previously unsaturated wick (dis- 
played schematically in Figure 3.5-B). Inspection of equa- 
tions (3-6) and (3-7) then reveals that if m (or dx/dt) were 
plotted as a function of 1/x, the resulting plot (Figure 
3.3-C) should be a straight line and the permeability and the 
capillary pressuro could be found from the intercept and the 
slopo, respectively. Data obtained by both methods failed 
to yield the expected linear relationship. For near hori- 
zontal wicks (sino » 0) the capillary pumping pressure and 
the permeability turned out to be constant and the two 
methods yielded values which differed by about 10%. For 
sino f 0 both methods indicated that A? c was proportional 
to (x sina) and the reciprocal of tho permeability was 
approx imately- proportional to (x sino) 2 . This behavior is 
perhaps best explained in terms of tho existence of an 
optimum capillary radius corresponding to an optimum capil-' 
lary flow area through which a maximum mass flow rate is 
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Figure 3.3 Determination .of the Capillary Pumping Pressure 
and Permeability, Farran and Starncr [57] 
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pumped by capillary action. Farran and Starncr showed an 
inverse variation of this radius with x sino;-and further- 
more, they observed that if the optimum radius exists along 
the entire wick, AP C should be proportional to (x sina) and 
1/K to (x sina) 2 . They hypothesized that if the pore size 
distribution is large enough at any particular position, x, 
enough optimum sized pores would be available to dominate 
the flow. However, a determination of the pore size distri- 
bution was not made for the test samples thus leaving the 
hypothesis unchecked. Finally it is quite noteworthy that 
Feldman [8] has collected and present* 1 o t ab J. e summarizing . 
permeability values for a variety of wicks x.. eluding the 
sintered metals and some compressible materials. 


Compatibility of Components and Life Tests 
The choice of suitable materials for heat pipe construc- 
tion is dictated by a eompatability criterion of the differ- 
ent materials. Many of the problems associated with long 
term heat pipe operation arc a direct consequence of mate- . 
rial incompatibility which usually, manifests itself in chem- ' ; ; 
ical reactions. In general, improper selection of components 
results in a gradual appearance of noncondensablc gases. For 
high temperature liquid metal heat pipes, improper material 
selection furthermore accelerates corrosion and dissolution 
of the wick structure. . Grover and his associates [3] . were vy :- ’V 
the first to encounter and describe the generation .of non- -i- ^ 
condensable, gases . In their experiment , which- is' described ; ' 
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earlier, the temperature profile along the heat pipe dropped 


suddenly at the condenser (see Figure 3.1). This was attri- 


buted to a pocket of hydrogen gas produced from impure sodium 


in the reaction, Nall Na *. Ih .• Andcen, tt al. (16) tested 

■> 2 


a water-brass heat pipe and experienced severe problems with 


noncondensable gas. No attempt however, was made to deter- 


mine the source of the gas. - Schwartz (58) noticed the occur- 


rence of a noncondensable gas in several water-stainless 


heat pipes. Samples of the gas were withdrawn from one of 


the pipes and their composition analyzed with a mass spectro- 


meter. The results of the chemical analysis indicated that 


the noncondensable gas was composed of over 97$ hydrogen. 


Schwartz hypothesized that the hydrogen was formed as a result 


of a chemical reaction between the iron in the stainless 


steel and water. He suggested that the problem of noncon- 


densable gas generation could be avoided, either by. choosing 


a heat pipe whose metal components range below hydrogen in 


the electromotive series' or if the metals have an electro- 


chemical potential above hydrogen, by using a nonreacting 


working fluid other than water. To test this concept, an ' 


ammonia -stain less steel heat pipe was built and operated con- 


tinuously for 3 months with no measurable sign of • nonconden- 


sable gas formation. Conway and Kelley [7j were also 


troubled by noncondensable gases in a. water- stainless - 


steel hoat pipe. Although no tests were made to determine 


the origin of the gas, it was certainly most likely 




that hydrogen was produced in the sane chemical reaction 
which plagued Schwartz. Crovcr [59] has suggested a possible 

I 

solution to the problems caused by hydrogen formation in low 
and moderate temperature heat pipes, lie recommended to 
fabricate the condenser end cap out of palladium allowing 
hydrogen to diffuse to the outside while retaining the work- 
ing fluids. The prohibitive cost of palladium, however, 
limits its use to only the most amply funded experimental 
programs and as yet, the concept has not been tested. 

Dcvcrall and Komme (f>0] have reported the successful opera- 
tion of a water-stainless steel heat pipe for over 3000 hours 
without accumulation of noncondcnsabl-c g3S. The stainless 
steel tube and screen were first degreased in acetone and 
then bright-dipped to guarantee clean surfaces which were 
subsequently degassed at <>000°C in a high vacuum oven. 

Since the other investigators who encountered hydrogen 
formation in water-stainless steel heat pipes reported no 
extensive cleaning procedures; it may be speculated that 
for the particular case of water-stainless steel, hydrogen 
formation is more a function of the techniques used in 
processing the. materials than of the materials themselves. 
Jeffries and Zerkle [27] have commented on work done. by 
Lyons who tested several fluids in capsules of aluminum 
alloy 6061 at temperatures from 1 55°F to 322°F for durations 
in excess of 500 hours. Strong evidence, of corrosion was 
found with methanol and ethanol; n-Butane (1S5°F) and 



Monsanto Cp-34 (321°F) showed moderate corrosion; whereas 
no corrosion was evident with n-Pcntanc (310°F), Benzene 
(310°F), Heptane (320°F), Toluene (322°F), Ammonia (159°FJ, 
Freon-11 (.156°F), and Freon-113 (15S°F). In addition, a 
capsule of 321 stainless steel was tested with water at 
520°F, Definite signs of gas evolution were evident during 
the test, which proved in a qualitative mass s pec trome trie 
analysis to be hydrogen. , 

The most severe compatibility problems arc encountered 

i- ■ 

when' heat pipes arc operated at elevated tempo rat; a ">s v A 
considerable effort was expended by workers at Los. -?.l amos 
{15, 25, 48, 59, 61] to study this problem. A heat pipe was 
constructed selecting tantalum as container and wick mate- 
rial, and silver as the working fluid. The pipe was' 
operated at 1900°C for 100 hours. F.xarai nation of the sec- 
tioned pipe revealed that the wick in the condenser had 
almost disappeared while the wick in the evaporator was 
clogged with tantalum. This behavior was attributed to the - 
dissolution of the wick in the working fluid. In this case, 
small amounts of structural material were dissolved by the 
working fluid in the condenser section, transported to the 
evaporator section and deposited there as the working fluid 
evaporated. The amount of tantalum deposited indicated a 
solubility of tantalum in silver of approximately 10 ppm. * 
Dissolution and transport of the wick material also occurred 
in an indium- tungsten heat pipe which was operated at 1900°C 
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for 75 hours. Recent tests have attested to the potential 
of a silver- tungsten heat pipe for high temperature opera- 
tion. Such a pipe has been operated for 1000 hours with no 
noticeable deterioration in performance. Burse and his 
co-workcrs [47, 62 r 63) have conducted an extensive experi- 
mental program to determine material compatibility at the 
temperatures employed in thermionic converters, i.e., 1000°C 

for the collector and lbOO^C - 1800°C for the emitter. 

i 

Special care W 3 S taken to select materials for which no 
known i ntermetal 1 ic compounds exist. For many material 
combinations, they found that the dissolution of the con- 
tainer or wick led to eventual heat pipe failure due to 
clogging of the wick in the evaporator. Other material 
combinations produced ultimate heat pipe failure due to a 
weakening of the wall caused by intergranular corrosion and 
wall penetration. It was also noticed for a particular 
thermionic converter-heat pipe system that the lithium work- 
ing fluid diffused through the wall at rates sufficient to 
significantly lower the power output of the converter. 

Three systems, however, proved to be promising for operations 
longer than 1000 hours at 16C0°C: h'/Li, K/Pb, and SGS-Ta/Tl. 

In the 1000°C temperature range, Na and Cs have operated for 
1000 hours in a h'b-lZr container with no significant corro- 
sion. Workers at RCA (10, 64, 6S, 66, 67, 68] have also 
investigated material compatibi li-ty at thermionic tempera- 
tures. lixamination of a lithium - TZM alloy heat pipe .*■ • 


i 

i 




employing a molyl donum wire wick revealed about 10 ppm of ‘ 
the alloy in the lithium working fluid after 600 hours of 
continuous unuegraded operation. Another Li-TZM heat pipe 
has been operated successfully for 1 ,400 hours while a K-.N'i 
heat pipe has supposedly been working continuously for 
26,000 hours. At lower temperatures, a water-copper heat 
pipe has accumulated 8500 hours of operating time with no 
degradation’ in performance. Iirnst et al. [60, 70, 71] have 
reported the successful operation of a Li-TZM heat pipe for 
up to 5000 hours, joisc deposits were noted in the evaporator 
and it was speculated that it was probably titanium oxide; . 
however, no attempt was made to determine the composition of 


the substance. Johnson 


investigated the compatibility 


at thermionic temperatures of the working fluids Ag, Ba, Ca, 
In, Li, I’h , and T1 with containers made of Cb-lZr, Ta-lOW, 
and TIM. The fluids (metals) were placed into reflux cap- 
sules which were heated for times up to 1000 hours. The 
capsules were then sectioned and examined by x-ray diffrac- 
tion analysis. Results indicated that indium is not suit- 
able for long term heat pipe operation. Only minor inter- 
granular attack was observed when calcium was used as work- 
ing fluid, while other material combinations showed varying 
degrees of attack and corrosion. 

In summary, for low temperature heat pipes, care must ' 
be taken when selecting components to avoid any combination 
for which a possible chemical reaction exists which could 
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lead to the formation of noncondensable gases. Cleanliness-. ; 

and material treatment techniques also seem to govern to a ; 

* *: _ 

large extent the ultimate compatibility of heat pipe mate- 
rials. For high temperature heat pipes, due consideration 
must be given to the formation of possible intermetal 1 ic 
compounds, solubility of one metal in another, and diffusion 
effects. Obviously a great deal of systematic research is 
called for in high temperature applications since very little 
is known about the crucial properties of metals at these 
tcinperature-sw 


jeratine Characteristics of Heat Pipes 


General 


In addition to the experiments concerned with basic 
investigations of potential heat pipe materials, many inves- 
tigators performed experiments to determine the operating 
cnaracteristics of heat pipes. In early investigations often 
little more than the successful operation of a heat pipe was 
reported. Shortly after the work performed at Los Alamos 
had been published, Bainton [73] reported the successful 
operation of two sodium-stainless steel heat pipes. Tempera- 
ture uniformity over most of the length of the pipe was veri- 
fied by infra- red photography . Workers at- RCA [66] related 
the successful operation of a lithium heat pipe tor- over 
9000 hours claiming that, as a thermal energy transfer 
device, this heat pipe operated 1000. to 10,000 times more 
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efficiently than the thermal conduction process in an cqulv- 
alent rod of metals such as copper or silver. Hall [74] 
tested a lithium-TZM heat pipe and verified that fluxes 
sufficient in magnitude to operate a thermionic convertc'- 
could be obtained, fluxes on the order of 40 watts per 
square centimeter were achieved for both a lithium and a 
sodium heat pipe. Bowman and Crain ( 75 ] operated a water- 

copper heat pipe at near ambient temperatures. The tempera- 

! ’ 

turc profile along the axis of the pipe was measured to con- 
firm heat pipe operation. An electrically insulating fluid 
was employed by Rasiulis and Dixon [11] in a radial heat 
pipe which could successfully dissipate 900 watts regardless 
of its orientation with respect to gravity. Although no 
mention was made concerning the uniformity of the flux 
across the condenser surface (sec figure 2.4), it is most 
probable that the flux around the circumference of the con- 
denser was a function of its orientation in regard to gravity. 

Investigations of Heat Transfer Limits 
Several investigators have experimentally determined 
the Kicking Limit of heat pipes as a function of vapor 
temperature and/or geometric parameters. Bohdansky et al. 

[76] measured the maximum possible heat flow in a /sodium-' 
niobium heat pipe in the temperature interval from 500°C to 
800°C, The pipe was SO cm long and had a 2. cm inner diame- 
ter. The capillary system consisted of SS grooves of 
rectangular cross section with a:,»idth of .4 mm and a depth 


of .46 nun. The pipe was heated with an RF coil and cooled 
through a variable resistance helium thermal bridge with 
cooling water on the other end. In subsequent experiments 
the inclination of the pipe was varied in order to change 
the lift height, h, between evaporator and condenser. Upon 
each variation in orientation the power input was increased 
until a hot spot appeared at the far end of the evaporator 
indicating that the wick was no longer capable of supplying 
sufficient fluid to this part of the evaporator. The 
temperature variation was measured with thermocouples which 
were mounted at the outside wall of the pipe. Bohdansky 
et al. (76) plotted the heat flow rate versus operating 
temperature shown in Figure 3.4, where the detrimental 
effect of lift height on the maximum 'heat transfer cart be 
recognised. Occurrence of maximum heat flows for each 
experiment performed at a constant Lift height is attributed 
by Bohdansky et al. to the decrease of the surface ?cns ion 
with increasing temperature. They even further illustrated 
the effect caused by elevation (lift) by replotting the 
maximum heat transfer values as a function of height as 
shown in Figure 3.5. Neal (28J in a qualitatively similar 
representation for a water-stainless steel heat pipe fitted 
with four layers of 10S mesh screen as <i wick structure 
confirmed the results of Bohdansky et al. It should be 
. noticed that the shape of the curves in Figure 3.4 corre- 
sponds to the form of the qualitative wicking limit curve 
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Figure 3.4 Heat Flow in a Na-Nb Heat Pipe as a Function 
of Temperature and Elevation, Bohdansky et al. [76 



H (cm) 

Figure 3.5 Maximum Heat Flow at 850°C as a Function 
of Elevation, Bohdansky et al. [76] . 



displayed in Figure 2.7. The curves can be compared quali- 
tatively because the wall temperature and vapor temperature 
arc closely related for heat pipes operating under steady 
state conditions. 

Cosgrove ct al. ( 37 ) investigated two water-brass heat 
pipes in which the wick structures consisted of packed monel 
beads which were held in an annulus between a retaining 
screen and the wall. For a particular wick structure and 
pipe orientation, the maximum heat transfer was considered 
reached when a hot spot began to form in the evaporator 
which was concurrently detected with thermocouples installed 
in the packed beads.. The vapor temperature remained rela- 

‘ A 

tivclv constant in all cxj>eriments in which the primary 
variables were the pipe orientation and the diameter of the 
ncncl beads. Figure 5.6 illustrates the effect of particle 
diameter, and consequently pore size, on the maximum heat 
transfer as a function of inclination. For a given particle 
diameter, the maximum heat transfer decreases with increas- 
ing elevation once more reaffirming the results of Bohdansky 
and co-workers depleted in Figure 3.S. Notice, that for a 
selected pipe. inclination the maximum heat transfer increases 
with decreasing particle diameter. This effect is caused by 
the increased capillary pressure resulting from the smaller 
pore sizes. Cosgrove could not explain why the curve with 
t lie smallest particle diameter' intersected with the other 
curves. If the pore sizes were decreased indefinitely 


>. eo 

2 

a 


particle diameter ,t» 

O I.T20 . UO*») 

A 1.210 <IO**» 

O .020 (10** > 

V .650 (»*•) 

$ .546 00**1 


COS (3 

p ig ure 3.6 Maximus Heat Flow as a Function of Inclination 
and Particle Diameter, Cosgrove, et al. [37] 

however, eventually an optimum pore size would he obtained 
at which the viscous drag in the capillaries would become 
dominant and the maximum heat flow should decrease again. 

The results obtained by Bohdansky and Cosgrove and 
their co-workers suggest that the capillary limiting curve 
depicted in Figure 2.'7 i;. in reality a family of curves, 
each depending on the pipe orientation in the gravitational 
field. Moreover, * the entire family of curves depends on 
the geometry of the heat pipe, and in particular on the 
characteristic?? of the wick structure. 










Figure 3.7 Heat Pipe Dryout, McSweeney £ 50 1 

While Bohdansky et ai. and Cosgrove et al, installed 
thermocouples in the wall and wick to signal the formation 
of a hot spot, McSweeney [50] has found that the vapor 
temperature is a much more sensitive indicator of wick 
dryout. He experimented with a Na-Stainless Steel heat 
pipe and monitored both the vapor and the wall temperature 
as a function of increased power input. From his data 
shown in Figure 3.7 it is evident that the vapor temperature 
at dryout changes more rapidly in response to power increase 
than the wall temperature. The vapor temperature decrease 
under dryout condition is somewhat puzzling to say the. 
least. Unfortunately the location of the thermocouple probe 
in the vapor space was not mentioned. However, it appears 
most probable that the pr.obe was located in the proximity 
of the condenser. Possibly the termination of the heat flow 
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caused by the dryout of the wick resulted in considerable 
temperature variations within the now stagnant vapor core . . 
with rising temperatures in the evaporator and decreasing 
temperatures in the condenser which was still cooled. It is 
entirely possible that this type of behavior occurs only in 
heat pipes with a fixed temperature boundary condition 
imposed on the condenser. If a floating temperature sink 
were employed the entire vapor temperature would most likely 
increase once the dryout condition is reached. 

In related studies, Shlosingcr [43] found that for com- 
pressible wicks specifically, the manner in which the wick 
is retained against the heat pipe wall may significantly 
affect its wicking characteristics. For instance, if a 
helical spring is used to retain the wick, the fluid may 
have to travel by capillary action over a much,longcr spiral 
path from the condenser to the evaporator than otherwise 
would be necessary with a different type of retaining struc- 
ture. This effect, of course, would not be present for T.ore 
rigid wicking materials or for the commonly used axial slots. 
For instance, Bussc ct al. [65] found that heat pipes with 
axial grooves serving as a wick often formed a hot spot on 
the top side cf the evaporator wiicn the pipe was operated 
in the horizontal position. They attributed this behavior 
to the missing interconnections between the parallel grooves, 
thus preventing any cross. flow between grooves in the . • 

evaporator. In the condenser, however, the excess liquid 
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tended to reside at the lower side of the pipe thus making 
the lower grooves the preferred paths for liquid transport. 
Their explanation was verified experimentally by shielding 
the lower side of the condenser of a H'/Pb heat pipe. The 
shield- forced the vapor to condense on the upper side of 
the pipe and subsequently no hot spot was formed. In a heat 
pipe in which the grooves are interconnected by circumferen- 
tial grooves for instance, local overheating should not 
occur in the horizontal mode of operation. A few investi- 
gators have studied the problem of heat transfer and boiling 
in wicks. In a research project not directly related to 
heat pipes, Allinghan and McEntire [77] measured boiling 
film heat transfer coefficients on a horizontal copper tube 
which was surrounded by a ceramic wick and immersed in a 
pool of water. For lower heat fluxes they measured values 
for the boiling film heat transfer coefficient in excess of 
those established under similar condi tions in conventional 
pool boiling. The higher values were attributed to an 
increase in effective heat transfer surface area and also 
to an increase in active nuclcation sites. 

At higher heat fluxes, however, the trend reversed 
itself and the values of boiling film heat transfer coeffi- 
cients decreased for wick boiling below the corresponding • 

V3 lues of normal pool boiling. The reason for the decrease, 
appears to be two- fold. First- the very presence of the • • : 
wick prevents agitation of the liquid otherwise so common 
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in pool boiling. Second the vapor escaping through the pores . J . ; 

impedes the liquid counterflow directed tov;ards the heating 

surface. The data were reduced and an implicit function of 

the boiling heat transfer coefficient as a function of the 

radial Reynolds number, based on four times the hydraulic ; 


radius of the capillary pores is plotted in Figure 3.8. 

Likewise, Anand [78] and Anand et al. [79] tested a 
water-stainless steel heat pipe which had a 100-mcsh stain 
less steel screen wick. • The wall and vapor temperatures 
were recorded for different axial heat flows and the wick 
boiling heat transfer coefficient , h^., was calculated. 
Anand claims that the data can be correlated by 



St = .0051 Pr/ 6 Np* 2 Re* 1 * ^ * . (3-8) 

However, the data were plotted using 

St ■ .00051 Pr* * 6 Np * • 2 Re' 1 • " * (3-9) 

I 

Moreover, the ordinate on their graph was in error by a 
factor of ten. Allingham and McEntirc corrected An and' s 
graph and presented it in Figure 3.9 along with their own 
correlation. Both correlations . show the same trend and are 
in relatively good agreement considering that one represents 
wick boiling on the outside of a tube, and the other one 
wick boiling in the interior of a heat pipe. The validity 



of Anand ’ s results , however, may be somewhat questionable ' 
in view of the mistake made in their graphical presentation. 
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3.8 Wick Boiling Heat Transfer Correlation 
Allingham and McEntire (77] 
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Marto and Hosteller [80] studied the problem of wick ' 
boiling using a so-called- everted heat pipe. A sectional 
view of their heat pipe is given in Figure 3.10. The unique 
feature of this pipe is a vapor space, enclosed in an annulus 
between an interior tube and a confining envelope. The wick 
consisted of four layers of 100 mesh stainless steel screen 
attached to the outside of the inner pipe. Heat addition 
and removal were accomplished using a resistance heater and 

* . i v : - ‘ ‘ 

a tap water cooling system both installed within the inner 
tube.' The outer envelope was made of glass to facilitate 
visual observation of wick boiling. The results obtained 
with water as working fluid arc demonstrated in Figure 3.11, 
from which it becomes apparent, that lower superheats were 
required under boiling conditions in a wick thin in conven- 
tional pool boiling. As the radial heat flux was increased, 
dryout of the wick occurred at the same flux value whether 
or not boiling was observed. The authors concluded that 
wick boiling could exist in a heat pipe with no detrimental 
effect on its operation. The system pressure was also 
varied for both water and ethyl alcohol and it was found 
that for a given heat flux, the superheat decreased as. the 
absolute pressure increased. In addition, for a given 
superheat the heat flux in the case of water was almost one 
order of magnitude larger than that obtained with ethyl .v 
alcohol. This disparity in heat, fluxes. is expected in view 
of the higher surface tension and latent heat of water 
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Figure 3.11 Observed Superheat vs. Radial Heat Flux, 
Marto and Mosteller {80] 



compared to the same properties of the alcohols. In a study 
involving much larger superheats, Langston and Kunz (13} 
measured the heat flux through several wick samples as a 
function of the superheat. They found that for superheats 
of the order of 15°F and larger the flux through the sample • 
wicks became much smaller than for boiling on a flat plate, 
thus confirming the results reported by Allingham and 



Menu tire. The samples used by Langston and Kunz included.;: 
sintered nickel powders and sintered nickel screens. 

i . • ■ 

In another series of experiments on boiling in heat 
pipes, Moss and Keiley[>13 employed a neutron radiographic 
technique to measure the liquid content (i.c., liquid thick- 
ness) in the wick of the evaporator ir. a co-planar heat pipe 
The wick was made of sintered stainless steel screen (1/4" 
thick) and the working fluid was water. Measurements proved 
that only under conditions of :cro heat transfer did the 
wick in the evaporator remain completely saturated. As soon 
as heat was supplied to the evaporator the 1 iquid interface 
receded into the wick reducing the extent of saturation in 
the wick. In addition, the data demonstrated that under 
normal operating conditions the degree of saturation of the 
wick in the evaporator was inversely proportional to the » 
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Figure 3.12 Vapor Blanket Thickness vs. Heat Transfer Rate, 

Moss and Kelley [313 

attempt to describe the heat transfer characteristics of the 
partially saturated wick. In a conventional model it was 
assumed that vaporization takes place at the liquid vapor 
interface. In the other model, however, the formation of a 
vapor blanket thickness as a function of the heat transfer 
rate and contact angle is shown in Figure 3.12. For heat 
fluxes smaller than 15,000 Btu/hr ft* the second model more 
closely predicts the measured values than the more conven- 
tional first model. The results led to some allegations 
that the study did not pertain to heat pipe operation since 
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it is generally believed that vaporization should take place 
' * ’ - * •* ' ’ ' • 
at the liquid vapor interface during successful heat pipe 

operation. The authors countered these allegations by insist- 
ing that their measurements indicated beyond any doubt that 
the previously accepted idea of heat transfer through wicks, 
is in error. At the present time, however, there are more 
data available supporting the model of conduction through 
the wick, with vaporization taking place at the liquid vapor 
interface, than results sustaining the observations made by 
Moss and Kelly. For example, Ferrell and Ailcavitch {41J 
measured the heat flux through packed beads saturated with 
water. Their data for 30-40 mesh monel beads arc displayed 
in Figure 3.13. It is seen that the data fall very close to 
the line predicted by assuming pure conduction through the 
saturated bed. The curve obtained from the conduction model 
was calculated under the assumption that heat flows by con- 
duction through a thin liquid-bead layer in contact with 
the heating surface. The thickness of the liquid layer was ; . 
determined by the location of the minimum pore diameter in 
the head configuration. It is somewhat surprising that such 
good agreement was found between their theory and experi- 
ments since in the experiments the bed of beads was com- 
pletely flooded to a level well above the upper surface of 
the bed; hence, no such thin. liquid-bead layer existed in 
actuality. Recent work by the same authors has included : 
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heat transfer measurements in a similar apparatus, except 
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Figure 3.13 Experimental Results for Surface Covered with 
30-40 Mesh Monel Beads Ferrcl 1 and Alleavitch [41] 


that now the bed was not flooded and the liquid was drawn to 
the evaporator section by capillary action identical in 
manner to an operating heat pipe. The more recent data 
showed excellent agreement with the conduction mechanism 
discussed above. . • 

Similarly, Phillips [293 and Phillips and Ilinderroann 
[S33 measured the heat flux through a composite wick of 
nickel foam and stainless steel screen using water as a 
working fluid. The fluid was supplied through 3n artery 
and moved by capillary forces. A typical sample of their .. 
data is illustrated in Figure 3.14. It is again obvious 
that conduction was the mode of heat transfer for low values 
of AT. The wick exhibited a hysteresis effect after 



Figure 3.14 l-xper incnta 1 Results for Distilled Water and 
Nickel l-'oam, Phillips C 29 3 and Phillips tail llindermann 15 

nucleate boiling was first observed. In addition, the maxi- 
mum heat flux decreased with decreasing chamber pressure. 

The reduction was attributed to the significantly increased 
sice of the vapor bubbles formed during nucleate boiling at 
decreased pressure; hence causing a premature burnout due 
to vapor blockage in the wick. The blockage occurred 
because the vapor was forced to vent through the top of the 
wick by purposely sealing its sides. The same effect was 
observed in experiments performed by Costello and Rcdcker 
[ $ 2 ] . They concluded that proper venting of the vapor was 
necessary if the full capabilities of the capillary supply 


system were to be utilized. It is interesting to note that 
at the present time, no consideration is given to proper 
- vapor venting of the wick in heat pipe design. 

In addition to the capillary and boiling limits which 
restrict the heat transfer capability of all heat pipes, 
those transporting heat of the order of kilowatts instead 
of watts, arc often limited by the sonic and/or entrainment 
limits. , These limits are often encountered during startup 
procedures from near ambient conditions where the initial 
vapor pressure is very low and the resulting velocities in 
the vapor core are consequently very high. Kemme [17] 
investigated the sonic limit using several different liquid 
metals as working fluids. The heat pipes were heated by an 
induction coil and cooled through a gas gap with a water 
calorimeter. The use of different mixtures of argon arid 
helium in the gap allowed heat pipe temperature variations 
at a constant heat input, or heat input variations at a 
constant heat pipe temperature. Figure 3.15 illustrates 
data obtained during the startup of a sodium heat pipe. The 
dashed line indicates the sonic limiting curve based on the 
vapor temperature existing in the evaporator exit.. The heat 
flow was increased in discrete steps and the pipe was allowed 
to reach steady state before the heat flow and wall tempera- 
ture measurements iverc made. The evaporator exit temperature 
followed the sonic curve until it reached 560°C. For tern- . 
peratures lower than 560°C, the flow in the^ condenser section 
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Figure 3.15 Startup Behavior for Sodium Heat Pipe 
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consisted of continuum flow at the entrance and free mole- - 
cule flow at the far end of the condenser. . Hence, only" that 
part of the condenser in which continuum flow existed con- 
tributed significantly to the heat removal from the pipe. 

As the heat flux increased, eventually the entire condenser 
region was in the continuum flow regime and once this 
occurred, the heat removal area of the system remained * : 
fixed, so that a further increase in heat input now resulted 
in a larger temperature rise at the evaporator exit than was 
previously possible. The subsequent vapor density increase 
allowed the velocity at the evaporator exit to become sub- 
sonic. Further experiments indicated the occurrence of- 
supersonic velocities at the condenser entrance. Figure 3.16 
shows results obtained from steady state measurements at a 
constant heat input of 6.5 kw. The condenser temperature 
was adjusted by varying the concentration of the gas mixture 
confined in the gap described earlier. Curve A describes a 
condition in which the vapor velocity remained subsonic 
throughout the heat pipe. As the pressure in the condenser 
was decreased, the velocity became sonic at the evaporator 
exit (curve B). The existence of subsonic flow between 
curves A and B is evidenced by the changes in the condenser 
temperature and pressure which were transmitted to the 
evaporator where corresponding changes occurred. When the 
pressure in the condenser was further .decreased,- as- shown . 
by curves C and D, the evaporator pressure conditions ~ ~ 
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remained constant; : The vapor velocity did not .decrease. - . 

immediately upon entering, the condenser but continued to % 
expand and became finally supersonic followed by a more, 
abrupt pressure recovery "han was evidenced for subsonic 
flow. Figure 3.17 demonstrates tlie effect of different?, 
fluids on the sonic limit. Good agreement between ex-peri-*, 
mental results and theoretical predictions was achieved.* 

The sonic limiting curve (see Figure 2.7) is highly depend- 
ent on the working fluid and is also dependent on the pres- 
sure and temperature at which the heat pine is operating. 

In a similar study, Dzakowic, et *al. [83] Confirmed 
the resul ts obtained by Kcmme. A sodium- stai n less steel 
heat pipe with five layers of 60 mesh screen serving. as wick 
was employed to study the vapor velocity limit in heat pipe 
operation. Figure 3.18 is a plot of the axial temperature 
profiles obtained for two different heat inputs. It is 
interesting to see that with increasing heat flow the tran- 
sition from supersonic to subsonic flow at the condenser 
entrance took place at nearly the same heat flow rate and 
-temperature as measured by Kcmme (sec Figure 3.15). This 
emphasizes again the dependency of the sonic limit on the 
selection of the working fluid and operating conditions of 
the pipe and deemphas izes its dependence on a particular 
pipe geometry. Various startup tests were also conducted. 
The data are presented in Figure 3.19 along with the calcu- 
lated capillary (or AP) and sonic limits. Note the data of 
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Figure ,3.7 9 Heat Transfer Rate vs. Temperature at the 
Adiabatic Section, Dzakowic, et al.. [83] 


Dzakowic , ct al . fall to the left of the sonic limit curve 
while Kcmmc's data (Figure 3.17) fall consistently to the 
right of the curve. An explanation for this discrepancy 
may be derived from the fact that Kemnc plotted his data as 
a function of tbo maximum evaporator temperature, while 
Dzakowic et al. used^the temperature measured at the mid- 
point of the adiabatic section. It appears that perhaps 
the use of the evaporator exit temperature would move both 
sets of data closer to the predicted curve. Dzakowic et al. 
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attributed their observed temperature discrepancy to possi- 
ble supersaturation of the vapor arid to a lesser degree, to 
the uncertainty in values of sonic velocity and specific _ 
heat of the vapor. . . .. 'V 

The sonic limit is not the only factor which contributes 
to startup problems of heat pipes. Obviously the starting 
technique is also of major importance.’ Bussc, et al. [47] 
found it impossible to start a magnesium heat pipe unless 
the entire pipe was first preheated. Krnst, ct al. [69] -■ 

experienced that the method of orienting the heat pipe 
during startup had a significant affect on its final operat- 
ing condition. Figure 5.20 depicts the evaporator tempera- 
ture drop ns a function of difference in elevation between 
the evaporator and the condenser; For height differences 
of less than 4.4 inches, no effect depending on pipe orien- 
tation was evident. Once the pumping limit was reached, 
however, the temperature drop was greatly influenced by the 
method used to attain the elevation difference. The upper 
curve represents the situation where the pipe was first 
placed into the desired position and subsequently heated, 
while the lower curve represents the case where the pipe 
was first heated horizontally and then subjected to the • 
desired orientation. No satisfactory explanation was given 
for the difference in heat pipe behavior. ..-McSwccney [50] 
measured the temperature i n the vapor space . of a sodium heat V 
pipe which during the startup transient, period unexpectedly 
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Figure 3.21 Transient Temperature Profile Along Partially 
Frozen Heat Pipe, Shlosinger [43] 


Neal (28) and Shlosinger [43] studied the startup 
behavior of low temperature heat pipes with the working 
fluid initially frozen. The transient temperature profile 
along a water heat pipe in response to 50 watts of heat 
input is illustrated in Figure 3.21. Even such a moderate 
heat input of only 50 watts caused wick dryout and subse- 
quent overheating ir. the evaporator section before the 
entire pipe thawed out and could begin to operate in its' 
normal mode. With a heat input of 15 watts, however, the 
pipe thawed without wick dryout and normal operation began 
after approximately one hour. In an attempt to shorten the 
transient period, an auxiliary heat pipe which employed a 
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,, working fluid with a lower melting point was bonded to the 
primary heat pipe. The use of the auxiliary pipe greatly 
enhanced the thawing of the primary heat pipe without local 
overheating and the transient period of time was reduced by 
a factor one hair. Initially the pipe was frozen only in 
the condenser section and liquid was present in the rest of 
the pipe. Deverall et al. [84) investigated the startup 
behavior of a water' heat pipe which was initially frozen 
over its entire length. Ten watts of power were supplied 

I . • - ' • ■ • 

and the transient temperature profile, shown in Figure 3.22, 
was measured. From the shape of the developing temperature 
profiles, Deverall deduced the following sequence of events: 
First heat was transferred in association with vapor flew 
along the full length of the pipe. The vapor was formed 
by sublimation in the evaporator. Since the wick was still 
" frozen no liquid returned to the evaporator to replace the 
sublimated fluid. The lack of return fluid resulted in a 
rapid rise in evaporator temperature until it (sec stations 
1, 2, and 3} rose above the melting point. Water was then 
wicked into the heater area, and the temperature rise was 
temporarily arrested. Since most of the fluid was still 
frozen, not enough liquid was supplied to the heater.: Dry- „ 

out occurred and the temperature of the evaporator increased 
rapidly to above 100° C. When the entire. working fluid was 
: ■ finally melted, sufficient liquid was wicked into the - 
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Figure 3.22 Transient Temperature Profile Along Completely 
Frozen Heat Pipe, Deverall, et al. [84] 

evaporator and the temperature fell suddenly to 44°C. 

Normal operation began and the pipe became isothermal. 

Startup data were also obtained for a (nenfrozen) water 
heat pipe placed in earth orbit. Telemetry date indicated 
a frontal startup of the type illustrated earlier in Figure 
2.8-B. The startup procedure and steady state operation of 
the hc3t pipe in orbit were similar to those experienced in 
laboratory tests. 

Phenomena associated with the startup behavior of high 
performance heat pipes have been studied by Kemmc [8S, 86]. 
He found that it was virtually impossible to start several 




sodium heat pipes with wicks consisting of axial channels 
(grooves) running the length of the pipe. The difficulty 
was attributed to the attainment of the entrainment limit. 

As mentioned earlier, this limit is reached when a high 
velocity vapor stream shears liquid out of the wick and 
impedes a continuous supply of fluid to the evaporator. The 
entrainment limit is dependent primarily on the charac<.eris _ 
tics of the wick surface. Figure 3.23 shows test data 
obtained for three wick surface configurations. The maximum 
heat transfer was significantly improved hv covering the 
channels with a layer of screens . The experimental curve 
labeled 4 is a good example of the entrainment limit and its 
shape confirms the curve initially depicted in figure 2./. 

I n another experiment, the effect of "screen fit" was 
studied and is shown in Figure 5.24. The mechanical fit of 
the screen was very important in establishing the heat 
transfer limit of the heat pipe. The higher limit with the 
loose fitting screen was probably caused by the extra fluid 
paths available in the annulus between the screen and the 
outer wall. The curve calculated for the contacting screen 
in Figure 3,24 was based efi the capillary limit and not the 
entrainment limit. The capillary and entrainment limits are 
both dependent on the wick structure and it is impossible to 
theoretically or experimentally change one without also 
altering the other. This dependency associated with the 
fact that the entrainment limit, is gcncrully attained only 
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in a relatively narrow temperature range (between B and C in 

■ ■ " . • . < ' ■ . - • _ - 1 

Figure’ 2.7) has somewhat curtailed the investigation of this 

particular limit of heat transfer in heat pipes. 

’ Basic Studies 

The investigations discussed in the preceding section 
have dealt primarily with the limits to heat pipe operation 
and startup characteristics. A number of additional studies 
have yielded information contributing to the rapidly growing 
heat pipe technology. Dcverall, et al. [60, 84, 87, 88] 
operated a water-stainless steel heat pipe in an earth 
orbit and its performance was monitored bv telemetry at 
several tracking stations during fourteen revolutions. 

Results indicated that there was no degradation. of the heat 
pipe performance in a zero gravity field. In further tests 
[26, 89] a similar heat pipe was subjected to various sinu- . 
soidal and random vibrations to determine the influence of •. 
a vibratory environment on heat pipe performance. The 
experiments proved that vibration was not detrimental to 
heat pipe operation. On the contrary, vibration enhanced 
the wetting of the wick, forcing liquid into al] parts of ■ 
the wick structure,’ and thus actually improved heat pipe . 
performance. CaHmbas and Hulett [49} confirmed these 



results in a series of vibrational tests performed with a 
water- stainless steel, heat .pipe . ••• : ,y \ s -. 


/. Haskin [25). measured the total radial- temperature drop 
in the evaporator anJ condenser in a low temperature nitrogen 
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heat pipe. He found that for low heat loads (less than 36 
watts) the temperature drop in the condenser was larger than 
i?i the evaporator. The difference was assumed to be caused, 
by the extra liquid which accumulated in the condenser during 
operation. For larger heat loads, the temperature gradient 
in the evaporator became greater. This was attributed to 
partial drying of the wick and to the formation of a super- 
heated vapor film on the inner metal tube surface. Similar 
measurements were made by Schwartz [S8] with a water-stain- 
less steel heat pipe. For the same range of heat loads, but 
much higher temperatures obtained with different fluids he 
found that the temperature drop across the wall and wick in 
the evaporator was consistently higher than in the condenser. 
Hence, a higher condenser temperature drop is not the rule, 
but instead, the relation between condenser and evaporator 
temperature drops is dependent on the given geometry and the 
nature of working fluids. / • •"* . 

In yet another study, Rankcn and Kemrae [48 J measured 
the temperature variation. along the length of a lithium heat 
pipe operating at about 85Q°C as shown in "igure 3.25. The 
measured points have been fitted with a smooth curve which 
was corrected for the temperature drops arising from radial ’ 
heat flows. Vapor pressure values associated with measured 
temperatures and with the temperature minimum arc , also dis- ‘ 
played in Figure 3.25. The temperature reaches a minimum: 
between the evaporator and condenser. / In addition a pressure 
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recovery of .5 has occurred which compares favorably with 
theoretical predictions. to be. discussed in a later section. 
Rankcn and Kemmc compared data, obtained by dohdansky and Schins 
[90] using a lead- tantalum heat pipe, with predicted values'-* 
for both laminar and turbulent flow. Figure 3. 26. Blasius' 
turbulent flow- equation gave considerably better results 
than the laminar flow equation, which was expected, because 
the axial flow Reynolds' numbers were well over 1000 for 
this experiment. McKinney [35, 91] conducted extensive 
tests on a series of water heat pipes in a moderate tempera- ‘ 
turc range (temperatures up to 400°F). Among other con- 
clusions, which have already been amply discussed, he found 
that the magnitude of the radial Reynolds number had little 
or no effect on heat pipe operation. 

Tien [92] measured the axial temperature distribution 
along the outside of a water-ethanol heat pipe. He deter- 
mined the pressure inside of the pipe and compared this with 
the pressure which had to prevail if pure ethanol occupied 
the condenser and purt water occupied the evaporator. He. 
concluded that separation into pure components in a heat 
pipe is extremely difficult, if not impossible, to achieve. 
Instead he found that if the initial composition was rich ' 
in ethanol, the data attested to the existence. of a water- 
ethanol mixture in the evaporator while nearly pure ethanol 
(i.e. , the azeotropic mixture) occupied the condenser sec- : ; -• .- 
tion. All of Tien's .data, however, were obtained with the 
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pipe operating vertically with the evaporator below the con- 
denser. If the pipe contains excess liquid such an arrange- 
ment is usually referred to as a reflux condenser because 
gravity forces, instead of capillary forces, can always 
return the condensate to the evaporator. In such cases 
results would have to be viewed with some caution because 
the wick structure may or may not significantly alter the • 
liquid vapor equilibrium conditions. 


Heat Pi] 


ications 


Theoretically the heat pipe may be applied to almost a 
limitless number of thermal transport problems, which in 
general, can be subdivided into four broad topical categories 
(8, 9, 10, 20] depending on the particular feature of a heat 
pipe which is to he exploited. These areas of possible 
application arc: 1) Temperature Flattening, 2) Source-Sink.. 

Separation, 3) Meat Flux Transformation, and 4) Constant Flux 
Production. ^ 

The temperature equalising feature of the heat pipe has 
prompted numerous suggestions .and actual uses for the main- 
tenance of a desired constant temperature environment. Much : ‘ 
of the emphasis has been, and is still, focused on the prob- 
lems of thermal control of spacecraft. It is well known that 
large temperature variations may occur on the- surface of a : 
spacecraft resulting from nonuniform heating of the craf.t . , 
These temperature variations can cause a host of problems 
including undesirable thermal stresses. Kattoff [53 together" 
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with several other investigators [7, 93, 94, 95 3 recommended 
the use of long heat pipes, wrapped around the circumference 
of a spacecraft, to accomplish the necessary equalization of 
the temperature distribution. Naturally the evaporator sec- 
tions e-f the circular pipes would have to face the sun while 
the radiation cooled condensers have to remain in the shadow 
of the craft. According to Katzoff such an arrangement of 
heat pipes would reduce the temperature variations around a 
spacecraft with a circumference of ten meters from 275°K to ' 
44®K. Anand et al. [96, 973 reported the successful employment 
of two Frecn-11 heat pipes which reduced the temperature differ- 
ences between transponders located in different parts of a 
Geos II spacecraft. Deyerall utilized the isothermal walls 
of a heat pipe to measure the total hemispherical emissivity 
of variously prepared surfaces [25, 983. Several heat pipe 
containers were plated or sprayed with different materials, and 
the resulting emissivities were .determined over a wide range of 
temperature with an estimated. accuracy of ±21. Schretzmann 
[993 employed an isothermal surface of a heat pipe as a metal . 
source in a study of the effect of electromagnetic fields on ■ 
the evaporation of metals. Bohdansky and Schins [100] used • . 
a heat pipe for the determination of pressure-temperature re- ' 
lations of metal vapors at high temperatures and pressures. 
Feldman and Whiting [9]have suggested the construction of an iso- 
thermal flat plate for the installation of electronic components. 
They conceived a sandwich type plate, the interior filled with 
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many interconnecting honeycombed heat pipes, which would 
rapidly distribute any localized heat flow and maintain the 
plate at a uniform temperature. 

Another use of the heat pipe would allow separation of 
the heat source from the heat sink. Again, possible space- 
craft applications- of hear pipes have received, considerable 
attention since the heat source, for instance electronic 
components, is often located in the interior of. the craft 
and the waste heat must often be transferred over some dis- 
tance for ultimate rejection to outer space. Moreover, any 
sizeable temperature drop between the source and the radi- 
ator may induce a significant weight penalty because a 
larger radiator area will be required at lower temperatures 
to dissipate the same amount of energy. Since the heat pipe 
is both light weight and nearly isothermal it appears to 
provide ideal solutions to many thermal dissipation problems 
encountered in spacecraft. Many investigators [46, 101, 13, 
102, 103, 104, 61) have considered its application to such 
energy dissipation systems and have cited numerous advantages 
including greater heat transfer per unit weight, and some 
degree of meteor protection when used in parallel arrangement 
Werner and Carlson f 1 0 4 J have reported that heat pipes can 
operate sixty times more effectively as radiators than solid 
rods based on heat transfer rates per unit weight. Dcverall 
and Kemnc[23] support this claim by reporting that- a silver 
heat pipe is 520 times more effective than an equivalent 




solid tantalum rod. A team of investigators of RCA [20, 51] 

f * * 

has designed, constructed and successfully tested a space ' 
radiator composed of one hundred individual heat pipes. The 
system which weighs less than twelve pounds is capable of 
rejecting 50,000 watts of thermal energy at a temperature 
of 1420°F. McKinney [35] recommended the installation of 
cryogenic heat pipes around cryogenic storage tanks. Ii 
the heat pipes contained a working fluid having a lower 
boiling point than the stored fluid it should be feasible to 
transfer the heat, leaked froir the immediate surroundings 
into the storage vessel, to a remote sink for dissipation. 
Researchers at Los Alamos Scientific Laboratory explored a 
unique scheme ir. which the construction of a heat pipe plasma 
oven was proposed. The disposal of waste heat given up by 
electronic components at remote locations has already been 
mentioned several times [11, 49j. . 

Another feature of the heat pipe, which has generated ; 
much enthusiasm, is its ability for thermal flux transforMa- 
tion. Heat addition to and heat removal from a heat pipe 
are feasible across heat flow areas of different size. This 
potential for heat flux transformation has stimulated thermi- 
onic specialists to consider the conversion of low heat, 
fluxes generated by radioactive isotopes, for instance, into- 
a sufficiently high heat flux which is required for the y 
operation of a thermionic converter. "Lecfcr . [64 ] achieved a 
flux transformation with a flux concentration ratio cf . 
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approximately ten to one equivalent to an output flux density 
of 2S0 watts/cm 2 which is more than adequate to meet the 
requirements of a thermionic converter. Other investigators 
also reported measured flux conversions of a ratio of ten to 
one [107] and even twelve to one [74, 683. A number c?f 
papers have appeared in the literature [23, 64, 65, 67, 71, 

108 to 117] extolling the virtues of heat pipes When used in 
conjunction with thermionic converters. Numerous heat pipe- 
thermionic converter assemblies have been built and tested. 

The results have generated much optimism and it is believed 
that the heac pipe applied to thermionic converters wiM 
reach technical maturity in the near future. 

Heat pipes may also be used to "flatten" flux variations 
supplied by an unsteady heat source. Researchers at RCA [107] 
have developed a "classified" radioisotope powered heat pipe 
which supposedly maintains a constant thermal output flux inde 
pendent of variations in thermal input flux tor at least 
three half lives; A qualified, thermodynamicist might want 
to study these more exotic schemes in view of their com- 
patibility with the requirements of the Second Law of Thermo- 
dynamics. 

Two particular potential heat pipe applications have been 
proposed by investigators at the Lawrence Radiation Laboratory 
Hampel and Koopman [32 ] suggested the. utilization of the 
heat pipe concept for the. control of small, fast-spcctrum, 
high temperature reactors. Their scheme is based on the 




capability of the heat pipe to respond to a sudden increase 
in heat flux with an increase in evaporation rate. If the 
evaporator ends of a sufficient number of heat pipes, filled 
with special. working fluids having a negative void coeffici- 
ent, were installed in the reactor core, an increase in core 
heat flux would result in a decrease of the total mass con- 
tained in the core (evaporators). Such a hypothetical mass ' 
transfer device may eventually be used to provide some sort 
of reactor safety control. Werner (118] proposed to employ 
lithium heat pipes as tritium producers in the blanket 
structure of a. reactor, the scheme calls for the transport 
of the tritium within the heat pipe to an accessible proces- 
sing point outside of the blanket where it could be removed 
bv diffusion or equivalent means and then be used to 
replenish the tritium consumed in the core. 

Heat Pi pe Control 

Inasmuch as the heat pipe transfers energy between two 
points utilizing a continuous mass circulation, it is appar- 
ent that some degree of heat pipe control may be exerted by 
controlling the mass flow. Katzoff (5) suggested several 
concepts to accomplish thermal control witn the heat pipe 
serving as a variable thermal conductor. One technique 
involves the intentional introduction of a noncondensable 
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gas into the vapor space. - As discussed earlier, the gas 
tends to collect in the condenser where it forms a relatively 
stagnant gas zone which effectively eliminates any working 


fluid condensation. The length of the gas. zone, of course, 
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depends on the working pressure in the pipe; t.he zone length 
decreasing with increasing pressures. The gas layer can be 
exploited in several ways. For example, suppose it is 
desired to always furnish solar energy to an instrument 
inside of a spacecraft regardless of the craft's orientation 
with respect to the sun. The task could be accomplished by 
mounting the instrument in the center section of a heat pipe „ 
containing some noncondensablc gas. If the amount of gas 
was such that it filled slightly less than half of the pipe, 
then the instrument would always receive energy from the 
sun. The end of the heat pipe on the shaded side of the 
spacecraft should remain inoperative due to the blocking 
effect of inert gas and thus minimizing heat losses. The 
opposite problem could be solved by using a partially dry 
heat pipe. Now the dry part of th? heat pipe, and hence the 
inoperative part, could always face the sun, and the instru- 
ment which might generate heat could reject this - heat by . 
radiation from the shaded side of the spacecraft. to outer 
space. Besides Katzoff, Wyatt [1 193 and An and c t a 1 . 1 7 9 J have 
suggested the use of thermostatically, controlled valves and/*'"-' 
or bellows tc supply or withdraw the noncondensablc gas. 

Such an arrangement would allow the. effective condenser area ; 
to be varied independently of the operating pressure prevail-.’-, 
ing in the pipe. Katzofl* has aljo recommended a control 
technique which involves the interruption of the liquid flow 
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Figure, 3. 27 Schematic Cross Section of Variable Conductance 
Space Suit Shell, Shlosinger [42, 43, 44] 


in the wick. This control technique would be practical only 
in wicking structures in which the hulk of the flow takes 
place in an artery. A thermostatically controlled valve 
could he used to impede, or stop the flow in the artery, in 
response to the control requirements. Anand et al. [79] and 
Shlosinger et al. [44] have also considered a similar 
arrangement for the control of the vapor flow. Shlosinger 
[42,43] used this concept into the design of a variable 
conductance space suit. His design is illustrated in 
Figure 3.27. The design concept provides for both insula- 
tion when the valve is closed and for neat rejection when 
the valve is open during periods of high metabolic heat 
generation of the space suit wearer. Experimental results . 





were interpreted to consider such a technique as. feasible, 
but calling for ouch more research. Anand et al . [79 j" recom- 
mended a control technique involving the use of. two fluids 
whose pressure-temperature curve intersect at. the desired 
operating temperature. In anticipation of the difficulty . 
of finding such fluids, this method of control has received 
only minor attention. ; Workers at Honeywell [27} have studied 
the concept of bellow controlled feeder wicks which are ‘ 
either in contact with the heated surface thus providing a 
path for the condensate - return , of the contact is inter- 
rupted, thus preventing fluid return and hereby shutting 
off the heat pipe action. Heat pipes acting as "thermal 
switches" were built primarily in an effort to develop vari- 
able conductance walls. Conductivity ratios of 150 have 
been reported depending on the position of the feeder wick. 

Obviously a great many possibilities for heat pipe con- 
trol exist; however, most of these concepts remain still in 
the dreaming stage, and few controllable heat pipes have 
actually been built to date. Much more research will have 
to be devoted to this particular area of heat pipe applica- 
tion so that eventually some thermal problems can be solved 
with controllable heat pipes. . - ' 


Heat Pipe Theory - . . 

The heat pipe theory developed up to date is due, in 
large part, to a stimulus provided in a ..theoretical study : V- 
performed by Cotter HJ. He formulated the governing > 







104 


equations describing the processes taking place in the heat 
pipe and also developed a model to predict the capillary 
limit of heat pipe operation. His results have been used by 
other investigators who have often simplified or modified 
his equations to suit their particular assumptions and/or 
geometries. Because this has been done so extensively, 
Cotter's analysis will be presented in necessary detail. 

The description, of the analysis will be followed by a number 
of modified theories and by several other analytical- 
approaches, all of which are concerned with the capillary ‘ 
limit. Derivations for the predictions of the other limits 
of heat pipe operation are given subsequently and finally, 
several analyses which deal with specific problems conclude . 
the analysis section. . ' 

Cotter studied a heat pipe as shown in Figure 3.28. 

The capillary structure is assumed to have a pore radius, 
r c , and to be completely saturated with a working fluid. 

The radius of curvature of the meniscus surface is thought 
to be dependent on the distance, z, and the pressure differ- 
ence across a surface as a function of position is given by 


P y (c) - P.(z) 2o/r(z) = 2ocos0/r c 


(3-10) 


where r(z) is the local meniscus radius of curvature. In. 
equation (3-10), of course, it is assumed that the meniscus 
is represented by one radius of curvature only. If, in 
fact, a complex meniscus shape were formed, the term 
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Figure 3.28 Cylindrical Heat Pipe Structure, Cotter [ 4 ] 


( 2/r(z )} in equation (3-10) should be replaced by ♦ -j—j 
where r, and r 2 arc the two radii of curvature necessary to 
describe a three dimensional surface. If a heat pipe is to 
be operated in a gravity field, then the maximum length, 
r max , of the pipe is restricted by the lifting ability of 
the combined wick and liquid system. The maximum length is 
given by the well known relation 


’max 


cosB 


2ocos0 


p l8 r c 


(3 - 1 1 ; 


Cotter next considered the steady state operation of a 
heat pipe and applied the conservation of mass principle to 
arrive at the relation ^ 

n\ r (z) ♦ m^Cz) - 0 (3- 12 

where the two mass flow rates are both positive in the plus 


z-direction . 





The pressure gradient in the liquid was determined from 
the-Navicr Stokes equation for steady, incompressible, con- 
stant viscosity flow. By neglecting the inertial tern and 
by modifying the viscous terra with a dimensionless constant . 
to account for different wicks, the pressure gradient in the 
z direction was found to "be . 


dP. hu, m £ (z) 

~ji “ °ts cosB ■ 

B 


T i — 7S 5 

jV r vJ°* c r c 


(3-13) 


where the dimensionless constant, b, is defined as er^/K 
and has a value of approximately 8 for non-connected parallel 
cylindrical pores, «nd 10-20 for more realistic capillary 
structures with tortuous and interconnected pores. • The 
radial pressure gradient has been assumed negligible as will 
be the case for long thin pipes; that is, for those pipes 
for which r 1 » r 7 . 

V W . : 

The pressure gradient in the vapor was found by employ- 
ing the results of Yuan and Finkclstcin [120] and Knight and 
Mclntccr (121). These authors assumed incompressible laminar 
flow and uniform injection or suction at the vapor space 
boundary. The applicability of these results is based on 
the value of the radial Reynolds number, Rc* r , which is 
defined by * 


Rc 


°v r v v r 1 d * v 

a 2ir i r; ~Zz 


(3-14) 


It is positive for evaporation and negative for/*'ondcnsation. 
I Rc r [ << 1 , the vapor flow is dominated by viscous forces 


For 


) - \ 


--3 
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and the velocity profile approximates the usual parabolic 
shape for Foiseuillc flow. For this ease the vapor pressure" 


gradient is given by 




8 u y m 


rf (» * f R 'r - jf K * — ) ' »•«>' 


For | Re r | >> 1, however, the flow is qualitatively differ- 
ent in evaporator and condenser. For high evaporation rates, 
the velocity profile is not parabolic but is proportional to 
cos j * . which was verified experimentally by h’ageman 
and Cuevara [121J, and the pressure decreases in the direc- 
tion of flow. For high condensation rates, on the other 
hand,' the velocity profile is nearly constant across the 
vapor space with the transition to zero velocity occurring 
in a thin layer near the wall, and the pressure increases 
in the direction of motion due to a partial dynamic recovery 
in the decelerating flow. For { Ro ^ j «®, the pressure gradi- 
ent is found to be 




iv. 


1F V -S m v dm y 

■3? s ^TV 


(3-10) 


where S «• 1 for evaporation and S .= 4/n i for .condensation. 
Furthermore, Cotter hypothesized that for some situations - 
where the Rc r * 0 and the average vapor velocity is high,, 
as might exist in heat pipes with a long adiabatic section, 
fully developed turbulent flow may occur. - For such a case. 
Cotter recommended using ' 


■••Ik* 


■; :065S w v * Re 

p f J 
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in lieu of equation (3-15). , 

The relation between vapor and liquid pressures and the ■ 
mass flow rates is next given by a formula supplied by the 


kinetic theory. 


dm^ dm v 

“37 r ' “3z 
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(3-18) 


where a is a numerical factor which includes the probability 
of condensation of a vapor molecule and the surface rough- 
ness of the miniscus. The value of o is very nearly unity. 

Cotter next considered the conservation principle of 
energy and by neglecting radiative and conductive contribu- 
tions he arrived at 


Q(z) * h fg n v (z) 


(3- 19) 


Finally, Cotter coupled the heat pipe to the surround- 
ing environment by expressions accounting for heat fluxes • 
or imposed temperature conditions respectively 


JQO0 . 


2ir Vr^’ r n^ “ H(z. T n .Q) 1 (3-20) ' 


T(z,r D ) = T(z,r y ) ♦ ^ (j^.in JL * ^ fen *(3-21) 

where equation (3-21) relates the temperature on- the outer 
surface of the heat pipe to the temperature at the liquid 


vapor meniscus. 


Equations (3*10) through (3-21) are generalized heat 
pipe relations as presented by Cotter. 1 These equations can 
he solved if the vapor has a nearly uniform temperature, T„ 
throughout the vapor space. If equation (3-21) is solved 
for, Tp[ T ( 2 » r pj J » then 11 in equation (3*20) may he expressed 
as a function of length, z, temperature distribution in the 
vapor T y jV(z,r v )J , and the heat flow, Q. Also, if the heat 
flow through the pipe ends is denoted by F e (T y ) and F, (T y J,, 
then the average vapor temperature, T # , and 'the correspond- 
ing axial heat flux distribution , Q 0 (z) are determined from 
the relat ions C — 

' ■ ' dQ p ' : ' ... ' -V- ’ " ' ‘ 

-Xi “ h ( 2 . t # , q 0 ) . 

• (3-22; 

<M°) - M T o> ; V o ( 1 ) - F,(T 0 ) • 

With this approximation for the heat flux, the mass flow 

... - ** t 

rates may be determined from equations (3-12) and (3-19). 
The vapor mass flow rate can then be used in the appropri- 
ate equation, (3-1S) , (3-16), or (3-17) to determine the., 
pressure distribution in the vapor. . Also equation (3-13) 
may be integrated to yield the axial pressure distribution 
in the liquid. Finally , equation (3-18) can be solved for 
the vapor pressure of the liquid to within an additive 
constant. Cotter recommends using P(T ? ) for this constant. 
Since the vapor pressure is -a know; function of liquid ; 
surface temperature, T^f:) may be determined. . The ^ 




consistency of this procedure can be verified if the axial 


variation of is small compared to T 0 . 

Cotter used this procedure for the special case of 
constant heat addition along the evaporator and constant 
heat removal along the condenser. Thus he obtained 


Q«(*) " h fg™ v (z) 
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(3-23) 


where Q e ii> the total heat input to the evaporator. By 
employing equations (3-23) in (S-1S) and (3-16) and neglect- 
ing the term Refill the former, and assuming the vapor 
density, p y , to be constant in both, a straightforward inte- 


gration yields 
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(3-24) 


Likewise, integration of equation (3-13) yields 

. " P t Blc05S ♦ ^ CS-25) 

2 "[V r vjV :r c h f* 

The temperature difference in the vapor space can be 
obtained from equation (3-18) leading to 


lQ e /RT*/2xM 
'el 1 ' 1 .) 


(3-26), 
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Anil finally by considering the vapor as an ideal gas, appli 
cation of the Clapeyron - Clausius equation and neglecting 
the volume of the liquid phase results in 


RT fl 1 AP. 


AT » 
v 


The preceding equations were applied to a horizontal sodium 
heat pipe and the results of these calculations arc listed 
below. The left hand column indicates the pertinent values 
used and the right hand column shows the calculated param- 
ctcrs. ’ ‘ ....... 
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An inspection of t lie computed values explains why the heat 
pipe has generated so much enthusiasm among thermal engi- 
neers. A relatively large energy transport is accomplished 
with an almost negligible temperature drop.. 

Cotter next considered the maximum heat transport which 
is possible. in a heat pipe which is -limited .only by the 
pumping ability of a wick in association with a given work- 
ing fluid. The maximum pressure difference which can be. 


supported between the liquid and vapor is achieved when the 
meniscus radius in the evaporator achieves a mii:.imum value, 
i.c., equal to r». Therefore equation (3-10) becomes 

P y (t) - P t (2) 1 -- r ' 9 - 8 (3-28) 

which must be satisfied for all positions, z. Tn particular 

i 

for z ■ 0, equation (3-28) nay be rewritten with the aid of 
equations (3-24) and (3-25) to give the capillary limiting 
condition for the total axial heat flux 
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(3-29) 

where - 

P v (0) - P £ (0) P v (l) - AP v ^AP £ - P t (i)‘- AP £ - AP v (3-30) 

has been used to simplify the left hand side of (5-28) prior 
to the substitution of equations (3-24) and (3-25). 

Cotter next determined the optimum - -•.pillary pore 
radius, r c , in terms of Q c by considering equation (3-29) 
where the equality sign is used and the first term is simply 
-AP y . Now since AP y is directly proportional to ‘Q (or Q c *) 
the heat transported is a maximum if AP y is a maximum. Thus 
equation (5-29) becomes \ " 
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Kith equation (3-31), the pressure drop, LV yt becomes a 
maximum when 


(3-31 
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The maximum heat transport is determined by converting cqua 
tion (3-28) to the form - • f 


*p t - sr v 
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(3-33 


Now for Rc f << 1, and with no hydrostatic contribution to 
AP^, equation (3-33) is used in conjunction with equations 
(3-24), (3-25) , and (5-32) to solve for Q e , leading to 

- o?cos 2 8 r " (r ? -r ? ) 
n 2 s v l w v) 
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( 5 - 34 ) 


Obviously Q becomes a rcoximum when the term r ,, fr 2 -r*l - 

v ( w *vj 

assumes a maximum. Simple calculus shows that this situa- 
tion is achieved when ■ 
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Substituting equation (5-55) into (3-34) finally yields 

£ h fg qc°sO pP v P t e]V 2 
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Inserting equations (5-36) and (5-35) into (3-32) yields 
for .the optimum pore size for Re '«.l \ >V •’ ’. v • 


r • r 
c w 


*bu £ P 1 i/j ' • 


(3-37) ; 


The analogous expressions to equations (3-56) and (3-37) for 


Re >> 1 
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Cotter's, model for Re r >> 1 predicts that the' maximum 
extent of pressure recovery is fixed at 4/n 2 corresponding • 
to 40.5 percent of thc.drop that occurs in the evaporator - 
regardless of the amount of heat transferred. Therefore, 
the profile illustrated in Figure 3.29 represents the. vapor 
pressure distribution for any case which is dominated by 
inertial forces. It should be recalled that the profile in'. 
Figure 3.29 is dependent on all the assumptions made by..- . 

Cotter in his analysis and in particular, on the assumption 
that the vapor is incompressible. Parker and Hanson [24 3* have . 
written; a computer code in which they consider tho vapor, to •' 

be compressible, and consequently . treats the' vapor, density - 
as a variable. A comparison of ' their ' . results with - ' : 

Cotter's predictions demonstrates excellent agreement, for 
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the pressure profile in the evaporator section. However, 
a significant departure from Cotter's prediction was found 
to exist in the condenser section. Here, the extent of 
pressure recovery was dependent on the heat transported by 
the pipe. Figure '5.30 illustrates the predicted pressure ' 
recovery as a function of the heat transported for a partic- 
ular sodium he.=t pipe. The extent of pressure recovery is 
seen to be considerably larger than that predicted by ’ W 
Cotter for largo heat transport rates and hence the netiend 
to end pressure drop in the heat pipe will be smaller! When 
compressibility effects arc important. Cotter's analysis 
underpredicts the maximum heat transfer capability of the 
heat pipe, and it may be used to provide a conservative 
estimate of the capability of the pipe. Of course, an even 
more conservative estimate could be obtained by neglecting 
pressure recovery altogether. 

/ 

Frnst [54] has disputed the validity of the meniscus 
boundary conditions employed by Cotter for the special case 
of large radial Reynolds numbers, zero gravity, and incom- 
pressible vapor flow. ; It should be recalled that Cotter 
assumed a meniscus profile of the type shown in Figure 3.31. 
Using this profile, he went oh to optimize the capillary- 
radius and heat transfer rate ar.d arrived at equations 
(3-29), (3-38), and (3-39). Krnst has calculated the pres- ^ 
sure profiles based on these equations for the -special case 
of a sodium heat pipe operated at 700°C with r • l'cm. 
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1 • ■ 50 cm. and 1 ■ 50 cm.. His profiles arc illustrated ' 

c v a . w . 

in Figure 3.32 by the lower set of curves./ Figure 3.32 
shows that at the condition of maximum heat transfer pre-r 
dieted by Cotter, the pressure gradient in the liquid is 
smaller than that in the vapor. Consequently, the pressure 
of the liquid in the condenser is greater than the pressure 
in the vapor. Ernst has concluded that for such a situa- 
tion, a meniscus profile of the type illustrated in Figure, 
3.33 must exist. However, he regarded this meniscus profile 
as unrealistic and instead proposed the profile illustrated 

" ~ ; . ■ " V ■ * . 

in Figure 3.34 for the non-optimized capillary case. ‘ In 

this revised profile, pressure equality between the liquid, 
and the vapor is assumed to occur at the interface between 
the condenser and the evaporator instead of at the con- 
denser end as assumed by Cotter. For the special case of a 
.wire screen capillary structure, Ernst has derived the 
following equations for the pressure balance in the evapora- 
tor and condenser, the optimum mesh opening half-width, and 
the maximum heat transfer capability , respectively. * ' 


Evaporator: 


‘ b |i # Ql 


2ocos8 


Condenser: 
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(3-40) 


(3:41) 




Optimum mesh opening half-width: 



Maximum hc3t transfer: 



where r c * d(l+K) and K » d,/d has been used in deriving 
equations (3-40) through (3-43) and r* * (2/3)r*. has been 
applied to the derivation of equations (3-42) and (3-43).* 
Ernst also applied equations (3-40) through (3-43) to the 
same 700°C sodium heat pipe discussed earlier. His pres- 
sure profiles are illustrated by the upper set of curves 
in Figure 3.32. The pressure equality between the liquid 
and vapor in the condenser implies a meniscus profile of 
the type shown in Figure 3.35. A comparison of the revised 
maximum heat transport equation (3-43), with the expression 
derived by Cotter, equation (3- 38) , lead Ernst to believe 
that Cotter's expression tends to ovcr-prcdict the maximum 
heat transfer capability of the heat pipe. In particular, 
for the 700°C sodium heat pipe for which'.' the pressure pro-.;, 
files in Figure 3.32 were calculated; this over-prediction 
is approximately 401. The extent of over-prediction, of 
course, will vary for different heat pipe geometries. V 
Ernst also claimed that the maximum heat transfer capability 
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of a heat pipe can be increased if the'* wick in the evapora- . 
tor is. different fro.n the wick in the condenser. If the 
optimum mesh half-width in the evaporator and condenser are 
respectively . ' \ 
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(3-45) 


then the expression for the maximum heat transfer, given 
previously by equation (3-43)., becomes 
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(3-46) 


A comparison of equations (3-43) and (3.-46) for a typical 
case where c « .8 and 1 » l e accompanied by some calcula- 

tions shows that the use of different wicks in the evapora- 
tor and condenser may increase the maximum heat transport 
capability of heat pipes by over 155. J'he percentage ' 
increase varies from pipe to pipe depending on the evapora- 
tor and condenser lengths, and on the porosity of the 
capillary structure. 

Many investigators have suggested various alternative 
express ions for the pressure drop in the gas and liquid for 
different heat pipe configurations which may be used in 
equation (3-29) instead of the pressure drops predicted by 
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Cotter as given by equations (3*24) and (3-2S). Bohdansky et al 
(46), for instance, recommended the following expressions 
for a heat pipe which employs axial channels for the liquid 
return and for which zero g, laminar, incompressible flow 
is assumed in both the liquid and the vapor: 
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where K is a channel shape correct ion . term which is approxi- 
mately 1.3 for channels of rectangular cross section and 
depth equal to twice the width, liquations (3-47) and (5-48) 
Mere used to derive the following relations for the maximum 


heat transfer rate: 


u r v* h fg ( K*A ] »/j 
b a 1- [2v v > vj (-3-49) 


where ’a' is a channel shape factor which has a value cf 
about 2 for the channel discussed above, and 'A* is another 
dimensionless number defined by v ! i.'.v-: 


r ch ; 


(3-50) 


Notice that A is dependent on the number and spacing of the 
axial channels and has a maximum value of rc . Equation 
(3-49) is of limited importance to the design engineer, ? . 
because of the large number of dimensionless constants ,*' v . 
which must be employed, and the lack of a method for the 



determination of these constants. Busse ct al. [47j modified 
equation (3-49) so that all the dimensionless parameters are 
incorporated into one constant. His equation for the maximum 
heat flux is - . 


'max 


13 & r v ohf cosfl 
r * — - 


eff 


v '/> 


(3-SI) 


where 5 is this dimensionless constant. If the channel 
depth is 1,, the channel width is 1 2 , and the dimension 
between channels is 1 j , then 6 is approximately .9 for 
lj *.21j = 21, and is approximately .5 tor 1, = 1 1,. 


The term l c fj- is defined bv Busse as 
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and is the so-called effective length of the, heat pipe 
where q(z) is the average axial heat flux at any position, 
z, averaged over the pine cross section. The inclusion of 
equation (3-52) into equation (3-51) would, however, cause 
the term of interest, i.e., q m . to cancel. A more 
straightforward approach is to simply replace l eff in (3-51) 
by 1, the total heat pipe length.- . v, . 

Bondansky and Schins (90) recommended the following 
expressions for the vapor and liquid pressure drop for 
turbulent flow: ■ z -v.. ;."V -rfL- '■ 
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liquation C 3 - 5 3 > is obtained by a straightforward integra- | 


i - 

tion of equation (5-17). The expression for is obtained ", ; 

in a similar manner but has been modified to account for j. 



channel shape by the dimensionless K tern. As suggested j 

above, 1 ^ should perhaps be replaced by the total heat j 

- „ . . . ■ ■ . ■ ‘ ■* . - ( 

pipe length.' ‘ Equation (3-S4) is very likely of limited : 

value because the liquid flow will generally be laminar even j 

for very high heat, transfer rates. 

! 

Bussc 1 1 2 3 ] considered laminar vapor flow in a cylin- 
drical heat pipe which had an adiabatic section separating 
the evaporator from the condenser. For the case of constant 
heat addition and removal the vapor flow was described by 
the Navier-Stokcs equation which was solved by approximating 
the axial velocity profile by a fourth power polynomial of 



the radius. The analysis furnished a velocity profile which 
was relatively constant along the evaporator; that is, it 
approached the Poisciulle profile in the adiabatic section, 
and deviated cons idcrably. : from the Poise iul lc prof ile in the 
condenser. The pressure drop.in the evaporator is given by 


V 2 > - f> v 
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(3-55) 



where V is t the axial, velocity averaged over the cross sec- 
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Equation (3-55) may be approximated with an error of less 
than 11 by . ' . .'• -.*■ , . 


P v ( 2 > * V 0 > 


-4u V 1 T .61ReJ , 2 •-, - 

“ *— [ l+ * 6 IRc r + T:^F;J T^J (3 ~ 57) 


The pressure drop in the adiabatic section is Riven by 


~8p V z r Re Pq *>- ’] 

P v (2)> P/0) .-— U?- 1 * ; -J- (A-a) (A^aMj (3-58) 


where Re is the axial Reynolds number defined by 

Re = 2r V P../M -and 'a* is a correction to .the Poiseiulle 
v m v v . «. 

* ^ -.V 

velocity profile defined by 


a cxp(-22a/75) = A exp (- 22A/75)cxp(- 1 44c/Sr v Re) (3-59) 


Typical values of 'a* arc illustrated in Figure 3.36 for 

several radial Reynolds numbers. Equation (3-58) can be 

* » 

approximated with an error of about It by 


P v (z) - P v (0) = 


-3p v V m s .l06Rc r r-cxp(-302/r v Re)T - 

— TyT- 1 * :/r v kc i : 1 


j (5-60) 


The pressure distribution in the condenser section is given 

by , ; • •' ■/. ■ V . •= ; '/'.TV; . 

> v (o • |> V C=-1 C ) 

■ -v v-': ; (3-6i). 





Figure 3.36 Velocity Profile Correction vs. Position 
in. Heat Pipe, Busse C 1 23 3 

where now -'a' is. found by the solution of 
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(3-63 


Here a o is the velocity prpfile correction at the beginning 
of the condenser section. From Figure 3.36 it is' seen that 
0 - a o - * 655 * Busse plotted the dimensionless pressure 
distribution in the condenser section as a function of 
dimensionless length as illustrated in Figure 3.37. . Figure 
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Figure 3.37 Pressure Profiles in Condenser, Busse £123} 

3.37 represents the special case where the adiabatic section 
is so long that the flow at the beginning of the condenser 
has assumed a Poiseuille velocity profile, i.e., a Q = 0. 

It also shows that up to a certain radial Reynolds number a 
pressure minimum occurs at the end of the condenser and' that 
increasing condensation rates tend to move this minimum 
toward the beginning of the condenser; The dislocation of 
this minimum indicates a partial pressure recovery. The’ 
pressure drop of the vapor which should be employed in an 
overall pressure balance equation; such as equation (3- 29) , 
depends on the net heat transfer of the heat pipe. Busse 
has mentioned that for strong heating and cooling rates, a.', ., 
pressure equality exists in the vapor and the liquid phase' 
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in the proximity of the beginning of . the condenser:' ^The 

« - ‘ . ,* ■ _ » ■ •* n 

observation was substantiated by Bohdansky et al. [763; ..recall 
that Ernst also had noted this equality in pressures. For 
such a situation, only the vapor pressure drop in the 
evaporator and’ adiabatic section has to be taken into ' 
account for the calculation of the maximum heat flow. Thus 
a combination of equations (3-55) and (3-58) yields the 
total pressure drop in the vapor expressed by Bus sc as 


AP 




. ; (3-64) 


where 
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F * 


8a 


- ’ % 
tip - ' 
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Tor 


(3-65) 


The factor F can be approximated with an error of less than 

li by .. ‘‘ ‘ . ‘ ■’ • 


1 . 7 Re 


F = i * cx n 


-7.5 1. 


Rc r ^ , 


(3-66) 


and has values ranging from .61 to .81.' On the other hand,’ 
for low heating and cooling rates, i.e., small Re r , the ' 
pressure equality between the vapor and liquid will generally 
exist at the far end of the condenser. For this case knowl- 
edge of the pressure drop in the entire vapor. space is ’ 

required and Bussc combined equations (3-55), (3-58) , and ' 
(3-61) to arrive at the simple result- 'y:]: ’ 
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The relation was previously obtained by Cotter,' equation' 
(3-24), for a heat pipe with no adiabatic section, i.e., 

''here the. bracketed lengths in equation (3-67) were 
replaced by the. .total "heat pipe length. It is interesting 
to note that .-the length-ydf -the adiabatic section has more 
influence on the^p.^ossurc drop than the contributions caused 
by the long ths;V6f -the evaporator or condenser. 

Ilaskin (2SJ has dealt with a heat pipe containing an 
adiabatic section. He transformed equation (3-23) into 




2 Q e /i e 


V z > “ h f>v ( = ) > H 


; 0 < z < i e 

: l c ~ z - l c (5 * 68) 


a-z)Q e /i c ; 1 C < Z < i 


The pressure drop in the evaporator and condenser was 
obtained by substituting equation (3-68) into equation 
(3-15) and its subsequent integration. The pressure drop 
in the adiabatic section was found by assuming the existence 
of Poiscui 1 le flow. The net end to end pressure drop for 
RCj. <<c 1 is given by the sum of the individual contributions 


’°v y h f S 


(3-69) 


This expression is identical to equation (3-67) which was 
derived by Sitissc. For the case of Re f >> 1, Ilaskin used 
equations (.3-68) and (3-16) for the condenser-’and evaporator 
and applied equation (3-17) to the adiabatic section and 
obtained . . • -V r: •• 




Notice that both equations (3-69) and (5-70) reduce to 
Cotter's results, equation (3-24) , for the special case of 


0. The equation for the pressure drop in the liquid 


was found to be identical to that obtained by Cotter, equa- 
tion <3-25), with the exception that the length, 1, in the 
viscous term is replaced by ^1 * l^j. 


Schwartz fSR] has used Cotter's analysis and. its exten- 
sion by Uaski n to display the liquid and vapor velocity and 
the vapor mass flow rate as a function of the heat load for 
a water-stainless steel heat pipe. The wick consisted of 
two layers of 100 mesh screen and the pertinent dimensions 
5f the pipe were; = 3 in., l c = 3 in., l" = 3.5 in., and 
r v = .189 in. The results are presented in Figure 3.38, and 
arc typical for a low temperature heat pipe. The vapor 
velocity decreases with increasing heat load due to the 
large vapor density variation in the temperature range of 
interest . - ' 


Horner [118] has employed Cotter’s pressure balance, 
equation (3-29), for >> 1, and obtained in a computer 
calculation the axial heat . flux and Mach number as a 
function of length for a 1500°K lithium heat pipe with 
r D = .5 cm and r y =..4 cm. The results of these computa- 
tions are shown in Figure 3.59 illustrating- the high v T‘ v 











velocities aiul fluxes which may be obtained in a high tem- 
perature heat pipe. Figure 3.39, of course, is based on 
the equations for a capillary limited heat pipe and thus 
the flux rates may well exceed other heat pipe limits such 
as the entrainment or boiling limit. Werner and Carlson 1 1 0 S 3 
mollified equation (3-29) for the special case of a heat pipe 
with a capillary structure consisting of axial grooves 


covered by a single Ijyer of screen material. If the 
intervening wall between grooves is negligibly thin at the 
inner radius, then the pressure drop in the liquid may be 


ivrittcn as 


\ = PjjglcosB ♦ 


3u» Q l- 


I r r i wh r p„ r r ' 
l Vv v j fg t v c. 


(3-71) 


where fr -r] is the groove depth and r, is the capillary 
V w * t c 1 J 

radius (groove half width). Following a procedure identical 
to Cotter’s and assuming no gravity effect, Werner inserted 
equation (3-71) into equation (3-29) and arrived at the 
following expressions for the optimum capillary radius and 
the maximum heat transfer rate: * 


c,opt AC 


(3-72) 


Q max = A 17 


C,Optf 


(3-73) 


(1 -4/tr 2 ' 

JCT'V «. Ti ’ 


*£ 1 


— , C = 2ocos0 . 





In addition, he found that the maximum axial heat flux is 


obtained when r y /r w * 5/6 a value approximately 2.51 higher 
than the one found by Cotter (/2/3) . Notice that Werner 


assumed the same site for the capillary radius and the 


groove half-width. Consequently the screen mesh size and 


the groove width can not be selected independently Ham- 
pel and Xoopman [32] rectified this limitation by decoupl- 
ing the channel dimension from the pressure supporting mesh 
dimension . The liquid flow in the channels was treated 
simply as flow through an annulus and corrected adequately 


to account for the actual area available to the flow. In 


addition, a scaling factor ’a’ was introduced to account 


for the additional pressure drop induced by the channel 


configuration. For the case of Rc r • >> 1 and zero gravity 


effect, Hampel derived an equation similar to equation 


(3-29) : 


(l-4/n i )Q 2 •; . 2a Q 1 

JpT V h. ! ; 7 ( — r ~ • . r a ' 

v v ■ ' , * cp 0 »t h r r fr -r ] r 2 ♦ v z - — <£ 

fg V [ w vj ^ w v |ln|r w /r y j J 


2ocosC 


r w'V 


;; (3-74) 


where * d * is the screen mesh opening half-width and 'e* is 
the ratio of the active channel area to the total circum- 
ferential area of r y . By differentiating equation (3-74) 

with respect to r , the resul ting -ratio , r/r , was : found / 

v .. ... • • • y vr • - 

to no longer be a constant, i.e., V27J or 5/6 , but to depend 
on the heat flux, the operating temperature, and the ■ • 
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Figure 3.40 Optimized Relationship Between Vapor Radius and 



Wall Radius for Lithium Heat Pipe, Mcmpel and Koopraan [32] 


thermodynamic properties of the working fluid. Typical 

values for this ratio, r.,/r , are illustrated in Figure 3.40 

* . w 

for a high temperature lithium heat pipe. 

Anand et al. [14] and Anand [96] employed advantageously, 
equation (3-56). The total length, 1, was replaced by 
(l e + l c j and the condenser length, l c , was replaced by Q e /C 
where C is a condenser parameter defined by ■ 


C • 2h 


k w K; t w) 
( r P /r ») 


(3-75) 


This procedure led to a quadratic term in Q with the solution 

r, f, p„p, '|i/a Cl 

0._^ - U C» r.V hri ocose It c-l • — r-\ - 
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Equation (3-76) illustrates how a restraint on the axial heat 
transport capability is imposed by the radial heat flow in 
the condenser section. The condenser parameter, C, may be 
varied by wick flooding, introduction of noncondcnsablc 
gases, or by manually changing the surface area. Equation 
(3-76) reveals the possible applicability of these tech- \ ' . 
niques for heat pipe control. 

Broscns (110) started with equation (3-19) and con- 
sidered both the liquid and the vapor flows as laminar, 
steady and incompressible, i.c., 1’oiscuiiTe flow. The wick 
structure was thought to consist of *n’ cylindrical capil- 
lary tubes and pressure equality was assumed at the fur end 
of the condenser, for the case of zero gravity it ml perfect 
wetting of the wick, the expression for the maximum heat 
transfer was found to he 


[■■’ « i> f£ 

4-r t n 


0-77) 


The same result was independently arrived- at by Schlitnior and 
Wassnor 1124 1. Ilrosens optimized ctpiul ion (3-771 with respect 


to the capillary radius and . obl.-iimn 

3 v £ i/*. 

r c,opt r v n~v~ 


(3,78). 


from which resulted 


' 0 r v* h fg 


‘max /opt io 1 v 


v c f opt 


(5-79) 


The values predicted by equations (3-77)- and (3-79) tend to 
he much larger than those measured experimentally. The 
deviations arc caused by the many simplifying assumptions 
made during the derivation of these equations. 

Frank ct al . [30.) used Poiscuille flow for the l iquid and 
a modi r icd Pojseuillo flow for the vapor. He gave the 
vapo. pressure drop as • . 

8 “v ‘r V 


AP. 
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v it r. 


(3-80 


where If- is the effective flow length and extends from the 
mid-point of the evaporator to the mid-point of the con- 
denser, and the function, is giv?n by 




* = 
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00494 Re 




Re £ 2200 
Re > 2200 


(3-81 


Notice that is a discontinuous function at Re = 2200 (bv 
approximately S0$), so some reservations as to its applica- 
bility in this range are indeed justified. Prank made use 
of equation (3-80) to derive an expression for the maximum 
heat transfer capability of a- horizontal heat pipe: 





d'j' in the preceding expression is the hydraulic diameter of ’>•; 
the capillary "pores. Frank [26] next discussed the, applica- 
tion of equation (3-82) for a grooved heat pipc : and by 
assuming perfect wetting, i.b., a = 1, and by replacing 
by A, /A obtained for maximum heat flow •• • ' 
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. (3-85) 


Referring. to Figure 3.41, a number of auxiliary variables 
are defined. The mean radius of the grooves is;** . 


r m » r 
m 


K .2 v 2 


(3-86) •;;; 

while the dimensionless pitch of the grooves is given by • 


8 « 


w ’ 


(3-87) 


and the number of grooves is 


2 7T r 


m 


8 w 


(.3-88) 


The aspect ratio of the grooves is defined by 


6 

a “ w 


(3-89) 


and finally a depth ratio is given as 

6/2 


4- - - 


‘m 


■ - v ^ ' (3-90) 


ry inserting equations (3-86) through (5-90) into ' (3- 8S) the 
following expression is obtained: . . -■ • 
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Figure 3.41 Sketch Illustrating Design Variables 
in Grooved Heat Pipe, Frank et al. [50] 


= (I+V)“ 


(3-91) 


where 


S - — - 1 

f l 


; C = *N/B1 f ' (3-92) 

and 

s - 1 (3-93) 

f l . ; . , •• • 

hquation (3-91) was next extremized to yield the optimum • 

value of Q /r ,* and if a bar is used to denote the opti- 
max w ^ »* ; 

mum value, then the result is 


value of Q /r , J and if a bar is used to denote the opti- 
max w »•> : 


^max = c »* . ~ 

(]^) 3 1 '♦ 7 *. 


where is given implicitly by: 


fl-») (2-») (l*7i{t) ; 16 _v ^ 

r II ■ —2^2 6 v „ V 

<J> (lOip -•!-<>) * • 


(3-94) 


(3-95) 



















and o is given by ' 


- 2 ( 10 * - 1 . - * ) 


(3*96) 


In addition, since the optimum heat flow generally occurs 
when the vapor flow is turbulent, may be written as 


♦ * .0O494Re 


r fo u n 3 /«* 

V - .00494 111 (3-97) 
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The solutions of equations (5-94), (5-95), (3-96), and 
(3-97) were obtained by Frank and arc reproduced in Figures 

•• ' • j . .v " v . - . . 

3.42; 3.43, 3.44 and 5.45, respectively. The- procedure for 
optimizing the ratio Q m „ v ./r' J is an iterative one. Firs*, 
the effective flow length,. 1^, the operating temperature, 
the working fluid, and the pitch, 8, are selected. The 
pitch has a minimum value of unity and it is advisable to 
make it as small as possible. In general, the minimum ; 
value of the pitch will be imposed by machining and strength 
requirements. Once these values have been selected, an 
initial value of <2> v - 1 is assumed and may be found from 
.Figure 3.44. Using this value of tj;, the optimum value of. 
Q max / r w 3 can be taken from Figure 3.42 and 4> y from Figure 
3.4S. If ♦ <, 1," the vapor ‘flow is laminar and no further 

iterations are required. If, on the other, hand , -■* 1 , 

then the prevailing vapor flow regime .is turbulent and the 

’ r . ■ .* •- : 1 V 

procedure should be repeated using , the new value of 4> v to 
enter Figure 3.44. The optimum value of a can be found. 
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Figure 3.43 Optimum Value of a, Frank et al. [30] 








Figure 3.44 Optimum Value of Frank et al . 1.303' 

from Figure 3.43 once the iteration procedure has converged 
to constant values of ♦ . Frank claims that this converg- 
ence is rapid and may generally he accompl i shed by approxi- 
mately three iterations. Completion of this procedure 
finally yields the optimum values of Q max / , ' w ' the aspect 

ratio, a. This implies that Tor ii given heat pipe radius 
the maximum heat flow may be found, or converse ly , for a 
given heat flow, the minimum pipe radius may be determined 
and, moreover, for both situations, the optimum groove shape 
(as give. i by the aspect ratio). Only the capillary limit to 
the maximum heat flow has been taken into account. in the 
iteration procedure and other limits, such ns the entrainment 
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Figure 3.45 Graph for Determining 4 y , Frank et al. [30] 


limit which is especially important for open grooves, may 
well pose additional restrictions and warrants further 
consideration.- 

McKinney [35, 91] discussed a heat pipe having an adia- 
batic section. He found for fee r << 1 the pressure drop in 

the vapor to be identical to that reported by Haskin, equa- 
tion (5-69). The resulting expression for the maximum heat - 
transport rate was given as ’ 


.. r 2 2 

. ■ [2acosd „ " r w r v 

It : — p t glcosB “rr"2~ j "i ' 
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(3-98) 


where 
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F : or the case of Re r >> 1, McKinney employed Cotter’s results 
for the evaporator and condenser and assumed Poiscuille flow 
in the adiabatic section finding the heat flow 


Q 

X max 


t E,V «( 


P a lcos6 - 


2ccos0 


(3-100) 


where 


(l-4/ir ? ) 


(3-101) 


Sp 1 

^v h fg r v- * h fg M r w' r v] 


' (3-102) 


Notice that for Re >> 1; McKinney worked with the Poiscuille 

r . 

flow assumption for the vapor in the adiabatic pipe section 
whereas Mask in assumed fully developed turbulent flow for 
the evaluation of the pressure drop. For intermediate axial. 
Reynolds numbers, therefore, these two techniques may both 
be employed to encompass the actual pressure drop and the 
resultant heat transfer limit. McKinney with the aid of 

• t . * ' • ^ *’ 

equations (3-98) and (3-100) developed a computer program ‘ 
for the graphical disolay of v versus the ratio r/r - 
for various adiabatic lengths, pc rntcab i li tics , wetting ‘ 
angles, and temperatures. In all cases,; the maximum heat ; 

transfer tended to maximize in the neighborhood of r^/r w = ; 3. 
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The reason for the large deviation of his value from the 
values obtained by Cotter and Werner, i.e., /2/3 and 5/6 
respectively, cannot be explained without a more detailed 
study of the problem. 

In all of the preceding analyses the pressure drop 
in the vapor space in one way or another was considered. 

Since for low temperature heat pipes this pressure drop may 
generally bo assumed negligible, a quick estimate of the 
capillary limited heat. flux can be obtained by regarding 
only the pressure drop in the liquid. Several authors whose . 
•results arc discussed below have used this simplification. ' 
Phillips [29] and Phillips and Hindcrmann [S3] applied Darcy's 
law to the liquid flow and arrived at equation (3-25) where 1 is 
replaced by l c + 21 a + l c j from which he obtained for 1 = 0: 


Kgl KA w 

prj 


p i h f g 

v t 


KA W g cosB (3-1.03) 


In addition, he recommended the use of the following expres- 
sion for heat pipes which employ a bypass or arterial type 
of wick: . . \ 
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(3-104) 

The first term in equation (3-104) represents the pressure 
drop in the artery while the second, term pertains to the 
pressure drop associated with the flow of liquid to and from 
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the artery: in" the circumferential direction. It is note- : ;. : 
worthy that: the effective length utilized in the pressure ; " 
.drop calculation equation' (3-25) , is essentially V v - ' 

.l e /2 *♦ l a ♦ l c /2 and that l’ a * 0, which reduces , it to half 
of the total heat pipe length. Feldman (8) and St re clccrt 
and Chato (125, 126] have used the entire heat pipe length 
in their pressure drop calculations; consequently their : V. : 
resultant maximum, heat transfer rate is one half of that 
calculated with equation (5-103). A similar analysis has 
been performed by Neal' (28]. ? Langston and Kunz (13, 55] 
considered mass, momentum, and energy balances based on ah 
elemental thickness of wiick in the condenser. By assuming 
an infinite meniscus radius at the far end of the condenser 
and a minimum value at the condenser-evaporator interface, 
he obtained the limiting heat flux 


Q = h K 
X max max 


2|> t h f g ° A V c 
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cosB 


( 3 - 105 ) 


where the minimum meniscus radius has been evaluated from 




Y - - 

min : p g h 


( 3 - 106 ). 
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The subscript, KR, refers to the temperature conditions at • 
which the maximum Kicking rise is. measured. Implicit in 
the derivation of equation (3-105) is the assumption that 


■ ■> r-v>*s- r . 






e * 1. Cosgrove ct al. [37, 33] extended the analysis of 
Langston and Kunz to include the effect caused by an adia- 
batic heat pipe' section and the temperature . variation of 
the fluid properties. The maximum heat transfer is then 
given by . . . 




(3-107) 


where T a is the temperature of the adiabatic section and is 
equal to -the saturation temperature of the fluid, and 7 C is 
the condenser temperature and is taken as the average of the 
saturation and sink temperatures. Notice that for the 
special case of 1 = 0 equation (3-107) reduces to equation 

(3-105) . 

The above discussion has dealt solely with the capil- 
lary limit to heat pipe operation. As mentioned earlier 
cnii type of limit is especially important for low tempera- 
ture applications where relatively low vapor velocities and 
heat fluxes prevail. For high vapor velocities, on the 
other hand, the sonic and entrainment limits become import- 
ant. Levy [127] performed a one-dimensional compressible 
vapor flow analysis on a control volume basis restricted 
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to. the vapor space. Two models were used to relate the 


thermodynamic properties in the vapor. First byftreating 
the vapor as a perfect gas, the sonic limiting heat transfer 
rate was found to be vi 
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(3-108) 



This condition is reached when the vapor flow chokes at the 
downstream side of the evaporator. The second model 
described a ' single-component equilibrium two-phase saturated 
vapor and the analysis to which it was applied yielded a 
complex tranccndcntal equation for the limiting heat trans- 
fer rate. Equation (3-108) therefore nay he used to obtain 
the theoretical limiting curves which were illustrated 
earlier in Figure 3.17. Levy compared the limiting heat 
transfer rates predicted with both models for a particular 
sodium heat pipe. These limiting rates are illustrated in 
Figure 3.46 as a function. of temperature. Curve A was . 
obtained with equation (3-108) while curves B and C repre- 
sent the two phase model solutions. Curve B was calculated 
using the temperature at the upstream end of the evaporator 
and curve C was obtained using the temperature at the down- 
stream end. Also displayed are the wicking limit derived 
from Cotter's fundamental equations . and experimental data \ 
provided by Kcmmc'(SS). Relatively good agreement- i s dis-- • 
ccrnible between the sonic limiting curves and the experi- 
mental data for temperature less than 600°C. Above' that . 
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Figure 3.46 Comparison of Perfect Gas Model and Two Phase 
Model for Sonic Limit to Experimental Data, Levy [127] 

temperature the sonic limit curve greatly overprcdicts the 
measured maximum heat transfer rates and, in fact, the 
measured rates were probably limited by the pumping ability 
of the capillary system. It becomes also apparent from the 
agreement between curves A and B that the perfect gas model 
(equation (3-108)) is useful to estimate the heat transfer 
rate required to achieve choking within the evaporator 
section. • 1 / 

The theory of the entrainment limit has received little 
attention, the reason being the dependency of this limit on 


the detailsof the geometry and the. interfacial shear stress 
distribution. Cotter [19] and Kemme [86] claim that, for 
Weber numbers greater than unity, the possibility of entrain 
ment exists. The Weber number is the ratio of inertial 
force to surface tension force written. as 

p V . 2 £' - ' - . ‘ 

< Weber Number = ■ — ^ (5-109) 

where 1' is a characteristic dimension associated with the 
wick surface 1 . An estimate of the entrainment limited heat 
flux may be established by equating the Weber number to 
unity. This assumption, together with the energy equation 
[equation (5-19)), yields • -v ’ ; ;;V 

Kemme claimed that for screen wicks, the characteristic 
length, l * , is very nearly equal to the screen wire diameter 
and that it probably depends to some extent on the wire 
spacing. ' 

The heat transfer limit associated with boiling within 
the wick did not receive much attention cither.. This type 
of limit is difficult to predict since it requires, among 
other properties, a thorough- knowledge of the cavity dimen- 
sions in the wick and of the effective thermal conductivity 
of the saturated wick. The .boiling limit was illustrated ■ 
in Figure 2.7 and has to be considered as Qualitative in 




nature according, to Dcvcrall [18]. Notice that the limiting 
heat flux decreases with increasing temperature. Neal [28]- 
took the superheat which is necessary for the incipience of 
nucleation in the Kicking into consideration and related the 
superheat to the temperature difference existing across the 
wick in the' evaporator and obtained for the boiling limited 
heat flux: • ‘ , 


2 r X e k w q T sat f_2 cos0 c ^ 
ln r K /r v L r b 2,1 E r v. 


(3- 111) 


Marcus [128] contributed yet another relation for the boil- 
ing heat transfer limit . 
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(3-112) 


He further recommended to evaluate the effective wick con- 

l " ‘ 

ductivity with - 


k w ■ II . 


(3-113) 


Obviously the above expressions arc quite different from the 
qualitative limit given by.Uevcra.il since, both Neal and 
Marcus conclude that the limiting heat transfer increases 
with increasing vapor 'temperature. Hence the boiling limit 
curve, depicted in Figure 2.7, should have a positive slope 
instead of a negative one." A great deal more' experimental 
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and theoretical effort must be expended before the boiling: ‘ 
limit can be treated with sufficient confidence. i,'. 

The .foregoing discussions pertained to the numerous 
predictions of the maximum heat transport capability of heat 
pipes in view of the wicking, sonic, entrainment, and boiling 
limit respectively. Only a few analyses have been undertaken 
on other aspects of heat pipe technology. 

Lyman 3nd Huang 1129] studied • two dimensional 
liquid flow and heat conduction within the w'ick near the 
condenser entrance. Assuming constant pressure and a ‘ 
constant rate of condensation in the condenser he computed' : 
the temperature distribution in the wick. The results of ; ; 
his analysis are displayed in Figures 3.47 and 3.48.. The 
numbers on the isotherms represent the temperature above 
the coolant temperature, T^, in units of ^Q/bk^j while the 
numbers on the adiabatic curves give the fraction of the 
heat flow, Q, which passes through the portion of the wick \ 
to the left of t-he curve. Figure 5.48 illustrates the 
dimensionless temperature distribution in the midplane and 


at the surface of the- wick.' For wick matrices of low con- 
ductivity large temperature gradients are possible in the- 
wick at the junction between the adiabatic section and the 
condenser. Such temperature gradients have been observed 
qualitatively . in several experiments with' low temperature, 
water heat pipes. Unfortunately in general , the thermo- 
couples were not placed sufficiently close to accurately . 








Figure 3.47 Isothermal and Adiabatic Curves in Condenser 

Wick, Lyman and Huang C 1 2 9 3 
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Figure 3.48 Temperature Distributions on Midplane and 
. .. Surface o£ Wick, Lyman and Huang. [129] 
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verify the steepness of the gradient. It should be crapha-.. ' 
^sized that this analytical solution is applicable only to 
condensers with a fixed temperature boundary condition, c.g., 
calorimeter cooling. It is also evident that major conden-'-.. 
sat ion and the associated heat flow into the wick occur very 
near the entrance to the condenser section. The. effective 


wick conductivity was calculated by assuming parallel heat ^ 
conduction through the liquid and wick material as described 
by equation (5- 1,13) . Corring and Churchill (150] and Nissan 
cf al. (131] have suggested various techniques for the 
measurement and the computation of the ‘effective thermal 
conductivity of other wick materials. . ■ 

Brcsslcr and l'.'yatt (132] solved the differential equa- 
tion for the velocity during the transient capillary rise of 
a liquid in grooves of various geometries.' The mean velocity 
is plotted in terms of dimensionless groups shown in figure 
3.49. , The constant C equals [2 csc 3 (S/2)- 2 esc 2 (B/2) ] , 

(2J, (8(tr - 2 ) / it 2 ] and T* is [ -esc 2 (B/2) ] , [2], (2], for tri- - 
angular, semicircular, and square cross section grooves. 



respectively. For a given groove geometry, fluid properties, 
and temperature di fference between the wall and liquid sur- 
face, [t w - l' s , the mean velocity v may be determined from - - f f 

Figure 5.49. The total heat flux at steady state can then. : \ 
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Figure 3.49 Calculated Mean Velocity in Grooves, 

Bressler and Wyatt [132] 

Dressier used the above method to investigate the effect 
which groove geometry has on the maximum heat transfer. He 
found that a vertex angle of approximately 30° led to the 
highest heat transfer rates among triangular grooves. . 
Furthermore, it was found that square grooves are charac- 
terized by- the highest heat transfer rates per unit groove 
width when all grooves have an optimized depth. The results 
were different, however, when all grooves were compared at 
the same depth. Many more comparisons of groove charac- 
teristics may be made depending on the specific requirements 
to be evaluated. 

In another study, Galowin and Barker [133] employed a 
two parameter, fourth order velocity profile in the, Karman - 
Pohlhauscn boundary layer integral method to determine 
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velocity and pressure fields in a two-dimensional heat pipe. 
The vapor was assumed to be incompressible while the injec- 
tion and suction velocities at the wick surface were con- 
sidered as small. For the case of 1 = 1^ = L and uniform 

c c 

injection and suction rates with a velocity V f the velocity 
and pressure distributions wcr~ found to be: 

For 0 < z < L: • * 
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(3-116) 


For L < z < 2L: 


v (z) •= ft] f? . ll 
m UJ 2 [Rj (- I.; 


(3-117) 


P( z) ■ HO) - fli ) 1 
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(3-118) 


Approximate solutions were also obtained for the case where 
the injection velocity obeys a ramp function. 

Miller r.nd Holm [134] considered the possibility of " 
using model heat pipes to predict the performance of differ- 
ent prototype heat pipes. .A material preservation scheme 
was employed in which the same working fluid and wick mate- 
rial *?as used in both the model and the prototype. This 
implies that the thermal conductivities of the wall and 
wick, the emi trance of the condenser surface, and the 





permeability of the wick have to be preserved; With a 
starred quantity representing the model to prototype ratio 
of that quantity, the modeling equations are 



" ( T v - T oJ* " . 

. . (3-119) 


• t£ - (q*) l/ V(i *) ,/2 . V 

(3-120) 

and 

. '■ ' q* ,** N* A* . V .••• 

(3-121) 


Experimental verification of the above equations showed that 
prototype thermal behavior could be predicted from the model 
behavior to within 10°F over a temperature range of 140 to 
330°F for a pair of. water heat pipes. Modeling equations 
for another scheme which preserves the heat flux from model 
to prototype were also presented but an experimental veri- 

* 

fication was not attempted. 






tered in a more conventional single component heat pipe. 

In addition to variations in velocity, pressure and tempera- 
ture, concentration changes may occur throughout both the 
vapor and liquid phases. If it is assumed that local equi- 
librium conditions are present between the liquid and vapor, 
then the addition of a second component, in effect, provides 
an extra degree of thermodynamic freedom in the pipe; thus 
both the pressure and temperature must be specified to com- 
pletely determine the state of the mixture at a given loca- 
tion. Additional modes of mass transport, such as diffusion 
must also be taken into consideration. 

The phenomena which occur in the vapor core of a two 
component heat pipe are described by a set of eight coupled, 



nonlinear, second order, partial differential equations., 


If the vapor behavior may be regarded as two dimensional, 
then these equations include the momentum equations in each 


coordinate direction, the energy equation, two mass flux 


equations, an equation of state, and two continuity equations. 
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The latter two equations are required because the total mass 
and the mass of each component must be conserved which leads 
to three continuity equations, however, only two are indepen- 
dent. The eight equations mentioned above represent the 
relationship between the eight unknowns: 


u(x, y) 
v(x, y) 
J’ix(x. y) 
jivCx, y) 


pi(x, y) 
p(x, y) 
T(x, y) 
P(x, y) 


The solution of these equations requires a complete specifi- 
cation of the physical properties as well as a sufficient 
number of boundary conditions to define the problem. A 
consideration of the equations indicates that ! twenty boundary 
conditions are necessary. The velocity boundary conditions 
are easily specified by applying the" no-slip condition at 
the end walls and wicks and by assuming some sort of injec- 
tion velocity distribution at the wick surface. This - - 
injection velocity distribution is of course coupled to 
the external conditions applied on the container wall; 
e.g. it depends on whether constant temperature or flux 
conditions are utilized for the ultimate source and sink.- 
The temperature gradients normal to the end walls may also 
be assumed to be zero if the end walls are treated as adia- 
batic surfaces. The commonly applied condition of symmetry 


about the pipe centerline is not applicable due to the 
presence of the body force term which tends to produce an 
asymmetric density and temperature distribution depending 
on the pipe’s orientation with respect to gravity. Moreover, 
the density boundary conditions at the wick surface are not 
uniquely determined. In fact, to establish only the ratio 
of component density to mixture density (or mass fraction) 
at the wick surface requiros a prior knowledge of both the 
local temperature and pressure as well as an experimentally 
determined equilibrium curve relating pressure, temperature, 
and mass fraction. Considering the above difficulties, it 
is not surprising that the equations have not been solved. 

The prediction of complete component separation as 
presented by Cotter [4] is based on the assumption that only 
axial variations of temperature and concentration occur. 

His prediction 'is based solely on the criteria of mass con- 
servation and a consideration of the equilibrium curve for 
the binary two phase mixture. The characteristics of binary 
equilibrium curves are discussed in great detail by 
Bosnjakovic and Blackshear [135]. A typical set of such 
curves for the binary system water-methanol is illustrated 
in Figure 4.1 using data from references [136, 137]. Each 
enclosed curve represents the temperature and concentration 
characteristics at a liquid vapor interface for a given 
pressure. For- a particular pressure, any point below the 
closed curve represents a liquid phase only, while any point 
above the closed curve indicates that only the vapor phase 


exists. For points which fall inside of the curves, however, 
both phases are present. The compositions in the liquid and 
vapor are now determined by the intersection of the prevailing 
temperature line with the lower and upper branches of the 
curve, respectively. This implies that whenever a two com- 
ponent, two phase mixture with an equilibrium diagram of the 
type illustrated in Figure 4.1 exists, the vapor phase is 
always richer in the more volatile component than the liquid 
phase. If this is the case, then Cotter argued as follows. 

The same total amount of mass must be moving in opposite 
directions in the vapor core and in the liquid saturated wick 
at any given axial position to satisfy continuity. However, 
if the vapor is richer in one component, then more of that 
component is entering the condenser end of the pipe than is 
leaving. The net transfer of the more volatile component to 
the condenser end would continue until the two components 
are completely separated. Since the pressure in the vapor 
core is nearly constant, a complete separation into pure 
components would manifest itself by significant temperature 
variations from one end of the pipe to the other. In fact, 
where each pure component exists, the temperature would tend 
toward the saturation temperature of that component as deter- 
mined by the prevailing pressure. The result would be a 
temperature distribution along the pipe that consists of 
two plateaus with some transition region between them. 

The one dimensional argument discussed above r.ppcars to 
be highly idealistic. Since the velocity and pressure, fields 




are obviously two dimensional (due to the injection and suc- 
tion at the wick interfaces), there would seem to be no 
logical reason to assume that the concentration distribution 
is also not two dimensional. If this is the case then the 
vapor must be richer in the more volatile component only in 
the immediate vicinity of the liquid vapor interface. The 
conservation of species requirement would then mean that at 
each axial position the following condition must be satis- 
fied. 

I . " , 

•••• f o (y)u (y)dy .« f cp (y)u (y)dy (4-14) 
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vapor wick 



An infinite number of component density and velocity distri- 
butions exists which satisfy the above equation. It is 
readily apparent that situations may exist where a partial 
separation of components is possible. The exact nature of 
the resulting component distribution, however, must be deter- 
mined either experimentally or from! the solution of the 
eight equations mentioned earlier. 

Since a solution to the equations is virtually impossible 
due to an insufficient number of boundary conditions, an 
experimental solution to the problem was instigated. The 
experimental apparatus employed and the results obtained 
with it in the investigation are described in subsequent chap- 
ters. .. V . ... 'V;, 




CHAPTER V: EXPERIMENTAL APPARATUS 

Introduction 

The experimental program consisted of three phases. 

1 

During the early stages of the investigation it was decided 
that some effort should be committed to the study of Kicking 
materials. It was fclti that an adequate knowledge of vari- 
ous types of wicks was a prerequisite to the design of heat 
pipes. Of special interest was the maximum height to which 
liquids could be raised by capillary forces. Of equal 
importance was the magnitude of the velocity at which fluid 
could be delivered at any chosen height in the wick. This, 
velocity can then be used to calculate a conservative limit- 
ing mass circulation rate - assuming that all of the liquid 
is removed from the wick at a particular height. The pro- 
duct of this mass circulation rate and the latent heat of 

i - * 

vaporization provides a measure for the maximum heat trans- 
port capability of a heat pipe. The validity of using a 
dynamic velocity measurement technique to predict the steady 
state mass flow withii. a saturated capillary structure is 
somewhat questionable. Wetting difficulties associated with 
the moving liquid interface most certainly impede the flow. 
The difference between the rising and falling equi librium ■ 
heights discussed earlier in the literature survey was also 
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attributed to such wetting problems. However, a qualita- 
tive comparison of capillary induced mass flow rates :>ay 
certainly be made using this technique. . ^ ^ 

Hence an apparatus was designed and constructed which 
allowed in one experiment the measurement of both the maxi- 
mum wicking height and the liquid velocity as a function of 
height. Later a conventional cylindrical heat pipe was 
constructed and operated with a single component fluid. 

The construction of and experimentation with a co-planar 
heat pipe concluded the experimental phase of the investi- 
gation. The dimensions of this rectangular heat pipe were 
such that it should, however, be referred to as a vapor 
chamber. Two side walls were made transparent to allow 
visual and optical observations of the wick and the vapor 
core. Experiments in this heat pine were performed with 
both single and two component working fluids. 

IV i ck I nves t i gat i on 

For this phase of the investigation an apparatus was 
developed which allowed the advancing liquid front in a 
porous sample to be continuously monitored. The apparatus 
consisted of a rigid balance beam to which the test speci- 
men was attached at one end and compensating balance weights 
at the ether (see Figure S.l). The beam, was supported on a 
sharpened tool steel fulcrum which could be adjusted accord- 
ing to the weight of the test specimen. .. Once the dry speci- 
men was balanced with the compensating weights,* a Daytronic 
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(1S2A-1) load coll was rigidly attached adjacent to the 

compensating weights. The liquid level in the pan beneath 

the specimen was then slowly raised by adjusting the height 

of a larger liquid reservoir which was connected to the pan 

• . 1 

with a flexible syphon 'tube. This method of liquid height 
control eliminated any surface disturbances in the pan 
beneath the specimen and also guaranteed a constant liquid 
level as the specimen began to transport liquid from the. 
pan. The rising liquid exerted a force on the beam which 

■ - ii . , 

was sensed by the load cell and converted into an electri- 
cal signal which was suitably amplified and finally recorded 
with' a Honeywell strip chart recorder. Kith known porosi- 
ties of the test specimens and measured distances between 
fulcrum and the supporting point of the specimen and between 
fulcrum and the attachment of the load cell respectively, 
the strip chart recording, could easily he converted to a 
transient liquid height versus time graph. The maximum 
wicking height was simply the height at which no further 
variation of the signal could be detected. J * ' 

Eleven porous metal specimens acquired from the General 
Electric Company were tested. The porous metals are marketed 
under the trade name of "Foam Metal". They had various ,-' ;, 

porosities and were fabricated from two different materials, 
nickel and copper. Of the eleven specimens, only the two 
nickel samples of the highest density yielded repeatable 
data. Tt is believed that the' copper specimens oxidized 
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and that their wetting properties decreased to such an 
extent that essentially no water could be "wicked**. The 
maximum wicking height obtained with one of the nickel 
samples was 16. S cm. 

Cylindrical ilent P j pc 

Next a conventional cylindrical heat pipe was con- 
structed. The heat pipe was composed of three sections 
bolted together resulting in an overall length of 71. j cm 
(see figure S.2). The two end sections were made of 2.54 
era O.D. brass tubing while the center section was made of 
2.51 cm O.D. stainless Steel tubing. The lengths of the 
evaporator, adiabatic, and condenser sections were respec- 
tively 2.. I cm, 20.3 cn, and 2<».7 cm. The sections were 
separated from each other by bake] ire rings which were 
scaled with sheared annealed copper gaskets. The bnkclitc 
rings as well as the low thermal conductivity stainless 
steel adiabatic section reduced the axial heat conduction 
through the container wails. The brass walls of the con- 
denser and the evaporator tended to reduce the radial 
temperature drop. 

The condenser was surrounded by a 5.08 cm O.D stain- 
less steel tube. Cooling water supplied from a constant 
head supply tank was pumped through the annulus . (see Figure 1 
5.3). Since the pump was a constant displacement pump, a 
bypass line to the condenser was required to regulate the 
water flow through the condenser annulus. By observing the 
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Fig re S.2 Cylindrical Heat Pipe 

flow meter in the main line and by adjusting the valves in 
both lines, any given flow rate from 8.7 liters per minute 
to 178 liters per minute could be easily maintained. 

Nine chromcl-alumcl thermocouples were placed firmly 
against the exterior of the container wall at various axial 
positions. Thermocouple potentials were measured with a 
potentiometer by selecting any desired thermocouple with an 
isothermal switch. The reference junction was immersed in 
a distilled water-ice bath. . 

Four Ogden hose clamp type band heaters were installed 
firmly around the evaporator. The power input was con- 
trolled with a Powcrstat and monitored with a wattmeter. 

The maximum power input was 1000 watts. 

Small access tubes were provided on both ends of the 
heat pipe. The tube at the condenser was utilized for 
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Figure 5.3 Schematic Diagram of 
Cylindrical Heat Pipe Facility 
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■fluid addition or removal, and if desired, for the evacua- 
tion of the pipe with a mechanical vacuum pump. In opera- 
tion It also served as a pressure tap connection to a 
differential pressure transducer. The access tube on the 
evaporator was connected to the other side of the differ- 
ential pressure transducer and also to an absolute pressure 
transducer, both manufactured by the Pace Kngineering Com- 
pany. The pressure signals were amplified with a Pace 
carrier demodulator, model CD- 10, and finally recorded with 
a four channel Sanborn recorder. Prior to operation the 
differential transducer was calibrated against different 
heads of water acting on one side of the diaphragm, while 
the other side remained exposed to the atmosphere. Pull 
scale deflection on the Sanborn recorder was obtained for 
pressure differentials ranging from *1, psi in the evaporator 
to ♦ ! psi in the condenser, based on repeated calibrations, 
the accuracy obtained was t .02 psi for any given pen posi- 
tion. The absolute pressure transducer was calibrated with 
a licisc pressure gauge and full scale pen deflection was 
obtained for absolute pressures ranging from 0 to 45 psia. 
The accuracy for any given pen deflection was t .4 psia- . 

The heat pipe was mountedon a pivot which allowed the 
pipe to he operated in any orientation with respect to 
gravity. Although two wicks were employed during the experi 
rocr.ts, for the majority of the experiments a wick consisting 
of four layers of 150 mesh stainless steel screen with a 








• . t 


i 

t 



• l • _ ' •• ' ’ .. ••• 17 * 

wire diameter of .0026 inches was used. The screen was hold 
firmly against the interior of the container wall by a coiled 
steel spring. Trier to assembly, the pipe interior, as well 
as tl.e screen wick and spring were thoroughly degreased by 
repeated rinses with trichlorcthy lone , ethanol * and dis- 
tilled water. In addition, the stainless steel r-crccn wes 
passivated with a dilute nitric acid solution. The other 
wick used consisted of a silicon oxide cloth which is mar- 
keted under the trade name of Rcfrasil by the Mitco Company. 
The particular forn of this cloth used was a braided tubular 
sleeving with a nominal inner diameter of 1.91 cm. The 
sleeving was retained in the heat pipe in an annulus between 
the container wall and a smaller diameter perforated teflon 
tube. 

The pipe was prinari ly operated in the horizontal or 
near horizontal orientation both with and without insula- 
tion. The insulation consisted of approximately 3 inches 
of vermiculitc which completely surrounded the entire pipe. 
This insulation material was chosen because it could be 
simply poured into a retaining structure which surrounded 
the heat pipe and was likewise, easily removable. 

The energy transported from the evaporator to the con- • 
denser was calculated in the following manner. The pipe was 
first evacuated and if desired, insulated. The power input - 
was selected, and once steady state operating conditions \ 
were achieved nine axial temperatures were recorded. This" 
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procedure was repeated for several power inputs. The heat 
pipe was divided into nine segments, each segment containing 
one of the thermocouples. The initial evacuation of the 

i . • ■ 

pipe insured that all energy added was either lost from the . 
surface of a particular segment, or conducted to the next 
segment, i.c. no energy was transferred through the vapor 
core. Energy balances were made for each segment and the 
. heat loss from each segment was computed. This loss and. 
the measured temperatures were then used to calculate- the 
thermal resistance between each segm' ='..*• and the surrounding 
environment. Since the pipe was rut. n: various power 
levels, and consequently various temperatures, these thermal, 
resistances were established as a function of temperature 
in each segment. For the actual heat pipe test, the tem- 
peratures were again inc-asured and again an energy balance 
was made, this time however including energy transfer within 
the vapor space. The heat loss from each segment was found 
by employing the temperature dependent thermal resistances 
obtained with the evacuated pipe. The heat transfer rate 
due to heat pipe action was considered to be the sum cf the 
energy transfers to the vapor core in the evaporator. For., 
the uninsulated case, the accuracy of the above method was- 
very good. During the determination of the external thermal 
resistances as a function of temperature, a^-minimum of v 
95.4$ of the energy input could be accounted for, while for 
the insulated cases, a minimum of 37.4$ could be accounted 



, The working fluid in all tests was distilled water. 

The pipe was charged initially with a measured quantity of . 
water after which it was placed in a vertical position with 
the evaporator below the condenser. The heaters were 
turned on and steam and air were ejected from the pipe at 
the top where the steam was condensed, collected and 
measured. When the collected condensate indicated that the 
desired amount of water remained in the pipe, a valve was 
shut and the pipe was ready for operation. The vacuum pump 
was employed to periodically withdraw noncondcnsablc gases 
which occasionally accumulated in the pipe. The experi- ' 
mental data obtained with this heat pipe will be presented 
and discussed in chapter VI. 
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Meat Pipe 


A rectangular geometry was employed so that two sides 
could be made of glass, allowing visual observations of the 
wick, and possible optical measurements of phenomena occur- 
ring in the vapor core. The dimensions of the vapor space 
were 35 cm from end to end, by 22.7 cm from window.to window,' 
by 14 cm from wick to wick (see F : igurc 5.4). The entire : 
structure was made of stainless steel with .655 cm plates 
used for the longer side walls and 2.54 cm plates employed 
for the ends. The plates were welded together and the 
sides for the mounting of the windows were machined and . 
ground flat and parallel. The windows were 1.43 cm thick 
and were held in. place with two brackets Which fastened to 




174 



a. Thermocouple Prebet 4 Aectongulor Brett Tube* 

b Strip Heotert • Oteet Plota 

e wick 





the pipe with three brass strips which were screwed against . 
the .635 cm stainless steel .plates. 

' j ' ' - ' ' ' ' 4 * ’ 

Packing glands were provided on the condenser wall for 
nine thermocouple probes. Seven were symmetrically placed 
on the centerline between the two windows and two were 
placed close, to the windows to check the two-dimensionality 
of the device. Two - types of probes were used. The upper- 
most and lowermost probes on the centerline were purchased 
from the Conax Corporation and were 1 mm in diameter. The 
other probes were made "in house" and were l.S diameter 
and 40 cm long, and were of the copper- cons tantan. type. 

The packing glands contained a teflon bushing which was 
compressed against the probe shaft to provide a leak proof 
seal and yet permitted movement of the probes to any desired 
location. The two smaller diameter probes were located in 
thin grooves cut into the underside of the wick material.. 

The thermocouple probes were calibrated against a 
mercury in glass thermometer which could be read to the 
nearest .2 °F. A thermostatically controlled bath pro- 
vided constant temperatures for calibration from 55 °F to. 

151 °F. For convenience, a least square fit straight line 
was fit through the calibration points and extended for all : 
temperatures measured during the experiments. The accuracy , 
of this technique within the calibration, range was i .5 ‘F.'.V 
For temperatures considerably ■•.-in ‘excess of 151 °F, the error 
was estimated to be somewhat greater due to the nonlinearity 




of the thermocouples; however, it was anticipated, and later 
justified, that such high temperatures occurred only under* 
dry out conditions in the wick and hence an accuracy of 
t 1 1/2 - 2 °F was adequate. V 

The wiring diagram. for the temperature measurement . 
system is illustrated in Figure S.5. Eleven thermocouple . 
junctions were used in the experiments. During the early 
tests two' dimensionality of the device was ascertained and 
subsequently only seven thermocouples were employed in the 
heat pipe, all located midway between the windows. The 
remaining thermocouples were used to measure the ambient 
air temperature, thermostat bath temperature, and the inlet 
and outlet temperatures of the condenser coolant. All 
thermocouple voltages were measured against a single refer- 
ence junction which was immersed in a distilled water-ice 
mixture. All wire connections were made in an isothermal 
enclosure, and solid copper wires housed in a shielded 
cable led to an isothermal switch and ultimately to a 
digital volt meter of the type Vidar 510 Integrating Digital. 
Voltmeter (DVM) . The inaccuracy of this instrument accord- 
ing to the Vidar Manual is "typically less than i .07! at 
full scale including the effects of non-linearity, 24*hour 
drift, 5 °C temperature change, 10! line voltage change 
and attenuator inaccuracy all added in rms fashion" The . 
drift problem was reduced by re-calibration of the DVM 
which was done immediately prior to any data acquisition. 










The condenser surface of the heat pipe (vapor chamber) 
was maintained at a constant 'temperature by circulating 
cooling water through a bank of two times eight thin walled 
brass tubes, shown in Figure 5.4. This particular design 
was selected because it allowed the condenser area to be ' 
easily controlled by merely interrupting the flow in any 
number of tubes. The water was supplied from a constant " 
temperature bath and the flow rate was measured with a flow 
meter as depicted in Figure S. 6. Between the flow meter, 
and the actual condenser surface the flow was divided into 
sixteen currents by two manifolds placed in series. The * 
first manifold divided the flow into two streams, one. for 
each side of the heat pipe, and the second manifold which 
was mounted on the heat pipe support subdivided the two 
flows into eight parts each. The opposite manifold arrange 
menu was used to recombine the flows after they had passed 
through the condenser tubes. The two thermocouple probes 
used to measure the condenser inlet and outlet temperatures 
were placed inside of the "one to two" manifolds. For the ■ 
purpose of clarity the Vonc to eight" manifolds as well as 
the sixteen individual condenser tubes arc not shown in 
Figure 5.6. The constant temperature bath which supplied 
the condenser coolant fluid was a Lauda NBSD Thermostat. 

The thermostat consisted of a reservoir, circulating pump, . 
and a temperature control unit. The reservoir, was filled 
with distilled water and connected to. the manifolds forming 
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a closed loop, through which the fluid was circulated hy a 
Duplex pump,. The reservoir temperature was maintained 
constant by bucking the heat losses with an adequate heat 
input. The heat loss was controlled either by a continuous 
flow of cold tap water or by an auxiliary Freon refrigera- 
tion unit while the heat was supplied by an electrical 
resistance heater installed inside the hath. - This type of 
temperature control avoids the fluctuations normally 
encountered when the resistance heater is periodically 
cycled to maintain a desired temperature. Such temperature 
fluctuations were easily sensed by the thermocouples in the 
manifolds and since the fluctuations were out of phase with 
one another (due to the different positions in the flow- 
path), an accurate determination of the inlet and outlet 
condenser temperatures was difficult when the on-off tempera- 
ture controlling mode was used. ! 

The evaporator was heated with ten Ogcicn strip resist- 
ance heaters. Five heaters were located on each side of 
the pipe as shown in Figure S.4. All heaters were connected 
in parallel and the power input was controlled and measured , 
by an AC variac and a wattmeter, respectively. Each heater 
had a max-mum power rating of ISO watts at 120 volts, hence 
limiting the maximum power input to the heat pipe to 1500 
watts. As with the condenser, this design was considered 
desirable since the evaporator area could be varied simply - 
by disconnecting the leads to various heaters. 
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The gage pressure in the test section was measured with 
a mercury manometer which was connected to the condenser for 
horizontal operation of the pipe and to the evaporator for 
vertical- operation. The pressure could be read within 
t .5 mm Hg. A mercury barometer located near the test-..; 
section indicated the prevailing atmospheric pressure so 
that the absolute* pressure in the test section could be 
easily determined. .,'-1. " 

Access’ tubes were installed in the evaporator end wall 
for the charging of fluid to the pipe and for a pipe . 

connection to a mechanical vacuum pump. The procedures 
used. to charge and operate the pipe were as follows. The 
pipe was first evacuated to remove as much air as possible. 

It was then placed at a moderate angle fion the horizontal 
with the evaporator on top. About 800 cm 1 of fluid was 
forced into the pipe from a buret by atmospheric pressure 
resulting in a pool of liquid at the bottom.. The top wick 
was saturated by rotating the condenser into a top position 
and applying power to the lower bank of resistance heaters 
while only the top bank of condenser tubes was cooled. 

This caused a preferred condensation on the upper wick and,.: 
resulted in total saturation of both wicks. .Enough fluid 
was added so that a slight amount of excess 1 iquid • rema i ned 
on the lower wick after both wicks had been saturated. 

This procedure was repeated occasionally during the experi - , 
ments and especially after a relatively, high heat flux test 
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for which a partial drying of the top wick was often experi- 
enced. The majority of cases investigated were steady 
state, although several tests were made for which transient 
temperatures were recorded in response to a step change in 
power input to the evaporator. -«'■■■.. ; T * - ...... 

For the steady state modes , the temperatures were 
measured at 3 cm intervals starting 1 cm from the evaporator 
end wall., The condenser inlet and outlet temperature, as 
well as the room and hath temperature were recorded every 
time the thermocouples in the pipe were read, and the 
average values were used in subsequent data reductions. 

The gage pressure in the test section and the atmospheric 
pressure were read midway through the temperature readings 
so that an average value could be used if these values 
changed. A complete set of data required approximately one 
hour of time due to the repeated movement of the probes by 
3 cm increments and the subsequent waiting period for the 
temperature readings to stabilize. .7... ... / • r 

The energy, transported by the heat pipe was determined 
calorimctrically. However, because of the presence of heat 
losses from both the pipe and the calorimeter tubing and 
the effect of heat conduction along the, walls, the following 
reasoning in ascertaining these losses was used, h’hen the 
coolant flow is at a, temperature different from. room tempera- 
ture, the outlet temperature, T ut , will be different, from . 

the inlet temperature, T. n , due to the heat exchange with 
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room air. As shown in Figure S.7, if the ambient temperature 
is less than the coolant temperature (given by Tj n ) pad the 
power input to the heaters is zero, then T out is less than 
T. „ and vice versa. In addition, if for a particular differ- 
cnce in coolant and ambient temperatures the power input is 
greater than zero, then the additional energy conducted 
along the heat pipe wall will further increase (T ^ -T. n ). 

The procedure adhered to for the calibration consisted 
of the evacuation of the pipe, and the subsequent measuring 
of the calorimeter temperature difference, (T out - T^ n )„, 
for various coolant tenperatures and heat additions. The 
coolant flow rate was maintained constant at 200 kgn/hr. 

The resultant curves which arc illustrated in Figure 5.7 
were used as reference curves in the experiments. That is, 
any increase in temperature difference (T - T. ) above 
these curves was considered to be due to the mass transfer 
and phase change process taking place inside of the heat 
pipe. For example, if the coolant temperature was 25 °C 
above room temperature ond the heat input to the evaporator 
was 30C watts, the calorimeter temperature difference , 

expected would be approximately -1.5 "C. However, if a • 
temperature difference of;»-r>’5 °C was measured during the 
experiment, then a difference of 1. °C would be used, in .* . 



the calorimetric, equation to calculate the energy ; transport'*' .. | 

in the heat pipe. . Since during normal, heat' pipe operation ' ? 

the presence of fluid in the evaporator wick caused a low.or 
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temperature there than in the absence of fluid (such as was 
the case during calibration) , the reference lines in Figure 
5.7 would probably be situated somewhat lower and conse- 
quently the actual heat pipe energy transfer was somewhat 
greater than predicted by this technique. 

Experimental data obtained with the rectangular heat 
pipe will be presented and discussed in chapter VI. 
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. CHAPTER VI: EXPERIMENTAL RESULTS 

Balance Beam 

T'.ic balance beam described in chapter V was used to 
measure the transient Kicking characteristics of eleven 
porous "Foam Metal" specimens. A particular test was ini- 
tiated by allowing the lower surface of the sample to touch 
a liquid pool which was slowly elevated beneath it. The 
sample weight, and hence the liquid height, was continuously 
monitored for approximately thirty minutes. The sample 
surface was left exposed to the environment so that a small 
amount of liquid evaporation undoubtedly took place. How- 
ever, for runs of relatively short duration this effect may 
be neglected. This, of course, is not the case for long 
duration runs; for example, up to 2 months have been reported 
in other studies where a saturated environment must be 

employed. ■ . ■ - ' ' 

„0f all the specimens tested, the best Kicking rate was 
obtained with a 403 dense nickel sample. The transient 
height versus time curve for this sample with water as test 
fluid is presented in Figure 6.1.' Also illustrated is the : 
curve for dense nickel "Feltmctal" as supplied by the 

manufacturer. The cirvcs illustrated in Figure 6.2 were 
obtained by numerically differentiating the height versus 
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time curves and replotting as a function of height. These 
curves are useful in comparing the relative merit. of various 
wicks for potential. heat pipe use. From these curves it is 
easy to see why the M Feltmetal u was ultimately chosen for 
use in the coplanar heat pipe. 

For all nickel samples tested, the data were repeatable 
over a time span of approximately one month. For the copper 
samples, however, the wicking performance decreased with 
increasing time. The transient rise of water in a 451 dense 

it. 

copper sample is illustrated in Figure 6.3. The higher 
curve was measured shortly after the sample was received 
from the manufacturer while the lower curve was measured 
twenty-five days later. The change in performance was attri- 
buted to oxidation of the wick which was readily apparent 
from the discoloration of the copper. 

The significance of the data obtained with this device 
lies not in the actual values measured since none of the - 
samples tested were used in subsequent heat pipe designs , 
but in the method used to measure them. A device has been 
designed and tested which continuously and accurately moni- 
tors an advancing liquid front within a porous sample. This 
is a distinct improvement oveT spaced visual sightings of •; 
the advancing liquid since often the liquid line in the wick 
is difficult to see. To overcome such difficulties organic ; ! - 
dyes are often added to the liquid, or litmus strips are 
taped to the sample ^surf ace . Such techniques often influ-:. 
ence the data to such an extent, as to make it almost - • '< 
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meaningless. The balance beam, moreover, may be easily ad- 
justed to different sensitivities by simply moving the ful- 
crum. For long duration runs the device nay be adapted to 
provide a saturated environment around the sample by enclosing 
the entire beam within a ploxiglas box. 



C yl indrical Heat Pipe 

In all of the tests with the cylindrical heat pipe 
distilled water was used as a working fluid. Prior to the 
actual heat pipe tests, the pipe was operated dry at several 
power inputs. This, in effect, calibrated the pipe so that 
for subsequent tests with a fluid in the pipe, the actual 
heat transport of the pipe could be determined using the 
method discussed in chapter V. With a stainless steel screen 
as wick, and water as a working fluid, several tests wore 
conducted at different power inputs with the pipe oriented 
horizontally (0 “ 0°), and insulated. The measured axial 
wall temporature distributions are presented in Figure 6.4 
as a function of power input. The power values in parenthe- 
ses indicate the power actually transferred by the heat pipe 
while the first value represents the power input. Tempera- - 
ture profiles for the insulated case were quite similar even 
though the losses from the surface of the pipe were much 
smaller. The rather abrupt temperature change at the con- •• 
denser entrance is probably due to the presence of a non- 
condensablc gas. Carnesale, et al.- [36] and Marto and 
Hosteller [801 also observed this behavior in their low 













tenporature water heat pipes. As with the other experimental 
investigations,, the thermocouples were not placed close 
enough together to verify Lyman and Huang's [129] theory 
which predicts largo temperature gradients in the wick at 
the junction between the adiabatic section and the condenser 
for wick matrices of low thermal conductivity. It is inter- 
esting to note that tho wall temperature behaves as expected 
for a pipe with a floating temperaturo source and a fixed 
temperature sink. For the tests depicted in Figure 6.4 the 
absolute pressure within the pipe varied from 11 to 28.5 
lbs/in 2 and closoly paralleled tho saturation curve as 
expected. The differential pressures indicated that the •< 
evaporator pressure was above the condenser pressure but all 
values were very nearly zero and no trend could be detected. 

In Figure 6.5 tho effect of pipe orientation at a 
constant power input is presented. The power actually trans- 
ferred as indicated by the valuos in parenthoses, however, 
was not constant. It is evident that the temperature in 
the evaporator begins to increase for angles from the hori- - 

zontal as small as 3 degrees. When the angle was increased' 
to 8 degrees (corresponding to approximately a 11.5 cm > v 

difference in elevation between the evaporator and condenser) 
the evaporator temperature began to approach that obtained 
for a completely dry pipe, hence indicating almost complete 
wick dryout. This apparent dryout was caused by the inability 


of the screen wick to supply fluid to the evaporator at c 
rate equal to that of the evaporation process.- 
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Figure 6.5 Axial Temperature Distribution 
as Function of Inclination 







The heat pipe was also operated with a Refrasil wick. 
Several tests were made and temperature profiles qualitative 
ly similar to those. of Figure 6.4 were obtained. Continued 
testing, however, indicated a significant decline in heat 
pipe performance as evidenced by the high temperatures 
moasured in the evaporator for e/en relatively low heat flux 
tests. Subsequent inspection of the wick indicated a high 
degree of deterioration. The fiber was semi-dissolved and 
when squeezed produced a gritty milky white fluid which 
seemed to contain many fine fibers and powder. It was con- 


cluded that Pefrasil cloth is nqf, a desirable heat pipe wick 

Coplanar Heat Pipe > 

The experimental program utilizing the coplanar heat 
pipe consisted of 67 tests. Either pure water or pure 
methanol or some mixture of water and methanol was used as 
a working fluid in all tests. This combination of fluids 
was chosen because 1) they were easily attainable. 2) they 
are completely miscible, 3) their vapor pressure character- 
istics are such that a large temperature difference could 
be expected if the components completely separated at con- 
stant pressure', and 4) their binary equilibrium curves do 
not exhibit an azeotropic point. Of the 67 tests, ail but 
6 utilizod a relatively cool condenser sink temperature of 
approximately 63°F. This temperature varied about 3 degrees 
in either direction depending on the ambient and tap water 
temperatures. The remaining 6 tests were conducted with 



elevated sink tenperaturos (as high as 145®F). For the 
tests with the "cool" sink temperature, 45 were conducted 
with "normal" heat pipe orientations, i.e., either horizon- 
tal or with the evaporator above the condenser. Five differ- 
ent mass fractions of methanol (0.0, 0.25,0.50, 0.75, 1.00) 
wero employed at each of three orientations, i.o., horizontal 
(0 - 0°), o » 45°, and vertical (e - 90®), for each of 5 
power inputs to the evaporator (50, 100, 250 watts). All 
of these experiments were conducted under steady stato 
operating conditions in the pipe. In addition to these tests, 
2 trials were conducted with the evaporator below the con- 
denser (i.e., as a reflux condenser), 10 trials were per- 
formed to determine the influence that stratification ha- on 
the steady state temperature profiles, and 4 tests 'ere con- 
ducted to study the transient behavior of the pipe after a 
step change in power input to the evaporator. The effect of 
stratification and the procedures adhered to for data reduc- 
tion is discussed in Appendix A. 

For a series of constant composition tests, the pipe 
was charged in the following manner. First it was evacuated 
and then approximately 800 cm 3 of the desired mixture were 
pressure fed into the pipe. A vacuum was again applied 
until vigorous boiling of tho excess liquid took place. 

This assured that saturation conditions were present in the 
pipe. The pipe was then sealed and the upper wick was 
saturated as explained in chapter V. The amount of fluid 
added allowed for complete saturation of both wicks as well 
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as some excess fluid which resided on the lower wick. Each 
of the 4S "normal" tests were begun by setting the power 
input and pipe orientation to the desired values. All 
thermocouple probes were positioned so that the probe tips 
were 1 cm away from the evaporator end wall (x/L ■ 1/35). 
Steady state conditions were considered present when the 
measured temperatures varied less than . 1°F per hour. The 
achievement of steady state usually took between 5 and 8 
hours after the operating conditions were changed. For a 
particular test, the 1 temperature field within the pipe was 
determined by measuring the output signal of each of the 5 
probes in the vapor space as well as of the 2 probes under 
the wicks at 3 cm intervals along the pipe axis. In addi- 
tion, the room, thermostatic bath, and calorimeter inlet and 
outlet temperatures were recorded for each axial position 
of the probes. The temperatures were made nondimensional 
by subtracting the calorimeter inlet temperature and divid- 
ing by the difference between the maximum temperature and 
the calorimeter inlet temperature. This procedure was 
performed twice, once Using the maximum temperature measured 
in the vapor space only, and once including the temperatures 
under the wicks. The resulting dimensionless profiles were 
plotted by a Calcomp Plotter as a function of the dimension- 
less axial distance, x/L. ;'v> 

Several .typical temperature profiles for 8 * 0°, 8 “ 45°, 
and 8 = 90° are illustrated in Figures 6.6, 6.7, 6.3 respec- 
tively. These curves have in common a methanol mass fraction 
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of SOI and a power input to the evaporator of 100 watts. 

The profiles for the horizontal case (Figure 6.5) indicate 
that the vapor temperature changed very little in th:* axial 
direction. In fact for the centerline probe, with -he 
exception of the extreme right hand point the temperature 
range was between 82.6°F and 83.3°F. The temperature 
difference between the top and bottom vapor probes averages 
approximately 7°F and is due to stratification. As expected, 
the wick temperatures are higher in the evaporator and lower 
in the condenser. In addition, the top wick temperature in 
the evaporator is hotter than the lower wick temperature. 

This is to be expected because any condensation on the end 
walls or the windows is gravity fed to the lower wick while 
the upper wick may not be totally saturated. In addition, 
the upper wick is in contact with a higher temperature 
vapor than the lower wick. For the case where the angle of 
inclination 0 = 45° (Figure 6.7), the same trends are evident 
with one exception; now the vapor temperatures exhibit a 
definite variation along the pipe axis. This phenomenon is 
agfdn caused to some extent by the occurrence of stratifi- 
cation. The curves illustrated in Figure 6.8 are for a ver- 
tical orientation of the heat pipe and again display some 
degree of stratification in the axial direction. The vapor 
temperatures at any particular cross-section are now nearly 
constant with slight deviations occurring in the evaporator, 
section. The wick temperature? follow the same, trend as 
observed for the other pipe orientations. The. highest curve 
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in the evaporator section indicates a semi-dryout condition 

l - . 

for one of the wicks. This particular curve was measured 
under the wick which was designated as the bottom wick in 
Figures 6.6 and 6.7. This asymmetric condition was noted 
several times and means that one wick operated more effec- 
tively than the ether. :. v . ‘ 

The data obtained for other mass fractions and power 
inputs wore qualitatively though not quantitatively similar. 

i ..... • . . . ’ 

The phenomenon of wick dryout was observed several times and 
was especially prevalent for high methanol content and high 
heat input tests. ' 

An evaluation of the measured heat transfer rates leads 
to several conclusions. The average heat transferred for 
all tests was greatest for horizontal pipe orientations, less 
for 45° inclinations, and the smallest for vertical orienta- 
tions. In addition, the pipe efficiencies (as determined 
by the ratio of the energy transferred to the energy input) 
exhibited the same trend as a function or orientation. 

These average efficiencies were 40%, 22%, and 15.3% for hor- 
izontal, 45°, and vertical orientations, respectively. For 
the horizontal case, no correlation appeared to exist 
between the pipe efficiencies (averaged over the three power 
inputs) and the mass fraction of methanol in the pipe. For 
the 45° case, however, the average efficiencies were greater 
for the pure component tests than for the two component • 
tests with the minimum efficiency occurring at a methanol - 
mass fraction of 0.5. The vertical testr. indicated the same 
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/ behavior as the 45° .tests with the exception of the pure 
: methanol experiments where dryout often took place. 

When the heat pipe was operated as a reflux condenser, 
i.e., with the evaporator located below the condenser, no 
significant temperature variations were measured at differ- 
ent positions in the vapor core. The temperatures for a 
particular location, however, fluctuated up to 1°F. These 
fluctuations were probably caused by the presence of strong 

convection currents. The temperature profiles under the 

wicks were qualitatively similar to those obtained for hori- 
zontal heat pipe operation. 

The temperature profiles measured in the heat pipe when 
an elevated sink temperature was employed were similar quali- 
tatively to those obtained for lower sink temperatures. The 
higher temperatures, however, contributed to wick dryout as 
was evidenced by the lower heat transfer rates which were 

measured. ’ • ' .. " 

The pressures for the single component tests were found 

to be the saturation pressures corresponding to the wick 
temperatures in the adiabatic section. For a majority of 
tests, moreover, the pressure was determined by the wick 
temperature at the location where the temperature profile 
in the vapor crossed the temperature profile in the wick, , 
i.e., where the vapor and wick temperatures were identical.- 
The basic question as to whether or not the components r 
in a two component heat pipe separate will now be considered. 
- A conclusive verification of separation (or the'^lack^of it). ; 





was attempted using a Mach-Zehnder interferometer. The com- 
bination of the measured index of refraction changes and tho 
temperature field would then have indicated concentration 
changes. This technique, however, was not successful due to 
the inhomogeneity of the vapor condensation on the optical 
windows in all tests which made the attainment of reasonable 
interference patterns impossible. A referral to the multi- 
component, multiphase literature, e.g. Bosnjakovic and 
Blackshcar Cl3S3, nevertheless indicates that due to the 
isobaric nature of the vapor core, any separation must mani- 
fest itself in temperature changes. Moreover, these temper- 
ature changes should be dependent primarly on the difference 
in component concentrations between locations and should 
not be evident for heat pipes employing only single component 
working fluids. The behavior of two component heat pipes 
can thus be ascertained qualitatively by comparing their 
characteristics with those obtained for pure component heat, 
pipes. , . 

As evident from the . raw data, the temperatures within 
the vapor space for both single and two component fluids, 
were dependent to a large degree on the vertical direction. 
This behavior is explained by the fact that the vapor 
velocities . (on the order of 2 cra/scc) were insufficient to ■ 
completely overcome buoyancy effects which tended to form a 
stable stratified vapor." The resulting - temperature distri- ' 
butions within the vapor are thus a result of both the ‘ - 
stratification effect and tho operational characteristics - • 


of the heat -pipe. The stratification effect is essentially, 

f , . ;■ • 

a pure conduction process and may be accounted for using 
the technique described in Appendix A. The dimensionless 
temperature profiles illustrated in Figure 6.9 were obtained 
by subtracting the stratification profiles from the dimen- 
sionless temperature distributions which were measured when 
the heat pipe employed pure water as a working fluid. The * 
resulting curves indicate the behavior of the tenperaturc 
field within the vapor as compared to the behavior experi- > ‘ 
enccd for pure stratification. Several trends are dis- 
cernible from the figure. The corrected profiles for the 
pure water tests are consistently negative in the evaporator 
and positive in the condenser. This means that the center- 
line temperatures in an operating single component heat pipe 
are cooler in the evaporator and warmer in the condenser 
than the temperatures which would be observed if only pure ~ 
stratif ication accounted for the entire temperature field. 

Such a behavior indicates that the axial temperature distri- 
bution within the vapor space is tending toward an isother- 
mal condition as is to be expected for a single component- 
heat pipe. The higher temperatures in the condenser may be 
attributed to. the forced flow of the hot vapor from the 
evaporator which tended . to heat the vapor in the condenser. 

A similar trend was observed for the pure methanol tests for; 
the orientation where 6 *> 45°, although the behavior was not . 
as clearly pronounced. For the vortical tests the tendency ... 
toward isothermal izat ion for methanol was very . slight, -- 
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especially in the evaporator where dryout often occurred. 

However, the behavior in tho condenser qualitatively matched 

. . . 

that obtained for pure water. The behavior of the corrected 
profiles for the probes nearer to the wick exhibited the same 
behavior as tho centerline probes for both puro water and 
pure methanol. . .. 

Tho conduct of the two component tests may now bo com- 
pared to that obtained for pure water. To account for the 
pressure differences between various tests, the temperatures 
arc made nondincr.sional by subtracting the saturation tem- 
perature of methanol and dividing by the difference in the 
saturation temperatures of water and methanol. The satura- 
tion temperatures are evaluated at the measured pressure for 
each particular test situation. 

Tho dimensionless temperature profiles for the two com- 
ponent tests are plotted in Figures 6.10 through 6.15. The 
curves in each figure are lor a common power input and pipe 
orientation. In addition, the profiles are normalized with 
respect to the corresponding profile which was obtained for 
pure water. This means that the temperature distribution 
for pure water is represented in these figures by a hori- 
zontal straight line where the ordinate has a value of one. 
This normalization technique is used so that the behavior 
of the two component temperature profiles with respect to 
those for pure water are more vividly illustrated. It is 
not meant to imply that the pure water heat pipe operated in' 
a completely isothermal manner. . ' ‘ 
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In regards to Figures 6.10 through 6.13, several obser- 
vations may be made. In all cases the temperature profiles 
indicate a consistent behavior. That is, the temperatures 
are always higher in the evaporator and lower in the conden- 
ser. Moreover, with the exception of Figure 6.11, where the 
profiles are very close and even overlap, the tests with the 
higher water content working fluids exhibit higher tempera- 
tures. .. - 4 •• •: _ . - • •• ; 

The data (not shown) for a power input of 250 watts 
showed somewhat the same trends. However, the occurrence of 
dryout in the evaporator wick produced some deviations from 
the above profile shapes. This, of course, was due to the • > 
fact that there was no heat pipe action occurring in the 
dtied out evaporator. 

The effect of cross axial temporature variations for 
3 =» 90° and 0 ■ 45° are illustrated in Figures 6.14 and 6.15 . 

respectively. These curves are plotted with'respect to the ' ; 
stratification profile as measured for each particular probe. 
Each figure presents three sets of curves corresponding to ’• ' J 
C3ch of the two component mixtures. The two prof iles in . . • 7 . ‘ 

each set compare the corrected temperature profile near the \‘V- : 
wick with that obtained at the pipe centerline. It is clear - r J 7 . 
from the figures that the centerline temperatures are higher - 
than the near wick temperatures, in the evaporator and the . 7 

difference becomes much smaller as the condenser section is . 
approached . . 
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Construction details arc given for a low cost apparatus 
for the measurement of wicking properties. ' It permits an 
accurate and continuous indication of the position of the 
advancing liquid front within a porous sample. 

For the cylindrical heat pipe the axial temperature 
profiles as a function of power input behave as expected for 
a heat pipe employing a fixed temperature sink and a float- 
ing temperature source. The large temperature drop measured 
at the condenser end of the pipe is largely unaccounted for. 
Since the same behavior was observed by other investigators 
who used a fixed temperature sink, it appears that this phe- 
nomenon may be attributable tc this particular mode of opera 
tion and not to. the commonly accepted presence of a stagnant 
layer of noncondcnsablc gas. The performance of the cylin- 
drical heat pipe was significantly reduced (evidenced by the 
increase in evaporator wall temperature) as the inclination 
from horizontal was increased. The data showed that the 
maximum Kicking height for the stainless steel screen, as »t 
was installed in the pipe, was in the neighborhood of 12 cm. 
The number of screen layers used and the firmness with which 
they arc pressed against, the container wall certainly influ- 
ences the wicking characteristics of the wick. 

The data obtained with the coplanar heat pipe operated 
at ;j fixed sink temperature and with floating source tempera 
ftircs lead to the following conclusions. The temperatures 
along the wick display considerable variation from the 







evaporator to the* condenser. This heat pipe worked more 
efficiently and transported more energy when operated in a 
' hori zon t a 1 , pos i t i on as compared to a situation with the 
evaporator above the condenser. In addition, the low vapor 
velocities permitted stratification to considerably influ- 
ence the temperature distribution within the pipe (vapor 

chamber) . ' ‘ - V'." : . r • 

In regards *-o the behavior of two component heat pipes. 




the data indicated that under no conditions did complete 
separation of the components take place as deduced from a 
one dimensional argument. The two component heat pipe 
revealed the same temperature characteristics as a single 
component heat pipe when in a horizontal orientation. .This 
behavior mav be due to some mixing which took place in the 
liquid accumulated on the surface of the lower wick. Thus 
the entire vapor space, in effect, was in equilibrium with 
a constant composition fluid. ' . 

For the tests where the evaporator was located above 
the condenser, however, the excess liquid resided only at 
the condenser end of the pipe. For these cases, the axial 
temperature distribut-ions obtained- wi th a two component 
working fluid differed, to some extent, from those obtained 
with a single component working fluid. Such behavior indl- 
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it can bo argued that this effect was small for two reasons. | 
First, once the condensation formed on. the windows it 
appeared to remain relatively stable so that the condenser 
wick became the location of preferred condensation and thus 
played a dominant role in establishing the operating charac- 
teristics of the pipe. Second, the temperature distribution 
within the vapor core showed no changes in the cross axial 
direction (normal to the ivindows) while changes were evi- 
dent normal to the wicks. * . •.?<. 

The merits of the heat pipe have been somewhat exag- 
gerated. This is especially true for low temperature appli- 
cations in a gravity field where the selection of the wick 
is of primary importance. Most of the wicks used are incom- 
pletely described in the literature. Often no cleaning 
procedure is reported even though this is as important as 
the wick select j-'h itself. The use of the better "grooved", 
wicks requires sp< ;ial machining techniques which are not 
readily available to all investigators. ~r . 

The isothermal nature of heat pipes which is so highly 
touted in the literature does not always occur. In fact, 
for fixed temperature condensers, it appears to be the 
exception rather than the rule. 

While the heat pipe may be useful for high temperature 
and space applications, it shquld not be regarded as the 
ultimate solution to many ordinary heat transfer problems 
where greater economy, performance, and reliability are 
possible with more conventional heat transfer devices. 




. Rccoirjaondat ions 

• The study of the heat pipe literature leads to several 
recommendations. A coordinated materials research program 
is essential. In addition to assessing the desirability of 
potential wicking materials, this program should also pro- 

• " • . t 

vide for the investigation of the effect which various • 
cleaning procedures and treatment techniques have on the 

' < 1 -■ . ’ ’ < . * ■ j 

wicking characteristics of potential heat pipe WicTt. 

The influence, which the external boundary conditions 

' ' -.1 ‘ ■ ■ ' 

(i.e., cither fixed or 'floating temperature) have on the 
operational, characteristics of heat pipes is unclear, 

r" • • 

experimental study in which a given heat pipe could sequen- 
tially be subjected to the various types of boundary condi- 


tions would be useful. 

Additional work is needed on the theoretical predic- 
tions of the operating limits of heat pipes. In particular, 
the boiling and entrainment limits need more work. 

Several recommendations are in order regarding an ex-. ; 
pcrimental extension of the work reported in this study. The 
behavior of the two component working fluid could be more' 
easily detected by employing a heat pipe whose geometry 
dictated higher vapor velocities. ' If such a geometry was 

used, the consequence of stratification on the resul t ing 

* ^ 

temperature profiles could be neglected. ", 

The use of the interferometer to measure the concentra- . 
tion distribution within the vapor space still appears feasi . 
b.lc. However, major design changes in the test section arc 
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necessary in order to avoid the condensation of. the vapor on 
the inside of the windows. A design in which the windows 
aro thermally insulated from both the ambient environment 
and the walls of the pipe should solve tho condensation 
problem. Thermal insulation from the ambient conditions 
may be possible by using composite windows consisting of 
two pieces of glass with an evacuated space between them. 
Isolation of the windows from the test section walls, however 
appears to be very involved. * "•*?■ 

The chemical analysis of the vapor with a gas chromato- 
graph is desirable to ascertain the eventual separation into 
two components. Such a technique however, would virtually 
necessitate a direct connection of the chromatograph to the 
test section to insure that the vapor sample is not contarai- - 
nated. An investigation of two component mixtures whose 
equilibrium curves exhibit an ateotropic point could be 
undertaken. In addition, the effect of the various types of 
boundary conditions on the characteristics of a two component 
heat pipe should be investigated. 
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APPENDIX A: STKAT I F I CAT I ON 

Stratification had a significant influence on the 
steady state temperature field within the vapor space. It 
was especially important for the coplanar heat pipe because 
cf the low vapor velocities which were prevalent in the 
vapor core. 

Since the evaporator was located above the condenser 
for 0 » 45° and 0 ® DO®, stable operating conditions existed. 
That is, the lower density, higher temperature vapor tended 
to accumulate at the top of the vapor space in the evaporator 
while the higher density, lower temperature vapor resided in 
the condenser. Because. heat addition took place at the top 
of the vaper space, no major convection currents occurred 
and the primary node of energy transport (for pure stratifi- 
cation) was effected by conduction through the vapor space. 
Consequently, for a situation where no heat pipe action .. 
exists the temperature distribution within the vapor core 
must satisfy Laplace's equation under the assumption of 
constant thermal conductivity of the vapor. In reality, 
however, it may vary with temperature up to approximately 
ISt from end to end. : . .. q 

The temperature distribution within. the vapor-space was 
determined by utilizing the analogy between thermal and 




) .• 




electrical systems. The vapor space was modeled on a sheet 
of Tclcdcltos paper. It was .assumed that the vapor at the 
wick surface in the condenser was at some temperature T , 
and that the vapor at the wick surface in the evaporator was 
at another temperature, Tjj. In the electrical analog, this 
implies that these surfaces were maintained at a constant 
voltage. All other boundaries of the vapor space were con- 
sidered to be adiabatic and thur- were open circuits in the 
electrical analog model. The resulting lines of constant 
voltage (and hence the analogous isotherms) were determined 
using a DC Null Detector and arc illustrated in Figure A.l. 

The dimensionless axial temperature distributions corre- 
sponding to each of the probe locations are presented in 
Figure A. 2. These distributions arc, of course, highly 
idealistic resulting from the assumed boundary conditions, 
in actuality, the temperatures at the wick surface ftr both 
the evaporator and condenser were not constant Also the 
end walls undoubtedly were not adiabatic. "y ' j .. 

To verify whether conduction was indeed the dominant 
factor in establishing the stratified temperature distribu- 
tion, ten tests were conducted. Fach test employed a wick 
in only the evaporator end cf the heat pipe. A sufficient . 
amount of fluid was introduced into the coplanar pipe so 
that some excess liquid always remained in the condenser 
section. The pipe was evacuated until the liquid pool • v 
boiled, thus assuring that saturation conditions were ! 
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Figure A.i Pure Conduction Isotherms in Vapor Space 




ini t iaJ ly present . The heaters were then cner^i- .d anil the 
resulting steady state temperature distributions in the 
vapor space were measured. Since the absence of a wick in 
the condenser section prevented liquid from reaching the 
evaporator, the wicks rapidly dried out and the resulting 
temperature distribution was determined solely by the 

occurrence of stratification. As expected, the dimension- 

• * ' * * 

less temperature distributions were essentially independent 
of composition, pressure, and temperature. Two such distri 
hutions arc presented in Figures A . 3 and A. 4 for the center- 
line probe and for orientations of 0 = 90® and 0 = -15°. 
Similar distributions were measured for the other probes. 

The similarity between the measured distributions and 
those obtained with the electrical analog suggests that con- 
duction does indeed establish the temperature field. The- 
different values near the ends of the pine are probably 
caused by the idealized boundary conditions which were . 
employed in the electrical analog apparatus: . 

The experimentally determined stratification profiles 
for each probe were used for subsequent data comparisons 
with the stratification profile. 









The temperature measurement with a relatively high 
thermal conductivity probe in a medium of very low con- 
ductivity naturally leads to some errors. To assess the 
influence that the room temperature, T, imposes on the 

measurement of the vapor temperature, T , the following 

i , ! ■ • ' 

technique was used. The temperature actually measured at 

the probe tip, T, was related to T Q and using. the 

standard expression Cor the temperature distribution along 

a finite rod which protrudes from a heat source with heat 

losses to a constant temperature vapor environment as given 

by Jacob [138|. This expression is given as 


J3 r fl ♦ P)e' mL 
°o 1 * Pc' 2 " 1 * 


(B-l) 


where 



(B2) ■ 




In the expressions , I. is the length of the probe from the 
outsiJc of the pipe, C is the probe circumference, A is the 
cross sectional area, h is the convective heat transfer 
coefficient, and K is the effective thermal conductivity of 
the probe. The dimensionless temperatures calculated with 
the above expressions arc shown in Figure R.l as a function 
of both the probe length and the convective heat transfer 
coefficient. It is readily apparent that large errors are 
encountered when h is small. 

A conservative estimate of h was arrived at by assuming 
that the probes were orientated horizontally with free con- 
vection as the only mode of heat transfer to the surrounding 
vapor. Assuming a temperature difference of S C F in the' 
Grashof number and using an accepted correlation for the 
relationship of the Nusselt and Grashof numbers for free 
convection from a horizontal cylinder as given by Holman 
[1391, the convective heat transfer coefficient was found 
to be approximately 2.S Btu/hr ft °F. This value of h was 
used to calculate the actual vapor temperatures for all 
probe positions and all tests which were compared with the 
temperatures measured at the probe tip. It was found that 

the probe tip temperature was an accurate indication of the 

‘ * * * w * * * * 
vapor temperature with the exception of the measurements 

taken at the last position in the condenser. Kith the 

ex'ccption of this last point the readings were in general 

accurate to within ]°F. The measurements; taken at the 
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extreme probe position were often in error by ns much ns 
10°F. However, this occurred only for the high heat flux 
tests where the presence of forced convection certainly 
resulted in an increased value of h and thus smaller errors 
than are indicated by this approa€?». j 
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APPENDIX C: RECENT LITERATURE 


The articles discussed in this appendix were presented 
at the Space Technology and Heat Transfer Conference, Los 
Angeles, California, June 21-24, 1970. They arc presented 
as a supplement to the previously discussed literature so 
that the present report is as up to date as possible. 

Dcvcral 3 . 1 140 } has found that the construction of mer- 
cury heat pipes Tor high heat transfer rates is feasible, 
for operation between 200 and 360 9 C. Previously encountered 
wetting difficulties with mjrcury were virtually eliminated 
by the additions of magnesium and titanium. 

Schwartz [141] tested an ammonia-stainless steel heat 
pipe. The operating characteristics were compared to the 
characteristics obtained with a geometrically identical heat 
pipe employing water as a working fluid [58]. It was lound 
that the ammonia pipe was more efficient in transporting 


thermal loads than the water pipe up to an operating tempera 

* 

turc of approximately 90°F. Above this temperature the 
ammonia pipes relative advantage declined rapidly until dry- 
out, at which point the- water pipe was able to transport .< 
30 percent more energy. :/ - 


In a similar study. Waters and Kino [142] tested the' V 
ability of an ammonia heat pipe to-operatc for extended 


periods without failure by fluid loss or by degradation of 
, 1 . 
the energy transport mechanisms. The heat pipe used had an 

aluminum container and a stainless steel screen Kicking 
structure. Accelerated tine testing for both continuous 
heat pipe operation, and with alternating freeze-thaw cycles 
indicated no degradation in thermal performance. Subsequent 
metallurgical examination of the pipe revealed little mate- 
rial corrosion. The authors concluded that such a heat pipe 
should have a useful operating life in excess of 20 years 
when operated at about 80®l ; . . 

Heat pines in the cryogenic temperature range have been 
theoretically considered by Joy [143]. liquations for the 
optimum pore size, optimum wick thickness ratio, and maximum 
heat transport were derived. The effect of gravity was found 
to play a major role and must be taken into consideration in 
cryogenic heat pipe design. The equations result . in the 
selection of oxygen and a channel wick fer an optimum cryo- f 1 -' 
genic heat pipe design in the temperature range from -77 to 
90° K . . 


Chi and Cygnarowicz [144] also presented a theoretical 
analysis of cryogenic heat pipes. The influence of liquid 
property variations was found to be significant. . The theory 
compared favorably with the experimental results of ilaskin . 


Ferrell and Johnson [145] obtained experimental results 
for both the heat transfer coefficient and the critical heat 




flux through saturated beds of monel and glass beads. The 
liquid in the neighborhood of the heater was supplied only 
by capillary action. Various wick inclinations were employed 
in the tests. The conduction mechanism through a thin 
liquid-bead layer in contact with the heating surface as 
proposed by Ferrell and Allcavitch [41] was found to be 
substantially correct. ; 

Soliman, ct al . 1146) measured the effective thermal 
conductivity of both dry and water-saturated sintered fiber 
metal wicks. Correlations were found for the effective 
thermal conductivity in terms of the thermal conductivities 
of the solid and liquid phases and the wick porosity. Sub- 
stantial differences in the effective conductivity were 
found when measured either along or across the fibers and. 
this was attributed to the importance of the contact resist- 
ance between fibers. 

The effect of the working fluid, in cither the liquid 



or vapor phase, within the reservoir of hot reservoir gas- 
controlled heat pipes was investigated by Marcus and 
Flcischman (147). It was found that the presence of liquid 
in the reservoir at startup results in temporary pressures 
and temperatures in excess of design conditions. A perfor- 
ated non-wetting plug at the reservoir entrance. was found 
to eliminate this problem. ^ - . V 

Bliss, ct al ..[ 148) tested a flexible heat pipe subject 
to varying degrees of bend and various transverse and 
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Carlson and Hoffman I l 50 { studied the influence of mag- 


netic fields on heat pipes. Such effects were fount! to he 
important when electrically conducting working fluids (such 
as liquid metals) arc used and the pipe axis is not aligned 
with the magnetic field. For such cases, the presence of a 
magnetic field always reduces the heat flux capability of 
the heat pipe. However, this outcome may be reduced by 
designing the pipe with a compound wick structure with a 
larger liquid flow passage and a proportionately smaller 
vapor flow passage than for the no-rnagnctic field design. 
The equations necessary for such a design arc presented. 




